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Because the interaction is well-known, Coulomb excitation is one of the best tools for the investi- 
gation of nuclear properties. In the last 3 decades new reaction theories for Coulomb excitation have 
been developed such as: (a) relativistic Coulomb excitation, (b) Coulomb excitation at intermediate 
energies, and (c) multistep coupling in the continuum. These developments are timely with the 
advent of rare isotope facilities. Of special interest is the Coulomb excitation and dissociation of 
weakly-bound systems. I review the Coulomb excitation theory, from low to relativistic collision 
energies. Several applications of the theory to situations of interest in nuclear physics and nuclear 
astrophysics are discussed. 
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"I don 't know what they have to say, 

It makes no difference anyway - 

Whatever it is, I'm against it! " - Groucho Marx 
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I. WHAT IS COULOMB EXCITATION? 

Coulomb excitation is a process of inelastic scattering 
in which a charged particle transmits energy to the nu- 
cleus through the electromagnetic field. This process can 
happen at a much lower energy than the necessary for 
the particle to overcome the Coulomb barrier; the nu- 
clear force is, in this way, excluded in the process. 

Let V be the relative velocity of two nuclei at in- 
finity which determines the energy of relative motion 
E = where m is the reduced mass. The strength 

of Coulomb interaction can be measured by the Sommer- 
feld parameter 



hv 



(1) 



where Zi^2 are charges of the nuclei. For Z1Z2 > 137, or 
w <C c, the parameter ([T|) can easily reach values 77 S> 1 . 

This situation allows for the use of the semiclassical 
approximation: the Coulomb interaction is taken into 
account exactly in determining the classical Rutherford 
trajectory of relative motion R(i) where R is the distance 
between the centers of the colliding nuclei. Fig. [T] The 
relative energy E is assumed to be large enough so that 
we can neglect the feedback from the intrinsic excitations 
to relative motion. Then the trajectory is fixed by energy 
and impact parameter or defiection angle. 
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FIG. 1: Coulomb excitation occurs when a charged projectile 
passes by a target nucleus along a Rutherford trajectory. The 
coordinates used in text are shown. 



The classical distance of closest approach, 

2ZiZ2e^ 



a = 2an = 



(2) 



is larger than Ri 
Coulomb barrier 



i?2 at relative energy lower than the 



Ef 



ZiZ2e^ 
Ri + R2 



(3) 



The excitation is generated by the time-dependent field 
and the probability of the process is determined by the 
presence in this field of Fourier harmonics with the exci- 
tation frequencies 



Ef - E, 



(4) 



If the motion is too slow, the field acts adiabatically, 
the intrinsic wave function is changing reversibly and the 
probability of excitation is low. The corresponding adi- 
abaticity parameter is the ratio ^ of the time scales for 
the Coulomb collision, ~ ciq/v, and for the nuclear exci- 
tation, ^ l/uj: 



c 



ujao 

V 



(5) 



When C > 1, the situation is adiabatic and transition 
probabilities are small. 

The simplest treatment that one can give to the prob- 
lem is a semi-classical calculation, where the incident par- 
ticle describes a well defined trajectory, which is a classic 
hyperbolic trajectory of Rutherford scattering (see fig- 
ure [T|) . It has been proven that this treatment is valid 
in almost all situations studied in Coulomb excitation at 
low energies [AW75]. For high energy collisions, because 
the nuclear interaction distorts the scattering waves ap- 
preciably, a quantum treatment might be necessary for 
some observables, e.g. angular distributions [BN93]. For- 
tunately, at high energies, one can use the eikonal approx- 
imation for the scattering waves, which simplifies enor- 
mously the calculations. 
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The fundamental review paper [Ald56] contains a great 
deal of information on the subject and even now does not 
look obsolete, 40 years after its publication. However, 
in the last decades collisions between relativistic nuclei, 
with energies Eiab > 50 MeV/nucleon, have become a 
main tool of investigation in nuclear physics, in particu- 
lar for the study of nuclei far from the stability. Many 
new aspects of Coulomb excitation theory, such as the 
inclusion of transitions in/into the continuum, have been 
developed in the last two decades. It is thus timely to 
discuss the theory of Coulomb excitation from low to rel- 
ativistic energies. These notes, far from being complete, 
reviews the main aspects of the theory of Coulomb exci- 
tation. As the reader will notice, very little is said about 
the nuclear interaction as I want to emphasize the role of 
the Coulomb interaction in the excitation process. Also, 
nuclear structure and nuclear excitation models are dis- 
cussed only schematically and focused mainly on collec- 
tive properties, such as the giant and pigmy resonances. 

The readers are encouraged to study the Appendix sec- 
tion, where many basic quantum mechanics tools are dis- 
cussed. These tools will be used throughout the text. 

Beforehand, I would like to thank the JUSTIPEN 
(Japan-USA Theory Institute for Physics with Exotic 
Nuclei) for the financial support. Special thanks to Taka- 
haru Otsuka, Takashi Nakatsukasa and Naoyuki Itagaki 
for the invitation to participate in this school and for 
hosting my visit to Japan. 

II. INTERACTION OF PHOTONS WITH 
MATTER 

A. Electrostatic multipoles 

Electromagnetic multipoles appear in classical field 
theory as a result of the multipole expansion of the fields 
created by a finite system of charges and currents. We 
start with the system of point-like classical particles with 
electric charges located at the points r^. 

The electrostatic potential of this system measured at 
the point r is given by the Coulomb law. 

The function 



|r — r'l ^ J.2 _|_ ,j,i2 _ 2rr' cos 7 

depends on the lengths r, r' of two vectors and the angle 
7 between them rather than on the angles of the vec- 
tors r and r' separately. If r ^ r', this function has no 
singularities and can be expressed with the aid of the 
expansion over the infinite set of Legendre polynomials 
with the coefficients depending on r and r', 

— i— = f;PL(cos7)/L(r,r'). (8) 



Using the notations r< and r> for the smaller and the 
greater r and r', one can show that the expansion ([s]) 
takes the form 

The applications of the multipole expansion usually 
consider the potential ^ outside the system, i.e. at the 
point r with r > Vi. Then we can use the expansion ^ 
and the addition theorem 

^l(cos7) = E ^LM(n)r£M(n'), (10) 

M 

where n (n') is the direction of r (r'). With that we get 
^(^) = E ^^^^iM(n)Al(i?L,M). (11) 

LM 

Here the electric multipole moment of rank L, L ^ 

0. 1. 1s defined for a system of point-like charges i — 

1, 2, A as a set of {2L + 1) quantities 

A^(i;i,M) =^e,rfriM(n.), M ^ -L, -L+1, +L, 

(12) 

where the sum runs over all charges located at — 
{ri,9i,ipi) = (r^jiii). Exactly in the same way one can 
define, instead of the charge distribution, multipole mo- 
ments for any other property of the particles, for example 
for the mass distribution, Ci pi. 

In quantum theory, multipole moments are to be con- 
sidered as operators acting on the variables of the par- 
ticles. Containing explicitly the spherical functions, the 
operator A4{EL, AI) has the necessary features of the 
tensor operator of rank L. Introducing the charge den- 
sity operator 

p{r)=Y,e,5{r-r,), (13) 

i 

we come to a more general form of the multipole moment, 

MiEL,M)= J d\p{r)r'^YLMin), n=^. (14) 

In this form we even do not need to make an assumption 
of existence of point-like constituents in the system; for 
example, in the nucleus charged pions and other medi- 
ators of nuclear forces are included here along with the 
nucleons if p(r) is the total operator of electric charge 
density. As expected, we can separate the geometry of 
multipole operators from their dynamical origin. 

The lowest multipole moment L = is the monopole 
one. It determines the scalar part, the total electric 
charge Ze, 

A^(£;0,0)--^^e, = -^ f dh-p{r)^^Ze. 

(15) 
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The next term, L = I, defines the vector of the dipole 
moment 



d = 'Y^eiri = J(frp{r)r. 



(16) 



Taking into account the relation between the vectors and 
the spherical functions of rank L = 1, we obtain 



MiEl,M) 




(17) 



Subsequent terms of the multipole expansion deter- 
mine the quadrupole (L = 2), octupole {L = 3), hex- 
adecapole (L — A), and higher moments. 

In a similar way one can define magnetic multipoles 
A4{ML, M) related to the distribution of currents. The 
convection current due to orbital motion and the magne- 
tization current generated by the spin magnetic moments 
determine corresponding contributions to the magnetic 
multipole moment of rank L, 



M{ML,M) 



L+V 



(18) 

Here and 1^ stand for the spin and orbital angular 
momentum of a particle i, respectively; and gf are 
corresponding gyromagnetic ratios. (In this section, we 
measure all angular momenta in units of h and the gyro- 
magnetic ratios in the magnetons eh/(2mic)). 

The expression ( 18 1 vanishes for L ~ demonstrating 
the absence of magnetic monopoles. At L = 1, we come 
to the spherical components fiM of the magnetic moment 



M{Ml,M) 



"Mm, 



(19) 



(20) 



Higher terms determine magnetic quadrupole, L — 2, 
magnetic octupole, L = 3, and so on. 



B. Real photons 

1. Radiative decay 

The probability for the radiative (by emission of a pho- 
ton) transition for a nuclear transition i ^ f integrated 
over angles of the emitted photon and summed over its 
polarizations involves the same electromagnetic matrix 
elements as discussed above. We will not derive the equa- 
tions but only quote the results for the the probability 



for radiative decay [Ald56]: 
L[i2L + iy.\]^h 



LM 



A^(ML,M))^ ^} 



(21) 



where k — lo/c, and the multipole matrix elements for 
electric and magnetic transitions are expressed via the 
current operator j(r). 



M{ML,M) = / d?r{i ■ 1)$ 



(22) 



7W(EL,M) = J 



V-j -— 

K or 



K(r-j) $iM, (23) 



1 = — i[r X V] is the orbital momentum operator and the 
function ^lm arises from the partial wave expansion of 
the plane wave, 



LM 



.(2L + 1)!! 1 



1 



CK 



■ L 



jL{Kr)YLM{n) 



(24) 



where jLii^f) are spherical Bessel functions. Of course, 
for a given pair of states /, i and a given multipolarity L, 
only one term, either electric or magnetic, in (21 1 works 
if parity is conserved. 

An important practical case is connected to the long 
wavelength radiation: X ^ 1/k ^ R where i? is a size 
of the system. In nuclei R « XHA^I"^ fm so that the 
condition 



Ki? < 1 



(25) 



is equivalent io hio <^ 165 MeV which is usu- 
ally fulfilled. In the long wavelength approximation, we 
can use the limiting values of spherical Bessel functions 
3l{x) ~ /{2L -t- 1)!! at small arguments x. Then the 
magnetic multipoles ( 22 1 become 



A1(ML,M) 



c{L + l) 



dV(j(r).i)r^yiAf(n) (26) 



One can show that for the orbital current in Eq. ( 18 1 , this 



expression coincides with the orbital part of the static 
magnetic moment [EG88] . The presence of spin implies 
magnetization currents. In macroscopic electrodynamics 
such a current is c[V x m] where m stands for the magne- 
tization density. The analogous quantity for a quantum 
particle is 

m(r) = 5]gX'5(r-r.). (27) 

a 

The induced spin current is 

jspin = c^5:[VxsJ(5(r-ra). (28) 
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Being substituted into (26 1, this current reproduces the 



spin part of static magnetic multipoles as in (18 1 



In the long wav elen gth limit, the second term of the 
electric multipole (23 1 is smaller by a factor (kR)^ than 
the first one. In the first term V • j reduces to the time 
derivative of the charge density p'^^ with the aid of the 
continuity equation, 



V-j + 9/9"''/9i-0. 



(29) 



This time derivative can be expressed through the differ- 
ence of energies between the initial and final states gives 
huji^. Performing the expansion of the spherical Bessel 
functions we come to 



7W(EL,Af) 



LM 



(n) 



(30) 



which is a standard definition of the electric multipoles. 

Usually the angular momentum projection Mf of the 
final state is not measured. Then we have to sum the 
transition rate over all Mf defining the reduced transition 
probability 



B(EL;i^f)= J2 |(a^(EL,M) 



MM] 



(31) 



According to the Wigner-Eckart theorem [BM75], the en- 
tire dependence of the matrix element on the magnetic 
quantum numbers is concentrated in the Clebsch-Gordan 
coefficients of vector addition. Using instead the Wigner 
3j-symbols, we have for any tensor operator T^m 



{JfMf\TLM\W^ 



(-) 



Jf-Mf 



Jf 



L J,, 



-Mf M M, 



ifWTi 



L\W 



(32) 



where the reduced matrix element {f\\TL\\i) does not de- 
pend on projections. We can use the orthogonality of the 
vector coupling coefficients and perform the summation 
over M and Mf. 

The reduced transition probability is then related to 
the reduced matrix element of the multipole operator. 



B(EL;z^/) 



1 



2J, 



1 



if\\M{Em) 



(33) 



This quantity is convenient because it does not depend 
on the initial population of various projections Mi. Note 
that for the inverse transition induced by the same oper- 
ator, the detailed balance relation is valid. 



S(EL;/^z) 



2Ji + 1 
2J7TI 



B(EL;i^/). 



(34) 



The reduced transition probability determines the par- 
tial lifetime of a given initial state with respect to a spe- 
cific radiative decay (all Mf summed up). 



8n 



L+l 



Mf 



h L[(2i + 1)!!]2 



^B(EL; 



>f)- 
(35) 



Here the kinematic factors are singled out. They are 
associated with the geometry and phase space volume 
of the emitted photon. Information concerning structure 
of the radiating system is accumulated in the reduced 
transition probability. With the substitution EL^ML, 
the same expressions are valid for magnetic multipoles. 

Classically, the electromagnetic radiation emitted by a 
system is the result of the variation in time of the charge 
density or of the distribution of charge currents in the 
system. The energy is emitted in two types of multipole 
radiation: the electric and the magnetic. Each one of 
them is expressed as function of the corresponding mul- 
tipole moments, being the quantities which contain the 
variables (charge and current) of the system. If the wave- 
length of the emitted radiation is long in comparison to 
the dimensions of the system (which is valid for a 7-ray of 
< 10 MeV energy) the power emitted by each multipole 
is given by ([Ja75], chap. 16): 



Pe{LM) 



8tt{L+1)c 
L[(2i + 1)!!]2 



(-) \M{EL,M)\' (36) 



for electric multipole radiation and 
87r(i + l)c 



Pm{LM) 



(-j \MiML,M)\' 



i[(2L+l)!!]2 



(37) 

In quantum mechanics, the energy is not emitted con- 
tinually but in packets of energy hw. In a quantum calcu- 
lation the disintegration constant is the same as the num- 
ber of quanta emitted per unit of time when the power 
is given by the classical expressions (36 1 and (37 1. Thus, 



Xe{LM) 



Pe{LM) 

huj 

87r(L-f 



1) 



hL[{2L + 1) 



1112 



2L+1 



IMiELiM)]"^ 
(38) 



and 
Xm{LM) 



Pm{LM) 
87r(L- 



1) 



fiL[(2L-hl)!!]2 



2L+1 



\M{M, LM)\'^. 

(39) 



The decaying nucleus should also be treated as a quan- 
tum system. In this sense, the expressions for the mul- 
tipole moments M{EL, M) and M{ML, M) continue to 
have validity if we use for the charge and current densities 
the quantum expressions 



p/,(r)=vl/*(r)vl/,(r), 



j^(i') 



2im 



(40) 



(41) 



where the argument i (/) denotes the initial (final) state 
described by the wavefunction 'i'i (^'/). Equations (40 1 
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and ( 41 1 refer to a single nucleon with mass m that emits 
radiation in its passage from the state i to the state /. 
Thus, a sum over all the protons should be incorporated 
to the result, when we do the substitution of eqs. (40| 
and (|4l} in eqs. pSl) and ^ ([Ja75], chap. 16). 



Not only are the values of the magnetic multipole mo- 
ments are small in comparison to the electric moments of 
same order, but also the transition probability decreases 
quickly with increasing L, restricting the multipole or- 
ders that give significant contribution. To show this, it is 
sufficient to observe that the product {uj/c)'"Ai{EL,M) 
in pS] ) is at most equal to Ze{ujR/c)^, where R is the 
radius of the nucleus. For the energy that we consider, 
LoR/c is very small implying that, for the larger powers of 
L, the disintegration constant is also very small. These 
facts imply, in principle, that the electric dipole is always 
the dominant radiation. But the selection rules that we 
will see next can modify this situation. 



1 = ^ %.=RAr)/Vi7r 

FIG. 2: De-excitation of a proton to a level with L = 0. 



Let us take an example: if the initial state is a 3^ and 
the final a 2~, the possible values of L will be 1, 2, 3, 4 
and 5. But the change of parity restricts the transitions 
to El, M2, E3, M4 and E5, where El symbohzes an 
electric dipole transition {L = 1), etc. 



3. Estimate of the Distntegratton Constants 



2. Selection Rules 

The conservation of angular momentum and parity can 
prohibit certain 7 transitions between two states. The 
selection rules for the 7-radiation are easy to establish if 
we accept the fact that a quantum of radiation carries 
an angular momentum L of module L(L + 1) h and 
component z equal to Mh, where L is the multipolar 
order. Thus, in transition between an initial spin 1^ and 
a final spin 1/ the conservation of angular momentum 
imposes 1; = 1/ -I- L and, in this way, the possible values 
for the multipole order L should obey 



\h-If\<l<h 



(42) 



where = -^/Ji (/^ -t- 1) /i, etc. A special case is the 
transition 0+ ^ 0+: as multipole radiation of order zero, 
these transitions do not exist and they are effectively im- 
possible through the emission of a 7-ray. But in this case 
a process of internal conversion can happen, where the 
energy is released by the ejection of an atomic electron. 

Transitions between states of same parity can only be 
accomplished by electric multipole radiation of even num- 
ber or by magnetic radiation of odd number. The inverse 
is valid for transitions where there is parity change. Why 
this happens can be understood examining the definitions 
of the multipole moments. The functions that compose 
the integrand have definite parity and it is necessary that 
the integrand has even parity, otherwise the contribution 
in r cancels with the contribution in — r and the integral 
over the whole space vanishes. Let us look at the case 
of Eq. ( 30 ) : is always positive and the spherical har- 



monic 



i) is even if L is even. For Eq. (30 1 to be 
non-zero, should have the same parity of for L even 
and opposite parity for L odd. This justifies the transi- 
tion rule for the electric multipole radiation. A similar 
procedure applied to Eq. ( 26 1 justifies the selection rules 



for the magnetic multipole radiation. 



The use of Eq. ( 38 ) and ( 39 1 for the calculation of the 



transition probabilities in a real nucleus has the difficulty 
that wavefunctions that appear in eqs. (40 1 and (41 1 are 



not known. But, a prediction of their order of magnitude 
for the several modes can be done by assuming a single 
proton decaying from an excited state described by the 
wavefunction to a final state / with L = 0. as shown 
in figure |2] 

For an approximate calculation it is enough to do: 

Ri{r) — const. = Ri {r < R), 

R,{r) = {r>R) (43) 

and to use the same approach for Rf(r). The normaliza- 
tion yields immediately the values for the constants Ri 
and Rf : 



Ri — R 



f 



(44) 



where R is the nuclear radius. In this way, it is 
not difficult to calculate the electric multipole moment 
M{EL,M): 



M{EL,M) = ej 
which yields 

M{EL,M) 



3 Yl^{ 



i?3 



47r 



3ei?^ 
/4^(L-(-3)' 



-r^ dr dn, 
(45) 



(46) 



Thus, in this approximation, the disintegration constant 
in Eq. pSl) is 



2 .2 



E. 



2{L + 1) 
L[(2L + 1)!!]2 yL + 3j h \ hc 



2L+1 



R 



2L 



(47) 
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Number 
of cases 



-4 



3 - 

■exp/ -^c; 



logiri(A„p/A^^,J 



FIG. 3: Distribution of the ratio between the experimental 
and theoretical disintegration constants for transitions of the 
El and E2 type [Sk66]. 



III. LOW ENERGY COLLISIONS 

Coulomb excitation involves the same nuclear matrix 
elements as in radiative decay, but with different phase 
space factors. The reason is that Coulomb excitation is 
a process in which the excitation occurs when nuclei are 
outside their mutual charge distributions. Therefore, the 
Maxwell equation V • E = applies for the electric field 
inducing the transition. This means that only transverse 
photons fields, the same as for real photons, appear in 
Coulomb excitation [EG88]. In this section, we will prove 
this statement, and we will describe the semiclassical the- 
ory of Coulomb excitation for low energy collisions. 



Central collisions 



where we wrote explicitly the disintegration energy — 
hut. A similar calculation for the magnetic disintegration 
constant in Eq. ( 39 1 yields 



20(L+ 1) 
L[(2L+ 1)!!]2 \L + 3j h \hc ) 



2L + 1 



X R 



2L 



mcR 



(48) 



where m is the mass of a nucleon. For a typical nucleus 
of intermediate mass of ^ = 120, it is easy to see that 
^e/^m — 100, independently of the multipolar order 
L. Evidently, both constants decrease rapidly when the 
value of L increases. 

Disagreements of several orders of magnitude between 
the result of the calculation above and the corresponding 
experimental values can happen. In particular, if the ex- 
perimental disintegration rates are smaller than the ones 
predicted by eqs. (47 1 and (48), that can mean that Eq. 



According to Eq. (369) of the Appendix C, the proba- 



bility of exciting the nucleus to a state / above the ground 
state i is 



(49) 



with uj — (Ef — Ei)/h, is the probability amplitude that 
there will be a transition i f. The matrix element 



(50) 



contains a potential V of interaction of the incident par- 
ticle with the nucleus. The square of aif measures the 
transition probability from i to f and this probability 
should be integrated along the trajectory. 

A simple calculation can be done in the case of the ex- 
citation of the ground state J = of a deformed nucleus 
to an excited state with J = 2 as a result of a frontal 
collision with scattering angle of 6* = 180°. The pertur- 
bation V comes, in this case, from the interaction of the 
charge ZpC of the projectile with the quadrupole moment 
of the target nucleus. This quadrupole moment should 
work as an operator that acts between the initial and 



(liSl) is not very reasonable and that the small intercep- final states. The way of adapting ([11| is writing 



tion of the wavefunctions \E'i and \E' / decrease the values 
of A. Experimental values higher than predicted by eqs. 



V 



1 ZpC^Qif 



(47 1 and (48 1 can mean, on the other hand, that the tran- 



sition involves the participation of more than one nucleon 
or even a collective participation of the whole nucleus. 

Figures [3] and [4] illustrate the two situations: in figure 
|3]the experimental values of A for transitions of El multi- 
polarity are orders of magnitude smaller than that calcu- 
lated from Eq. (47 1. The opposite happens with the E2 



with 



(51) 



(52) 



where the sum extends to all protons. The excitation 
amplitude is then written as 



multipolarity where in most cases the experimental rate 
is larger than calculated; this is due to the fact that E2 
transitions are common among levels of collective bands, 
especially rotational bands in deformed nuclei. In fig- 
ure |4j on the other hand, what one notices is very good 
agreement between theoretical values and experimental 
ones for M4. This behavior is typical of transitions of 
high multipolarity. 



2ih 



dt. 



(53) 



At an scattering oi = 180° a relationship exists be- 
tween the separation r, the velocity w, the initial velocity 
vo and the distance of closest approach a = 2ao: 



dr 
dt 



( a\i/2 



(54) 
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an expression that is independent of the charge of the pro- 
jectile. It is, on the other hand, proportional to the mass 
of the projectile, indicating that heavy ions are more ef- 
fective for Coulomb excitation. 

The quadrupole moment operator Qif uses, as we saw, 
the wavefunctions and of the initial and final 
states. If those two wavefunctions are similar, as is the 
case of an excitation to the first level of a rotational 
band, the operator Qif can be replaced by the intrinsic 
quadrupole moment Q. The expression (59) translates. 



in this way, the possibility to evaluate the quadrupole 
moment from a measurement of a value of the cross sec- 
tion. 



FIG. 4: T(^haif-iife) > (R~ nucleaT radius, c = Ij as a 
function of the energy of the 7-ray for a series of transitions of 
the M4 type. One sees good agreement between the theoretical 
estimate (straight line) and the experimental data [GS51]. 



which is obtained easily from the conservation of energy. 
Besides, if the excitation energy is small, we can assume 



that the factor e*' 



in ( 53 1 does not vary much during 

Thus, 



the time that the projectile is close to the nucleus 
this factor can be placed outside the integral and it does 
not contribute to the cross section. One gets 



2if 



Zpe^Q 



if 



2ihvQ 



X 2 



dr 



r3(l -a/r)i/2- 



(55) 



The integral is solved easily by the substitution u = 1 — 
a/r, in what results 



(56) 



where the conservation of energy, E = 
ZpZTC^ /a was used, with mg being the reduced mass of 
the projectile-|-target system and Zt the atomic number 
of the target. The differential cross section is given by 
the product of the Rutherford differential cross section at 
180° and the excitation probability along the trajectory, 
measured by the square of a^/ : 



da 

dn 



dan 

dn 



|e=i8o° X loi/l^ 



(57) 



The Rutherford differential cross section is the classic 
expression 



and, at 



daji _ / ZpZT(? 
In ~ V 4:E 

180°, we obtain 



da 

dn 



e=i80° 



mgElQifl 

imz"^ 



(58) 



(59) 



B. Electric excitations 

The interaction hamiltonian responsible for the excita- 
tion processes in the nonrelativistic case can be written 
as 



H' 



d-\id'^r2 



pf(ri-Ri)pf (r2-R2) ^1^26^ 



ri -r2 



m 

(60) 

where charge densities of unperturbed nuclei depend 
on the distances from the corresponding centers. We 
subtracted the interaction between nuclei as a whole 
which determines the trajectory (time dependence of 
R = Ri — R2) but does not contribute to intrinsic exci- 
tations. We introduce the intrinsic coordinates for each 
nucleus Xi 9 = ri 9 - R1.2, 



d^xid^X2 



pf(xi)pf(x2) Z^Z^e" 
|R(i) + xi-X2| R{t) ' 



(61) 



and carry out the multipole expansion for a large distance 
between the centers, R 3> 2:1.2, 



H' = I d^Xid'x2pf'i^l)pti^2) 



An 



E 2L+ 1 flL+i(,) ^^M(-i^)^^V(R)-(62) 



L>0,Af 



Here X12 = xi - X2. 

This hamiltonian is rather complicated due to the cor- 
relations associated with the mutual excitation of the 
nuclei. It becomes much simpler if we are interested 
in the excitation of one of the partners only. Let the 
"projectile" 2 be not excited and we can neglect effects 
related to its structure as an extended object. Then 
P2^{x2) ~ Z2ed{:x.2), and the hamiltonian is expressed in 
terms of the electric multipole moments of the "target" 
h 



H' = Z2e 



in 



L>0,AI 



2L+1 R^+^t) 



. . ^. . . ^.^^^ 

The hamiltonian (63) is time-dependent since the trajec- 



tory is considered as a given function of time. According 
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to Eq. ( 369 1 the transition amplitude i ^ f with excita- 



tion by huj = Ef — Ei is the Fourier component for the 
transition frequency of the interaction hamiltonian taken 
along the unperturbed trajectory, 



aft 



(64) 



For the unpolarized initial nuclei and with no final po- 
larization registered, the transition rate is to be averaged 
over initial projections and summed over final projections 
of the target. 



2J, 



(65) 



MfM, 



the polarization state of the projectile is assumed to be 
unchanged. 

The trajectory enters the result via the time integral 



(co) 



dt 



1 



YLMmt))e 



(66) 



This Fourier component becomes small, cx exp(— const • 
C,) , if the trajectory changes at too slowly a rate compared 
to the needed transition frequency and the parameter 
of adiabaticity C > 1- We will discuss more about this 
integral later. 

The intrinsic matrix elements of multipole moments 



appear in the transition rate ( 65 1 in sums over magnetic 
quantum numbers 



^LM.L'M' 



2J, 



^ ^ M/,(EL,M)M*,(EL',M'). 



MfMi 



(67) 

The summation in Eq. (67) selects L' ~ L, M' = M and 



the result does not depend on M. 

J2 |M/.(EL,Af) 



* MfMi 

i3(EL; I ^ /) 
2L + 1 



^LL'^MM' 



SlL'^MM' 



(68) 



where i?(EL; i ~* f) are given by Eq. (33). 
The total excitation probability is therefore 



47rZ2e 



L>0 



i?(EL;»^/) 
(2L + 1)3 



Y.\lLM{u:f,)\\ 



M 



(69) 

From the viewpoint of the projectile the process is in- 
elastic scattering. The Coulomb trajectory defines the 
deflection angle 9 and the Rutherford cross section, Eq. 
(58). In our approximation, the trajectory is not influ- 



enced by the target excitation so that the inelastic cross 
section is factorized into the product of the Rutherford 



cross section ( 58 1 and the excitation probability ( 69 1 , 

dofi = dcTRWfi = ^ d(T/i(EL) (70) 

L>0 



where the cross section for the excitation of multipolarity 
L is equal to 



dn 



T:Z2ea \^B(EL;i 



/ 7rZ2ea \ 



hsm^{e/2)J (2i + l) 



M 

(71) 

where a is the distance of closest approach, Eq. ([2]). 

This theory can be extended, considering quantum 
scattering instead of classical trajectories, using relativis- 
tic kinematics, taking into account magnetic multipoles 
which become equally important for relativistic veloci- 
ties, and including higher order processes of sequential 
excitation of nuclear states. The last generalization is 
necessary for excitation of rotational bands and overtones 
of giant resonances (quantum states with several vibra- 
tional quanta). The mutual excitation of the projectile 
and of the target can also be studied. 



Estimates 



We can make a crude estimate of the cross section of 



Coulomb excitation. The trajectory integral (66), after 



changing the variable to dR — vdt, gives the dimensional 
factor /v. It has to be taken near the closest approach 
point ([2]) which is the most effective for excitation. The 
constants from the Rutherford cross section can be com- 
bined into the Coulomb parameter ([ij. As a result. 



cr(EL) 



2B(EL) , 



(72) 



where the function /l(C) depends smoothly on L but 
contains the exponential cut-off at large values of the 
adiabaticity parameter ^. We remember, Eq. (47 1, that 



in photoabsorption each consecutive multipole was sup- 
pressed by a factor (fci?)^. The situation for exciting 
higher multipoles is easier in the Coulomb excitation be- 
cause here 



(t(E,L+1) B(E,L + 1) 



a , 



(73) 



cr(EL) B(EL)a2 
This ratio is signiflcantly larger than (kR)^. 

D. Inclusion of magnetic interactions 

A more accurate description of Coulomb excitation 
for all scattering angles requires the correct treatment 
of magnetic interactions and the Coulomb recoil of the 
classical trajectories. Here we will discuss the role of 
magnetic interactions. Note that magnetic interactions 
induce electric transitions, too. And electric interactions 
also induce magnetic transitions. Thus, in electromag- 
netic excitation, both interactions mix and can only be 
isolated under special circumstances. We now show how 
magnetic interactions modify the results obtained in the 
last section. 
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Following the derivation of the electromagnetic inter- 
action presented above, the excitation a target nucleus 
from an initial rate | z > to a final state | / > is, to first 
order, given by 

a/. = - / dt e'(^/-^')*/'' < / I H„t I I > (74) 

where (see Appendix D) 

< / 1 n^nt \ i> = 



1 



Pfi<P{ic,t) - - j/,(r) • A(r,i) 



(75) 

with pfi{r) and j/i(r) given by Eqs. (40l and (41 1, re- 
spectively. 
Thus, 



dt 



(t)(r,t) d^r. 



(76) 



where v is the velocity of the projectile, and we used 



A(r,i) = ^0(r,t) 



(77) 



valid for a spinless projectile following a classical trajec- 
tory. 

The scalar and vector potentials at the target nucleus 
interior. 



r 

(w, r) = ZpC / 

J — c 



A{uj, r = 

c 



1 



■r'(t) 



dt 



r'it) 



■dt 



(78) 
(79) 



are generated by a projectile with charge Zp following 
a Coulomb trajectory, described by its time-dependent 
position r(t). 



We now use the expansions of Eqs. a9[ and (10 1 and 



assuming that the projectile does not penetrate the tar- 
get, we use r> (r<) for the projectile (target) coordinates. 

To perform the time-integrals, we need the time depen- 
dence, r'(t), for a particle moving along the Rutherford 
trajectory, which can be directly obtained by solving the 
equation of angular momentum conservation (see [Go02] ) 
for a given scattering angle in the center of mass system 
(see figure [ij. Introducing the parametrization 



rix) = flo [e coshx + 1] 



where 



ao 



ZpZre'^ 



one obtains [Go02] 

t - 



and 



1 



sin {d/2) 



ao 



[X + e sinh x] ■ 



(80) 



(81) 



(82) 



Using the scattering plane perpendicular to the z-axis, 
one finds that the corresponding components of r may 
be written as 

X = flo [cosh X + e] , 

y = flo - 1 sinh x, 

z = . (83) 

The impact parameter b in figure [l] is related to the scat- 
tering angle ?? by 



ao cot(z?/2) 



(84) 



and the eccentricity parameter e is related to it by means 
of 



62 

7^2"' 



(85) 



In the limit of straight-line motion e ~ b/ao S> 1, and 
the equations above reduce to the simple straight-line 
parametrization. 



y = vt, 



and z = 



(86) 



Using the continuity equation, V-j/i = —iujpfi, for the 
nuclear transition current, the terms in the expansion of 
the scalar and vector potentials mix up. After a long but 
straightforward calculation, one can show that the result 
can be expressed in terms of spherical tensors (see, e.g., 
Ref. [EG88], Vol. II) and Eq. ([76| becomes 



in 



ih ^ 2L+1 

LM 



(-1) 



M 



X \S{EL, M)Mfi{EL, -M) 
+ 5(Mi, M)Mf^{ML, ~M) 



(87) 



where 



Mf^{EL,M) = 



(2L+1)!! 



Ki+ic(L+ 1) 

yj/,(r).VxL[ji(«:r)yiM(f)]dV, 



X j ij,{v)-\.[jL{nr)YLM{i)]d\ . (89) 
The orbital integrals S{7tL, M) are given by 

L{2L-l)\\ d? 
xYLM[0'{t),(b'{t)] e'^'dt 



S{EL,M) = 



^-{r'{t)hL K(<)]} 



,L+2 poo 



cL{2L - 1) 
X YLM[e'{t)A\t)]e''^'dt 



„ , v'(i).r'(<)/iiK(<)] 



(90) 
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and 

S{ML, M) 



moc L{2L - 1)!! 

X YLM[0'{t),cf,'{t)]} 



e"^' dt 



v'{/iiK(i)] 

(91) 



where Lq is the angular momentum of relative motion, 
which is constant: 



Lq — flomof cot — . 

In non-relativistic collisions 

, ujr' V Lor' V 
nr = = <-<l 

C C V c 



(92) 



(93) 



because when the relative distance r' obeys the relations 
ujr' /v > 1 the interaction becomes adiabatic. Eq. (93 1 is 
long wavelength approximation, as discussed in connec- 
tion to Eq. ( 25 1 . Then one uses the limiting form of hj^ 



for small values of its argument [AS64] to show that 



S{EL, M) -- 
and 

S{ML, M) 



r'-''-\t)YLM{e'{t),^'{t)}e''^'dt 



(94) 



1 



Lm„c 



■dt 



(95) 



which are the usual orbital integrals in the non- 
relativistic Coulomb excitation theory with hyperbolic 
trajectories (see eqs. (II.A.43) of Ref . AW75]). 

It is convenient to perform a translation of the inte- 
grand by X — > X + i (7r/2) [AW75]. This renders many 
simplifications in the calculations of the orbital integrals 
S{ttL,M), which become 

S{EL,M) = ^ I{EL,M), 



S{ML, M) = - 



Leah 



[(2L+l)/(2L + 3)]i/2 



X [{L + rf - M2]1/2 cot {'d/2) I{ML, M) , 

(96) 



with 
Clm ~ 



2L+1 y/ {L - M)\{L + M)\ 

47r {L-M)\\{L + M)\\ 




(■_2^)(L+Af)/2 



(97) 



where the upper (lower) form is valid for L + M — even 
(odd) The (reduced) orbital integrals are given by 

/oo 
dx exp[-Cecoshx + <x] 
-C30 



(e -f i sinh x ~ ^/ — \ cosh x) 
(iesinhx+ 1)-^+*^ 
/(ML, M) = I{E, L + l,M) , 



M 



(98) 



where 



C = ujao/v. 



The square modulus of Eq. ( 87 1 gives the probability 



of exciting the target from the initial state | liMi) to the 
final state | IfMf) in a collision with the center of mass 
scattering angle i). If the orientation of the initial state 
is not specified, the cross section for exciting the nuclear 
state of spin // is 



dai^f = 



1 



4 2L 



- ^ |a/, I' dn, 

Mi-Mf 



(99) 



where a^e^d^l/A is the elastic (Rutherford) cross sec- 
tion. Using the Wigner-Eckart theorem, Eq. (32) and 



the orthogonality properties of the Clebsch-Gordan coef- 
ficients, one gets (for more details, see [BP99]) 



da,^f 47r2Z|,e2 



dn 



,2,4 ^ ^(y^^;^';;/^) |5(7rL,M)'^ 



{2L+lf 



(100) 

where tt = E or M stands for the electric or magnetic 
multipolarity, and B{TrL, Ij) are the reduced transition 
probabiUty of Eq. ( 33 ) . 



E. Virtual photon numbers 

1. Angular dependence 



Integration of (100) over all energy transfers E^ — fku, 



and summation over all possible final states of the pro- 
jectile nucleus leads to 



d(Jc 

~dn 



E 



da 



dn 



PfiEj) dEy , 



(101) 



where pf{E^) is the density of final states of the target 
with energy Ef = Ei + E^. Inserting Eq. (100) into Eq. 
( 101 1 one finds 



dn 



E 



da. 



dn 



(102) 

where cr^^ are the photonuclear cross sections for a given 
multipolarity ttL, given by 

L[{2L + iy.\y 



(103) 



The virtual photon numbers, UTrLiEj), are given by 

dn^L Zla L[(2L+l)!!f fa^^ 
^ = ^ (L+1)(2L+1)3 E I ^(-^' I 

(104) 

where k — E^/hc, and a — 1/137. 
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In terms of the orbital integrals I{EL,M), given by 
Eq. ( 98 1 , and using the Eq. ( 104 1 . we find for the electric 



niultipolarities 



dn 87r2 ^v' (L + 1)(2L+1)2 ^ 
{L-M)\{L + M)\ 



E 



[{L - M)\\{L + M)\^] 



I{EL, M) 



(105) 



In the case of magnetic excitations one obtains 

duML ^ ^ /C 2(L-i) [{2L+l)\f 

dfl 87r2 \' L(L + 1)(2L + 1)2 

X r'^^+^e^ (e2 _ 1) 

[{L + 1)2 - M2] (L + 1 - Af)!(L + 1 + Aiy. 



E 



L + M = odd 



[(i + l-Af)!!(L + l + M)!!] 



X I I{ML,M) p . 



(106) 



2. Impact parameter dependence 



Since the impact parameter is related to the scattering 
angle by Eqs. (84) and (85 1, we can also write 



4 dn^rL 



2TTbdb ale^ dVL 



(107) 



which are interpreted as the number of equivalent pho- 
tons of energy = hu, incident on the target per unit 
area, in a collision with impact parameter b. The impact 
parameter dependence of the Coulomb excitation cross 
section is 



dac _ \ ^ 

2Mb ~ ^ 

ttL 



dE. 
~E. 



^ n^L{E^,b)<j;'^{E^). (108) 
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The total cross section for Coulomb excitation is ob- 



tained by integrating Eq. (108 1 over b from a minimum 



impact parameter 6min, or equivalently, integrating Eq. 
( 102 1 over up to a maximum scattering angle i^max, i-e. 



^c = E/ ^N^L{E,)a;^{E,), 
where the total number of virtual photons is 

nOO 

N^l{Ej) ^ 2tt dbbn^L{Ey,b) 



(109) 



2tt 



desine- 



,dnT, 



dn 



(110) 



This condition is necessary for collisions at high en- 
ergies to avoid the situation in which the nuclear in- 
teraction with the target becomes important. At very 
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FIG. 5: Total number of virtual photons for the El multi- 
polarity, "as seen" by a projectile passing by a lead target 
at impact parameters bmin ~ 12.3 fm and larger (i.e., inte- 
grated over impact parameters), for three typical bombarding 
energies. 



low energies, below the Coulomb barrier, femin = and 

l9max = 1800. 

The concept of virtual photon numbers is very useful, 
specially in high energy collisions. In such collisions the 
momentum and the energy transfer due to the Coulomb 
interaction are related by Ap — AE/v ~ AE/c. This 
means that the virtual photons are almost real. One 
usually explores this fact to extract information about 
real photon processes from the reactions induced by rel- 
ativistic charges, and vice- versa. This is the basis of the 
Weizsdcker-Williams method [Fe24,WW34] (it should 
be called Fermi's method - see historical note later), used 
to calculate cross sections for Coulomb excitation, par- 
ticle production, Bremsstrahlung, etc., (see, e.g., Ref. 
[Ja75,BB88]). 



3. Virtual and real photons 



We have shown that even at low energies the cross sec- 
tions for Coulomb excitation can be written as a product 
of equivalent photon numbers and the cross sections in- 
duced by real photons. The reason for this is the assump- 
tion that Coulomb excitation is a process which involves 
only collisions for which the nuclear matter distributions 
do not overlap at any point of the classical trajectory. 
The excitation of the target nucleus thus occurs in a re- 
gion where the divergence of the electric field is zero, 
i.e. V • Ep (t) — 0, where Ep (t) is the electric field gen- 
erated by the projectile at the target's position. This 
condition implies that the electromagnetic fields involved 
in Coulomb excitation are exactly the same as those in- 
volved in the absorption of a real photon [EG88] . 
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4- Analytical expression for El excitations 



For the El multipolarity the orbital integrals, Eq. 



(98 1, assumes a particularly simple form. It can be per- 



formed analytically for M = 0,±1 [AW66]. Using these 
results, one gets the compact expression for the El vir- 
tual photon numbers 



driEi 

dn 



47r2 



(111) 



where KiQ is the modified Bessel function with imaginary 
index, 



r(zc+i)(2z)< 



dt- 



cos t 



(t2 4-z2)<+5 



(112) 



K'^^ is the derivative with respect to its argument. 

This result is not particular useful, as one still has to 
perform a time integration in the equation above. How- 
ever, as we will see later, the above formula will help us 
to understand the connection with relativistic Coulomb 
exctitation. 

In Figure [5] we show a calculation (with = — ex- 
citation energy) of the virtual photons for the El multi- 
polarity, "as seen" by a projectile passing by a lead target 
at impact parameters equal to and exceeding b = 12.3 fm, 
for three typical bombarding energies. As the projectile 
energy increases, more virtual photons of large energy 
are available for the reaction. This increases the number 
of states accessed in the excitation process. 



F. Higher order corrections and quantum 
scattering 

The results presented in this section arc only valid if 
the excitation is of first-order. For higher-order excita- 
tions one has to use the coupled-channels equations ( 366 ) 



of Appendix B, with the excitation amplitudes in each 
time interval given by Eq. ( 87 1 . Important cases of appli- 



cations of Coulomb excitation require a non-perturbative 
treatment of the collision process. We will discuss some 
of these cases in later sections. 

A full quantum calculation for Coulomb excitation uses 
scattering waves for the projectile (and target). The cross 
section in first order perturbation theory is given by Eq. 
(379 1 of Appendix B, with the transition matrix element 



of Eq. (3801 given by 



Vf^ / 1 n,nt I i > 



Pf^{r)(|)fla{r,r') - -j/i(r) • Apa{r,r') 



d^rd^r', 
(113) 

where (ppai^^, r') and A^a(r, r') being the electromagnetic 
potential generated by scattering waves Xa,i3{^') in the 
initial and final scattering states, a and /3, respectively. 
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FIG. 6: Cross sections for the Coulomb excitation versus 
incident beam energy for different collective states in the 
'^"S +^^^Au reaction, assuming a minimum impact parameter 
fcmin = 16 fm (from Ref [GlaOl]). 



Instead of the Coulomb gauge, i.e. with the Coulomb 
potential proportional to l/|r — r'|, it is better to adopt 
the Lorentz gauge, for which the scalar potential is given 
by [AW75] 



(r,r') 



'pexJ,-^*(r')^^7-^xL+^(r'), (1 



ZK|r-r'(t), 



14) 



and a similar expression for A(r,r'), with the expres- 
sion for the transition current replacing the product of 



the outgoing and incoming waves, i.e. 



4 



)* (+) 



In Eq. (114| one uses well-known expressions for 



Coulomb waves [BD04]. One also uses the expansion 



iK|r-r'| 



r r' 



47r i K ^ Jl (Kr< ) (? < ) /iL (Kr> ) Ylm (r > ) 



LM 



(115) 

where Jl (fiL) denotes the spherical Bessel (Hankel) func- 
tions (of first kind), r> (r<) refers to whichever of r and 
r' has the larger (smaller) magnitude. Assuming that the 
projectile does not penetrate the target, one uses r> (r<) 
for the projectile (target) coordinates. 

From here on, the calculation is tedious, but straight- 
forward. A detailed description is found in Ref. [AW75]. 
More insight into this calculation are not useful because 
on can show that the quantum treatment of the relative 
motion between the nuclei does not alter the results of 
the semiclassical calculations for ^ 1. Thus, we can 
safely use the machinery of the previous sections to have 
an accurate and reliable description of Coulomb excita- 
tion at low energies. For collisions at high energies, the 
distortion of the scattering waves due to the nuclear in- 
teraction yield important modifications of the angular 
distribution of the relative motion. We will discuss this 
in more details later. 

As the bombarding energy increases Coulomb excita- 
tion predominantly favors the excitation of high lying 
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FIG. 7: Schematic description of a nuclear excitation (solid 
line) followed by ^-decay (solid wavy line), (a) Tie dashed 
lines are transitions due to internal conversion (unobserved). 
The dashed wavy line is an unobserved 7-decay. (b) Direct 
emission of an observed gamma ray. 



states, e.g. giant resonances. This is shown in figure [6] 
for the cross sections of Coulomb excitation versus in- 
cident beam energy for different collective states in the 
40g _|_i97^,j reaction, assuming a minimum impact pa- 
rameter feijjin = 16 fm. Later we will discuss more about 
the excitation of giant resonances. 



IV. ANGULAR DISTRIBUTION OF 7-RAYS 

Coulomb excitation is a useful method to obtain static 
quadrupole moments as well as the reduced probabilities 
for several nuclear transitions. In order to identify the 
multipolarity of the excitation it is often necessary to 
study the de-excitation of the excited state by measuring 
a 7-ray from its decay (see figure [?]) . 

A detailed description of 7-ray emission following ex- 
citation is give in Appendix E. The angular distribution 
of the gamma rays emitted into solid angle fly, as a func- 
tion of the scattering angle of the projectile is given 
by Eq. m\, i.e.. 



wie^) = 1 



bl"^ WPfc(cos^^) 



(116) 



fc=2,4,. 



where Pfc (cos0^) are the Legendre polynomials. The co- 
efficients bZ'^ (z9) are related to the Coulomb excitation 



amplitudes n87h and to the B- values ( 33 1 for the transi 



tion from the excited state / to a final state g. 

Figure |8] shows the Coulomb excitation of sodium by 
protons. The yield of the 446 keV 7-rays is shown 
[Tem55]. Between the resonances due to compound nu- 
cleus formation one observes a smoothly rising back- 
ground yield which may be ascribed to Coulomb excita- 
tion. It is possible to determine the multipole order of the 
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FIG. 8: Coulomb excitation of sodium by protons. The yield 
of the 446 keV 7-rays is shown [Tem55]. The dashed curves 
correspond to the cross sections expected for L — 1 and 2 on 
the basis of the observed cross section for the excitation with 
a-particles. 



Coulomb excitation by comparing with the yield observed 
in the Coulomb excitation with a-particles [Tem55] . The 
dashed curves correspond to the cross sections expected 
for L = 1 and 2 on the basis of the observed cross section 
for the excitation with a-particles. The close agreement 
of the measured cross section with the theoretical curve 
for E2 excitation also confirms that the yield away from 
resonances is primarily due to Coulomb excitation. 

Figure |9] shows the 7-ray yield from Coulomb excita- 
tion and compound nucleus formation in ^^F bombarded 
with a-particles. The dashed curve shows the yields of 
the 114 keV 7-ray from the first excited state in ^^F and 
the solid curve shows the 1.28 MeV 7-ray from the first 
excited state of ^^Ne formed by an (a,p') process on ^^F 
[She54]. For bombarding energies below 1.2 MeV, the 
penetration of the a-particle through the Coulomb bar- 
rier is very small and the cross section for compound nu- 
cleus formation is small compared to that for Coulomb 
excitation. With increasing energy, the cross section for 
compound nucleus formation, ctcn, increases rapidly and 
soon becomes larger than the Coulomb excitation cross 
section, ac- However, even for Ea ^ 2 MeV, at which 
energy the average value of ctcn is an order of magnitude 
larger than ac, the yield of the 114 keV 7-ray is only 
very little affected by the compound nucleus formation, 
since the probability that the compound nucleus decays 
by inelastic a-emission is small. Finally, for Ea > 2 
MeV, the Coulomb excitation yield of the 114 keV 7-ray 
is overshadowed by the resonance yield from compound 
nucleus formation. 
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we make a multipole expansion of (1171, then we sep- 



FIG. 9: 'y-rays from Coulomb excitation and compound nu- 
cleus formation in ^^F bombarded with a-particles. The 
dashed curve shows the yields of the 114 keV 7-ray from the 
hrst excited state in ^^F and the solid curve shows the 1.28 
MeV j-ray from the first excited state of ^^JVe formed by an 
(a,p') process on ^^F [She54]. 



arate the intrinsic from the relative motion coordinates, 
and perform the time integrals for the trajectories (we 
named them orbital integrals). However, as shown in 
Ref. [AW79], it is better to first calculate the time inte- 
grals and then perform the multipole expansion. In fact, 
this seems to be the only way to get analytical expres- 
sions, except when one uses a more complicated approach 
as we will discuss later. 



The integral in Eq. ( 78 ) yields 



) = ^-^e-/^Ko m , (119) 



(r, Lo 



where Kq is a modified Bessel function and 



(120) 



V. RELATIVISTIC COLLISIONS 
A. Multipole expansion 

When the projectile has very high energies, e.g. Eiab > 
100 MeV/nucleon, there is very little deflection of the ion 
trajectory. The recoil by the target is also minimal. We 
thus assume that the projectile moves on a straight-line 
trajectory of impact parameter 5, which is also the dis- 
tance of closest approach between the centers of mass of 
the two nuclei at time < = 0. We will consider the situ- 
ation where b is larger than the sum of the two nuclear 
radii, R, such that the charge distributions of the two nu- 
clei do not overlap at any time. We will use a coordinate 
system with origin in the center of mass of the excited 
nucleus and with z-axis along the projectile velocity v, 
which is assumed to be constant. In this coordinate sys- 
tem the electromagnetic field from this other nucleus is 
given by the Lienard-Wiechert expression 



'(r,t) = 



"fZe 



^/{x - 6)2 + + 7^(2 - vty 



(117) 



We have chosen the x-axis in the plane of the trajectory 
such that the x-component of the trajectory is b. This 
expression reduces to the non-reletivistic Coulomb field 
of an low energy charge given hy = Ze/\ r — R(i) |. The 
appearance of the factors 7 = (1 — w^/c^)"^/^ g^j-g ^j^g 
retardation (see [Ja75]). The vector potential is 



A{r,t) 



(r,t), 



(118) 



with V = vz being a constant velocity vector. 

The procedure for obtaining the excitation amplitude 
in (76 1 should be the same as adopted before: first 



In Ref. [AW79] it was shown that the multipole ex- 
pansion of (119) is given by 



{r,cj) = ^M^iM(r,c.)r£M(f), (121) 



LM 



with 



WLMir, ^) = v/l67r(2£ + 1) ( ^ ^^j j ) ^ (2M - 1) 

M+1/2 / c 



- 1 VI 



L+M 



Ze ( c 



K 



(122) 



The quantity C^j{x) is the Gegenbauer polynomial 
[AS64], while jm{Kr) is a spherical Bessel function, with 
K = uj/c. For M < 0, 



WL,^M{r, to) = {-l)^'WLM{r, to) (123) 



After separation of the internal degrees of freedom, one 
gets [AW79] 



^/2L+lK^Mf^{^TL, -M) 



tLM 



.M , ^) (-l)-i^M f ^) 



(124) 



where the transition matrix elements are given by 
Mfi{TrL,—M) are given by Eqs. 



and (891 



The functions Gt^lm have analytical expressions 
[WA79]. For the El, E2, E3, Ml and M2 multipolari- 
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ties, they are given by 



Gew{x) 
Geii{x) 



-i-V'7r(a;2 - 1); 



Ge2i{x) 
Ge2o{x) 



/ 

-x^J 7r(a;2 - 1); 



Gm22{x) = -Gm2-2{x) = ^V^^^ ^ ^' 
Gm2i{x) = Gm2-i{x) ^ -xJ -; Gm20 = 0; 
^£33(2;) = -Ge3-3{x) ^ ^^^x{x'^ -l); 

Ge32{x) = Ge3~2{x) = -*;^\/ ^(3.t2 - l)Vx^ - 1; 
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^£31(2;) = -G£;3_i(a;) 



~J-x(15x^ - 11) 
lOSV 3 ^ 



105 



(125) 



From the excitation amplitude ( |124 1 one finds the total 
cross section for exciting the nuclear state of spin // in 
collisions with impact parameters larger than R given by, 



crf^ 



Ze 

he 



J2 n''^'^-'^B{nL,U^If) 



ttLM 



G 



.LAf (^)|%Af(e(^)), (126) 



where 



m 



LUR 



(127) 



and, for M > 0, 



to 



2M 



= ne 

where all Km's are functions of ^(6) = ujh/^v 



Kli+i - —Km+iKm 



(128) 
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FIG. 10: (a) A relativistic charged projectile incident on a 
target with impact parameter larger than the strong interac- 
tion radius. A sJcetci of the electric field generated by it is 
also shown. One of the effects of this field is to induce collec- 
tive vibrations of the nuclear charges, (b ) Two pulses of plane 
wave of light which produce the same effect on the target as 
the electromagnetic field created by the projectile's motion, 
(from [BB88]). 



numbers, as shown in Ref. [BB85]. The probability for 
exciting the nuclear state of energy If is obtained directly 
from Eq. (124 1 by using 



— ^ y 

* Mi, Ms 



afi 



(129) 



As with the low-energy case, we can rewrite the final 
result in terms of the reduced transition probabilities for 
photo-excitation. It can be cast in the form 



Pf^{b,E^) - PMb,E^) 

ttL 



(130) 



where a'!^^{E^) is the photonuclear absorption cross sec 
tions for a given multipolarity ttL given by Eq. ( 103 1 



The total photonuclear cross section is a sum of all these 
multipolarities, 



B. Excitation probabilities and virtual photon 
numbers 

The theoretical results for relativistic Coulomb exci- 
tation can be rewritten in terms of the virtual photon 



^7-E<(^7)- (131) 



The functions nj^L^Ej) are called the virtual photon 
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numbers, and are given by [BB85] 



V 



7 



2 ^0 r ' 



(132) 



(133) 



(134) 



where all Km's are functions of £,{b) — Lubf-fv. 

Since all nuclear excitation dynamics is contained in 
the photoabsorption cross section, the virtual photon 
numbers, Eqs. (1321, (133) and ( |134 1, do not depend on 
the nuclear structure. They are kinematical factors, de- 
pending on the orbital motion. They may be interpreted 
as the number of equivalent (virtual) photons that hit 
the target per unit area. 

The cross section is obtained by the impact param- 
eter integral of the excitation probabilities. Eq. (130) 



shows that we only need to integrate the number of vir- 
tual photons over impact parameter. One has to intro- 
duce a minimum impact parameter 5o in the integration. 
Impact parameters smaller than bo are dominated by nu- 
clear fragmentation processes. One finds 



nL 



E 



E~ 



N^L{E^)a;''{E^), (135) 



where the total virtual photon numbers N^^l^E^) 
2Tr J db b n-^L{b, Ej) are given analytically by 



Nei{E^) 



2Z^a ,c 



r-y 



^KoK, ~ {Kl - Kl) 



TT V 
TT W 



(136) 

2 , £ 

{Kl-Kl) + e{2-v^/^fKl , 

(137) 



^KoK, - I {Kl - Kl) 



(138) 



where all if^'s are now functions of ^(-R) = ojR/^v 

The usefulness of Coulomb excitation, even in first or- 
der processes, is displayed in Eq. (1301. The field of 



a real photon contains all multipolarities with the same 
weight and the photonuclear cross section, Eq. (131) is 



a mixture of the contributions from all multipolarities, 
although only a few contribute in most processes. In 
the case of Coulomb excitation the total cross section is 
weighted by kinematical factors which are different for 
each projectile or bombarding energy. This allows one to 
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FIG. 11: Left: Electric held of a slowly moving charge. Right: 
Electric held of a charge moving with v = 0.9 (in units of c). 



disentangle the multipolarities when several ones are in- 
volved in the excitation process, except for the very high 
bombarding energies 7 3> 1 for which all virtual photon 
numbers can be shown to be all the same [BB85] to 



njrL = -Z aln [ - 

TT 



(139) 



Since ^ = loR/jc ^ 1, we have a logarithmic rise of the 
cross section for all multipolarities with 7. The imping- 
ing projectile acts like a spectrum of plane wave photons 
with helicity m = =pl. Such a photon spectrum contains 
equally all multipolarities ttL. 



C. Historical note: Fermi's forgotten papers 

In 1924, Enrico Fermi, then 23 years old, submitted a 
paper "On the Theory of Collisions Between Atoms and 
Elastically Charged Particles to Zeitschrift fiir Physik 
[Fe24]. This paper does not appear in his "Collected 
Works. Nevertheless, it is said that this was one of Fer- 
mis favorite ideas and that he often used it later in life. 
In this publication, Fermi devised a method known as 
the equivalent (or virtual) photon method, where he 
treated the electromagnetic fields of a charged particle 
as a fiux of virtual photons. It is also interesting that 
Fermi published the same paper, but in the Italian lan- 
guage, in Nuovo Cimento [Fe25] (it is rather uncommon 
that the same paper is published twice!). Ten years 
later, Weiszsacker and Williams [WW34] extended this 
approach to include ultra-relativistic particles, basically 
restoring the Lorentz 7 factors in the right places. How- 
ever, it is indisputable that Fermi's papers [Fe24,Fe25] 
introduced the ingenious virtual photon method. 

A fast-moving charged particle has electric field vectors 
pointing radially outward and magnetic fields circling it. 
The field at a point some distance away from the tra- 
jectory of the particle resembles that of a real photon. 
Thus, Fermi replaced the electromagnetic fields from a 
fast particle with an equivalent fiux of photons, as shown 



schematically in figure 10 
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FIG. 12: Longitudinal and transverse electric Geld of a 
straight-line moving charge at high energies as a function of 
time. Both variables are gauged in appropriate units, as sown 
in the axis labels. 



Fermi's virtual photon method is based on the idea 
that, when v ^ c, where c is the velocity of light, the 
electromagnetic field generated by the projectile looks 
contracted in the direction perpendicular to its motion 
(see figures 11 and 12 1 and is given by 



= 

Bt — 



Zejvt 



(62+^2^2^2)3/2' 

Ze-fh 

(62+^2^2^2)3/2' 
V ^ 

- (g) Et 



and B, = 0. 



(140) 



where the z (T) indices denote the direction parallel 
(transverse) to the velocity of the projectile. 

When 7^1, these fields will act during a very short 
time, of order 



t = 



7W 



b_ 

7c' 



(141) 



and they are equivalent to two pulses of plane-polarized 
radiation incident on the target (see fig. 10 1: one in the 



beam direction (PI), and the other perpendicular to it 
(P2). In the case of the pulse PI the equivalency is exact. 
Since the electric field in the z-direction is not accompa- 
nied by a corresponding magnetic field, the equivalency 
is not complete for pulse P2, but this will not appreciably 
affect the dynamics of the problem since the effects of the 
field Ez are of minor relevance when v = c. Therefore, 
we add a field B = vEz/c to Eq. (1401 in order to treat 



also P2 as a plane-wave pulse of radiation. This analogy 
permits one to calculate the amount of energy incident 
on the target per unit area and per by Fourier transform- 
ing the Poynting vector S = E (g) B and calculating the 
intensity of the virtual radiation, b), with E-y — hio. 
This procedure is nicely explained in Ref. [Ja75] . 

Fermi associated the spectrum of the virtual radiation 
as described above to the one of a real pulse of light in- 
cident on the target. Then he obtained the probability 
for nuclear (in fact, Fermi was interested in atomic tran- 
sitions) by a fast charge, in terms of the cross sections 
for the same process generated by an equivalent pulse of 
light, i.e. 



P(6) = J I{uj,b)aj{Ey)dE^ = J n{u},b)a^{E^) 



duj 



where a^^E^) is the photo cross-section for the photon 
energy E-^, and the integral runs over all the frequency 
spectrum of the virtual radiation. The quantities n{uj, b) 
can be interpreted as the number of equivalent photons 
incident on the target per unit area. 



Following this procedure, Fermi obtained the Eq. ( 132 1 
for n(w, b) = n^i(w, b), without the 7 factors (Fermi was 
not interested in relativistic collisions in 1924!). The 7 
factors were found in the proper places by Weizsacker 
and Williams [WW34]. It is somewhat surprising that 
Fermi, Weizsacker, and Williams obtained the "exact" 
result of the virtual photons for the El multipolarity. 
Their method was completely classical and approximate 



(in adding the B^ field). To reach Eq. (1321 some quan- 



tum mechanics was used (e.g., the continuity equation 
for the nuclear transitions). Eq. (1321 is also the re- 



sult of a multipole expansion, which was not used in the 
classical prescription. In fact, the Eqs. ( |132]|134 l are 
an improvement of Fermi's method for higher multipo- 
larities, and were obtained in Ref. [BB85]. As we show 
later, a quantum mechanical description of high-energy 
scattering leads to the same expressions in Eqs 
1341). 



(132 



D. Spectrum of virtual photons 



In Eq. (132) the first term inside parentheses comes 



from the contribution of the pulse PI, whereas the sec- 
ond term comes from the contribution of the pulse P2. 
One immediately sees that the contribution of pulse P2 
becomes negligible for 7 ^ 1. The shape of the equiv- 
alent photon spectrum for a given impact parameter 
can be expressed in terms of the dimensionless function 
(j){x) — x^K^{x), if we neglect the pulse P2. In a crude 
approximation, </> ~ for x > I, and (/) — I for a; < 1. 
This refiects in a sudden cutoff of the virtual photon spec- 
tra, as can be seen from figure [l3j This implies that, in a 
collision with impact parameter b, the spectrum will con- 
tain equivalent photons with energies up to a maximum 
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FIG. 13: Equivalent photon numbers per unit area incident 
on ^°*Pb, in a collision with ^^O at 100 MeV/nucleon and 
with impact parameter 6 = 15 fm, as a function of the photon 
energy E = hui. The curves for the El, E2 and Ml multipo- 
larities are shown. 



value of order 



E^max 
^7 



(142) 



which we call by adiabatic cutoff energy. This means that 
in an electromagnetic collision of two nuclei the excita- 
tion of states with energies up to the above value can be 
reached. We can explain this result by observing that 
in a collision with interaction time given by Eq. ( 141 ) 
only states satisfying the condition T/At ^ 1, where 
T is the period of the quantum states, will have an ap- 
preciable chance to be excited. Otherwise, the quantum 
system will respond adiabatically to the interaction. In 
a collision with a typical impact parameter of 610 fm one 
can reach states with energy around E'max ~ 2O7 MeV. 
Among the many possibilities, we cite the following: for 
= 10 — 20 MeV (already small values of 7), excitation 
of giant resonances, with subsequent nucleon emission; 
for E = 20 — 100 MeV, the quasideuteron effect which 
corresponds to a photon absorption of a correlated N-N 
pair in the nucleus; and for Ej > 100 MeV, pion produc- 
tion through A-isobar excitation which has a maximum 
a,t E ~ 200 MeV. Also the production of lepton pairs 
(e+e^, qq (mesons), etc.) are accessible with increasing 
value of 7. 



by Eqs 
nE2 — nEi 
the limit 7 



In figure 14 we show n^/, (with Z = unity) as given 



( 132]|134 l 



as a function of luR/c. We see that 
njv/i for small values of 7, in contrast to 
1 . The physical reason for these two differ- 
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FIG. 14: Equivalent photon number per unit projectile 
charge, for El, Ml and E2 radiation, and as a function of 
the ratio between R and the photon wavelength. 7 is the 
ratio of the projectile energy to its rest energy. (Here, c — 1.) 



ent behaviours of the equivalent photon spectrum is the 
following. The electric field of a charged particle moving 
at low energies is approximately radial and the lines of 
force of the field are isotropically distributed, with their 
relative spacing increasing with the radial distance (see 
figure 11). When interacting with a target of finite di- 



mension, the non-uniformity of the field inside the target 
is responsible for the large electric quadrupole interaction 
between them. The same lines of force of an ultrarela- 
tivistic (7 3> 1) charged particle appear more parallel 
and compressed in the direction transverse to the parti- 
cle's motion, due to the Lorentz contraction (see figure 
111. As seen from the target, this field looks like a pulse 
of a plane wave. But plane waves contain all electric and 
magnetic multipolarities with the same weight. This is 
the cause for the equality between the equivalent photon 
numbers as 7 — > 00. 



VI. QUANTUM TREATMENT OF 
RELATIVISTIC COULOMB EXCITATION 

In contrast to sub-barrier Coulomb excitation, at rela- 
tivistic energies the effects of the nuclear interaction are 
visible due to the distortion it causes on the scattered 
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waves. Fortunately, at high energies, this can be treated 
in a simple manner by using the eikonal approximation. 
Our discussion will be general, with very little emphasis 
on the details of the nuclear interaction. 



A. The eikonal wavefunction 

The free-particle wavefunction 



(143) 



becomes "distorted" in the presence of a potential V{r ) . 
The distorted wave can be calculated numerically by 
performing a partial wave-expansion [Ber07] solving the 
Schrodinger equation for each partial wave, i.e. 



dr 



2 + k?{r) 



where 

ki{r) = 



2^ 



E - Vir) - 



2^r2 



1/2 



(144) 



(145) 



with the condition that asymptotically ^{r ) behaves as 
p3l. 



The solution of (144) involves a great numerical effort 



at large bombarding energies E . Fortunately, at large 
energies E a very useful approximation is valid when 
the excitation energies AE are much smaller than E 
and the nuclei (or nucleons) move in forward directions, 
i.e., 6* < 1 . 

Calling r = (z,b), where z is the coordinate along 
the beam direction, we can assume that 

^(r) = e*^-^0(z,b), (146) 

where (/> is a slowly varying function of z and b , so that 

|V2(/)| < fc|V0|. (147) 

In cylindrical coordinates the Schrodinger equation 



2^ 



V^V(r) + V"(r)V'(r) ^Ei;{r) 



becomes 



2ik e 



ikz 



d(j)_ 
dz 



dz^ 



ikz ^ I ^ikz 



2m 
7i2" 



V e**^^ = 



or, neglecting the 2nd and 3rd terms because of (147) 



dz 



whose solution is 
4> = exp 



Viv)<l> 



V{h,z')dz' 



(148) 



(149) 



That is, 



where 



ip{r) — exp{ikz + ix{h,z)} , (150) 



X(b,z) = -^ V{h,z')dz' (151) 

is the eikonal phase. Given V{r ) one needs a single inte- 
gral to determine the wavefunction: a great simplification 
of the problem. 

The eikonal approximation, in the same form as given 
by eqs. (150), can be obtained from the Klein-Gordon 



equation with a (scalar) potential]/. We will use this 
approximation in several discussions later on this review. 



B. Quantum relativistic Coulomb excitation 

Defining r as the separation between the center of mass 
of the two nuclei and r' as the intrinsic coordinate of the 
target nucleus, the inelastic scattering amplitude to first- 
order is given by [BD04] 



m 



ik 
2Trhv 



d^r d^r' 



n^r, r') I <^l+\r) ^v') 



(152) 



where <i>|^7'(r) and <l>|^^-'(r) are the incoming and outgoing 
distorted waves, respectively, for the scattering of the 
center of mass of the nuclei, and 0(r') is the intrinsic 
nuclear wavefunction of the target nucleus. 

At intermediate energies, AE/Eiab ^ 1, and forward 
angles, <C 1, we can use eikonal wavefunctions for the 
distorted waves; i.e.. 



W $r'(r) =exp{-iq.r-hix(6)} , (153) 



where 



1 

hv 



uT{z',h)dz' + ixc{h) (154) 



is the eikonal-phase, q = k'— k, [/^' is the nuclear optical 
potential, and xc{b) is the Coulomb eikonal phase. 



Xc{h) - ^^1^ ln(fc6) , 



(155) 



where Zp (Zt) is the projectile (target) nuclear charges. 
The Coulomb phase, as given by the above formula, re- 
produces the Coulomb amplitude for the scattering of 
point particles in the eikonal approximation for elastic 
scattering [BD04]. Corrections due to the extended nu- 
clear charges can be easily incorporated [BD04]. Here we 
have defined the impact parameter b as b = r x z. 
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In Eq. (1521 the interaction potential, assumed to be and 



purely Coulomb, is given by Eqs. (1131 and (1141. Per- 



forming the multipole expansion as in Eq. (1151, one 
finds [BN93] 



ttLM 



X flMiq) G^LM (^) (// Mf I M{7tL, -M) \h M,) 

(156) 



where the functions Gt^lm are given in Eq. (1251. The 
function rim(g) is given by [BN93] 

nMiq) = J dbb JMiqb) Km r^j exp{ix{b)} , 

(157) 

where q — 2ksm(d/2) is the momentum transfer, 6 and 
(f) are the polar and azimuthal scattering angles, respec- 
tively. 

In the sharp-cutoff approximation, one assumes 

exp{ix{b)} = 0, if b<R, 

= 1, if b> R, (158) 

where R is the strong interaction radius, or minimum 
impact parameter. Then the integral (1571 can be per- 
formed analytically [BB85]. 

Using techniques similar to those discussed in previous 
sections, one obtains 



dndE^ 



(^7) (159) 



where d driT^L / d^ is the virtual photon number given by 
[BN93] 



driT, 



LukV L[{2L + iy.lf 



dn " \-fv J (27r)3 (L+ 1) 

X \nM{q)\^- (160) 



M 



The total cross section for Coulomb excitation can be 
obtained from eqs. ( |159 1 and (160), using the approxi- 
mation dn ~ 2TTqdq/k'', valid for small scattering angles 
and small energy losses. Using the closure relation for 
the Bessel functions. 



dq qJMqxJAiiqx') — -6{x — x'), 

X 



(161) 



one obtains 

'^^'^ (E,) = Y.^.l{E,) < {E,) , (162) 



dE^ 



where the total number of virtual photon with energy Hlo 
is given by 

^-^('"^ = (^J!p(t+'i) ^ ^^^^"^ ' 



(163) 



dbbKl,{ — ] exp{-2 xiib)} 



gAiiuj) = 27r 

\ 'y^ /J V 7'^ / 

(164) 

where X/(fo) is the imaginary part of xi^), which is ob- 
tained from Eq. (154) and the imaginary part of the 
optical potential. 



If one uses the sharp-cutoff approximation, Eq. ( 158 1, 
one recovers the result in Eq. ( |128[ ) [BB85]. 

We point out that for very light heavy ion partners, the 
distortion of the scattering wavefunctions caused by the 
nuclear field is not important. This distortion is man- 
ifested in the diffraction peaks of the angular distribu- 
tions, characteristic of strong absorption processes. If 
ZpZta S> 1, one can neglect the diffraction peaks in the 
inelastic scattering cross sections and a purely Coulomb 
excitation process emerges. One can gain insight into the 
excitation mechanism by looking at how the semiclassi- 
cal limit of the excitation amplitudes emerges from the 
general result (160). 



C. Semiclassical limit of the Coulomb excitation 
amplitudes 

If we assume that Coulomb scattering is dominant and 
neglect the nuclear phase in Eq. (154), we get 



^Aiiq) - dbb Juiqb) Km — exp{ixc{b)} ■ 
Jo V7'f/ 

. . 

This integral can be done analytically by rewriting it 



as 



^M{q) = db 61+^2, ^-^^ j^^^ ^ (166) 

where we used xc{b) = 2r] ln(fc6), with 77 = ZiZ2e'^ /Hv. 
Using standard techniques found in Ref . [GR80] , we find 

^^M(g) - 2^"' ^ r(l + M + trj)T{l + trj)A'^' 

X F(1 + M + m;; l + iiy, 1 + Af ; -A^) , (167) 

where 

qjv 



A 



(168) 



and F is the hypergeometric function [GR80] . 

The connection with the semiclassical results may be 
obtained by using the low momentum transfer limit 
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and using the stationary phase method, i.e., 



2m 
'¥W) 



1/2 



G{xq) e^-^f^") , (170) 



where 



dx 



(xo) - 



and 



0"(xo) 



d^<t) 



{xo) ■ (171) 



This result is vahd for a slowly varying function G{x). 
Only the second term in brackets of Eq. (1691 will 



have a positive {b = > 0) stationary point. Thus, 
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where 

(j,{b) = ^qb + 2ri ln(fc6) . 

The condition 4>'{bo) = implies 

_ 2?] _ gp 
° ~ 9 ~ sin(0/2) ' 



(172) 



(173) 



(174) 



where oq = Z^Zte^ j [iv^ . 

We observe that the relation (1741 is the same [with 
cot(0/2) ^ sin""'^(6'/2)] as that between impact parame- 
ter and deflection angle of a particle following a classical 
Rutherford trajectory. Also, 
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which implies that in the semiclassical limit 
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Using the above results, Eq. ( 160 1 becomes 
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If strong absorption is not relevant, the above formula 
can be used to calculate the equivalent photon numbers. 
The stationary value given by Eq. (174) means that the 



important values of b which contribute to flm{q) are those 



close to the classical impact parameter. Dropping the 
rom I 



index from Eq. (1741, we can also rewrite Eq. (177) as 
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which leads to the semi-classical expressions given in Eqs. 



(132l-(134l 



For very forward scattering angles, such that A <C 1, a 
further approximation can be made by setting the hyper- 
geometric function in Eq. ( 167 1 equal to unity [GR80] , 
and one obtains 



2+2i77 



(179) 

The main value of M in this case will be Af = 0, for 
which one gets 



noiq)^2^^^T{l + tr^)ni + tr^) (^) 
= -772 22'" r(z77) T{iTj) C^) 
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and 



\no{q)\'^v' r-Y .\ — 

\ LO J rf' sinh (Trry) 
which, for 77 3> 1, results in 
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(180) 
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(182) 



This result shows that in the absence of strong ab- 
sorption and for rj ^ I, Coulomb excitation is strongly 
suppressed a,t 6 = 0. This also follows from semiclassical 
arguments, since 9 —^ means large impact parameters, 
& ^ 1, for which the action of the Coulomb field is weak. 

The results discussed in this section will be useful to 
study the corrections for Coulomb excitation at interme- 
diate energy collisions, which we shall discuss later. But, 
before that, let us remind ourselves what are giant reso- 
nances. 



VII. EXCITATION OF GIANT RESONANCES 



A. What are giant resonances? 



Figure 15 exhibits the excitation function of photoab- 
sorption of^^'^Sn around the electric dipole giant reso- 
nance at 15 MeV. The giant resonance happens in nu- 
clei along the whole periodic table, with the resonance 
energy decreasing with A without large oscillations (see 
Figure 16) starting at A = 20. This shows that the giant 



resonance is a property of the nuclear matter and not 
a characteristic phenomenon of nuclei of a certain type. 
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The widths of the resonances are almost all in the range 
between 3.5 MeV and 5 MeV. It can reach 7 MeV in a 
few cases. 

In the Goldhaber- Teller model [GT48], the photon, 
through the action of its electric field on the protons, 
takes the nucleus to an excited state where a group of pro- 
tons oscillates in opposite phase against a group of neu- 
trons. In such an oscillation, those groups interpenetrate, 
keeping constant the incompressibility of each group sep- 
arately. A classic calculation using this hypothesis leads 
to a vibration frequency that varies with the inverse of 
the squared root of the nuclear radius, i.e., the resonance 
energy varies with A~^^^. 

In the Steinwedel- Jensen model [SJ50] developed a 
classic study of the oscillation in another way, already 
suggested by Goldhaber and Teller, in which the incom- 
pressibility is abandoned. The nucleons move inside of a 
fixed spherical cavity with the proton and neutron den- 
sities being a function of the position and time. The 
nucleons at the surface have fixed position with respect 
to each other and the density is written in such a way 
that, at a given instant, the excess of protons on one 
side of the nucleus coincides with the lack of neutrons on 
that same side, and vice- versa. Such a model leads to a 
variation of the resonance energy with A^^^^. 

If one assumes a mixed contribution of the two models, 
obtains an expression for Eqdj^ as function of the mass 
number A [Mye77], 



GDR\ 



(MeV) = 112 X \a^/^ + {AoAy/^ 



-1/2 



(183) 



where Aq = 274. This expression, with the exception 
of some very light nuclei, reproduces the behavior of the 
experimental values very well, as we can see in Figure 
16 An examination of Equation ( 183 1 shows that the 



the contribution of the Steinwedel-Jensen mode is negli- 
gible. The latter mode increases with A but it only be- 
comes predominant at the end of the periodic table, at 
A = Ao.^ 

The giant electric dipole resonance arises from an ex- 
citation that transmits 1 unit of angular momentum to 
the nucleus (A/ = 1). If the nucleus is even-even it is 
taken to a 1~ state. What one verifies is that the tran- 
sition also changes the isospin of 1 unit (AT = 1) and, 
due to that, it is also named an isovector resonance. Gi- 
ant isoscalar resonances (AT = 0) of electric quadrupole 
(A/ = 2) [PW71] and electric monopole {Al = 0) [Ma75] 
were observed in reactions with charged particles. The 
first is similar to the vibrational quadrupole state created 
by the absorption of a phonon of A = 2, since both are, 
in even-even nuclei, states of 2'^ vibration. But the gi- 
ant quadrupole resonance has a much larger energy. This 
resonance energy, in the same way argued for the dipole, 
decreases smoothly with A, obeying the approximate for- 
mula 



EGQRiMeV) ^ 62^-1/3. 



(184) 



Gamow- Teller mode prevails broadly in light nuclei, while 



In the state of giant electric quadrupole resonance 
the nucleus oscillates between the spherical (supposing 
that this is the form of the ground state) and ellipsoidal 
form. If protons and neutrons act in phase, we have 
an isoscalar resonance (AT = 0) and if they oscillate 
in opposite phase the resonance is isovector (AT = 1). 
Figure [T7| illustrates these two possible vibration modes. 

The giant monopole resonance is a very special way 
of nuclear excitation where the nucleus contracts and ex- 
pands radially, maintaining its original form but changing 
its volume. It is also called the breathing mode. It can 
also happen in the isoscalar and isovector forms. It is 
an important way to study the compressibility of nuclear 
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FIG. 17: Four stages in the vibration cycle of a giant 
quadrupole resonance. In an isoscalar resonance, protons and 
neutrons vibrate in phase, while in an isovector resonance, 
the vibrations occur in opposite phase. For opposite phases, 
when the protons are at the stage (b), the neutrons will be at 
the stage (d), and vice-versa. 



matter. Again here, the isoscalar form has a reasonable 
number of measured cases, the location of the resonance 
energy being given by the approximate expression 



E, 



GMR 



(MeV) = 80A"^/^ 



(185) 



Besides the electric giant resonances, associated to a 
variation in the form of the nucleus, magnetic giant reso- 
nances exist, involving what one calls by spin vibrations. 
In these, nucleons with spin upward move out of phase 
with nucleons with spin downward. The number of nucle- 
ons involved in the process cannot be very large because 
it is limited by the Pauli principle. 

Another important aspect of the study of the giant 
resonances is the possibility that they can be induced 
in already excited nuclei. This possibility was analyzed 
theoretically by D. M. Brink and P. Axel [Ax62] for giant 
resonances excited "on top" of nuclei rotating with high 
angular momentum, resulting in the suggestion that the 
frequency and other properties of the giant resonances 
are not affected by the excitation. A series of experiences 
in the decade of the 1980's (see Reference [BB86a]) gave 
support to this hypothesis. 

A special case happens when the giant resonance is 
excited on top of another giant resonance. Understand- 
ing the excitation of a giant resonance as the result of 
the absorption of one phonon, we can view these double 
giant resonances as states of excitation with two vibra- 
tional phonons. The double giant dipole resonance was 
observed for the first time in reactions with double charge 
exchange induced by pions in '^^S [Mo88]. As first shown 
in reference [BB86] a much better possibility to study 




FIG. 18; Cross section for the excitation of the giant dipole 
resonance followed by ^-decay to the ground state in the re- 
action 

gle 0-y — u\j cLim t/y-y 

[Bee90]. The curve is predicted for Coulomb excitation 



20Spi^(i7 0,^''0'jo) at 84 MeV/nucleon for fixed 7 an- 
90° and (j>-y ~ 270°. The points are experimental 



multiple giant resonances is by means of Coulomb exci- 
tation with relativistic heavy projectiles. Later on, this 
was indeed verified experimentally and several properties 
of multiple giant resonances have been studied theoreti- 
cally (for a theoretical review, see [BP99]). 

Figure [T8| shows the cross section for the excitation of 
the giant dipole resonance followed by 7-decay to the 
ground state in the reaction 208pb(i7o^i7o'7o) at 84 
MeV/nucleon for fixed 7 angle 0^ = 90° and (j)^ = 270°. 
The points are experimental [Bee90]. The curve is pre- 
dicted for Coulomb excitation using the eikonal wave- 
function, as described in the previous section, and a re- 
duced transition strength calculated according the de- 
formed potential model [BN93] . The agreement with the 
data is excellent. 



B. Sum Rules 

It is useful to be able to estimate the total pho- 
toabsorption cross section summed over all transitions 
|i) — + I/) from the initial, for instance ground, state. 
Such estimates are given by the sum rules (SR) which 
approximately determine quantities of the following type: 

^1"^ - ^ E(^/ - ^'')"{ I (/i^i*) r + 1 (/i^^ f} ■ 

(186) 

Here the transition probabilities for an arbitrary pair 
of mutually conjugate operators F and are weighted 
with a certain (positive, negative or equal to zero) power 
n of the transition energy. Fo r a hermitian operator 
F — F^ the two terms in ( 186 1 are equal and the fac- 
tor 1/2 is cancelled. 

The exact result follows immediately for non-energy- 
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weighted SR, n = 0, based on the completeness of the set 
of the states |/), 



sI°\f] = -(i\F''F + FF''\ 



1 



(187) 



Often it turns out to be possible to get a good estimate 
for the expectation value in the right hand side of (1871, 
or to extract it from data. 

For the energy-weighted sum rules (EWSR), iS*-^', a 
reasonable estimate can be derived for many operators 
under certain assumptions about the interactions in the 
system. First, using again the completeness of the inter- 
mediate states I/), we can identically write down 5*^^^ as 
an expectation value in the initial state \i) of the double 
commutator 



sI'\f] 



[F,H],F^ 



(188) 



where H is the total hamiltonian which has energies Ei 
and Ef ds, its eigenvalues. Thus, we again need to know 
the properties of the initial state only. Now we choose 
the operator F as a one-body quantity depending on co- 
ordinates Ta of the particles. 



(189) 



Apart from that we assume that the hamiltonian does not 
contain momentum-dependent interactions. Then only 
the kinetic part of the hamiltonian contributes to [F, H\ , 
and the result can be found explicitly, 

ah a 

(190) 

where [..., ...]_|_ denotes the anticommutator. The outer 
commutator in (188) leads now to the simple result 



2ma 



^afa 



(191) 



As an example we take the charge form factor 

F = ^eae*('^-^"). (192) 

a 

The sum rule in this case is universal for any initial state 



(193) 



Taking k along an (arbitrary) z-axis and considering the 
long wavelength limit, kR <C 1, we get from the first 
nonvanishing term in the expansion of the exponent the 
EWSR for the dipole operator, = ryio(f), 



2ma 



(194) 



This is an extension of the old Thomas- Reiche-Kuhn 
(TRK) dipole SR in atomic physics. For a neutral atom, 
in the center-of-mass frame attached to the nucleus of 
charge Z (here m is the electron mass). 



2m 



(195) 



The atomic TRK SR is essentially exact (up to relativis- 
tic velocity-dependent corrections). 

In (1921 Ca are in fact arbitrary numbers. For intrinsic 



dipole excitations we have to exclude the center-of-mass 
motion. Therefore our z-coordinates should be intrinsic 
coordinates, Za ^ Za — Rz, where Rz = -^a/^- Hence, 
the intrinsic dipole moment is 



dz^ 
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This operator can be rewritten as 
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where protons and neutrons carry effective charges 
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z 



Now ( |194| gives the dipole EWSR 
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(201) 



where tojv is the nucleon mass. 

The factor (NZ/A) is connected to the reduced mass 
for relative motion of neutrons against protons as re- 
quired at the fixed center of mass. This result does not 
include the dipole strength related to nuclear motion as 
a whole. According to the classical SR ( 195 1, this contri- 
bution is (m Amj^, e Ze) 



5*-^ c.m. = -— 

^ ^ 2AmM 



(202) 



The sum of the global (2021 and intrinsic (201) dipole 
strength recovers the full TRK SR ^951), 



5f^[tot. dip] = 
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2mN V ^ 
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2mjv 



(203) 



In contrast to the atomic TRK case, the nuclear dipole 
EWSR (201) cannot be exact. Velocity-dependent and 



exchange forces are certainly present in nuclear interac- 
tions. Nevertheless, Eq. (201) gives a surprisingly good 
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estimate of the realistic dipole strength which is not fully 
understood. In a similar way one can consider SR for 
other multipoles but the results are not universal and in 
general depend on the initial state. 

The EWSR (201 1 is what we need to evaluate the sum 



of integral dipole cross sections for real photons over all 
possible final states |/). Taking the photon polarization 
vector along the z-axis, we come to the total dipole pho- 
toabsorption cross section 
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t'< — 
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This universal prediction, 



ZN 

cr7„t = 0.06— barn • MeV, 



(204) 



(205) 



on average agrees well with experiments in spite of crude- 
ness of approximations made in the derivation. One 
should remember that it includes only dipole absorption. 

For the E2 isoscalar giant quadrupole resonances one 
has the approximate sum rule [BM75] 
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{E^) ~ Q.22ZA^/'^ ^b MeV" 
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C. Coulomb excitation of giant resonances 

A simple estimate of Coulomb excitation of giant res- 
onances based on sum rules can be made by assuming 
that the virtual photon numbers vary slowly compared 
to the photonuclear cross sections around the resonance 
peak. Then 
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In figure [19] we show the Coulomb excitation cross sec- 
tion of giant resonances in '"^Ca projectiles hitting a ^^Sy 
target as a function of the laboratory energy per nucleon. 
The dashed Une corresponds to the excitation of the gi- 
ant electric dipole resonance, the dotted to the electric 
quadrupole, and the lower line to the magnetic dipole 
which was also obtained using a sum-rule for Ml excita- 
tions [BB88] . The solid curve is the sum of these contri- 
butions. 

The cross sections increase very rapidly to large values, 
which are already attained at intermediate energies. A 
salient feature is that the cross section for the excitation 
of giant quadrupole modes is very large at low and inter- 
mediate energies, decreasing in importance (about 10% 
of the El cross section) as the energy increases above 1 
GeV/nucleon. This occurs because the equivalent pho- 
ton number for the E2 multipolarity is much larger than 
that for the El multipolarity at low collision energies. 
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FIG. 19: Coulomb excitation cross section of giant resonances 
in *"Ca projectiles hitting a ^-^^U target as a function of the 
laboratory energy per nucleon. The dashed line corresponds 
to the excitation of the giant electric dipole resonance, the 
dotted to the electric quadrupole, and the lower line to the 
magnetic dipole. The solid curve is the sum of these contri- 
butions. 



That is, nE2 3> uei, for t; <C c. This has a simple expla- 
nation, as we already discussed in connection with figure 
[TT] Pictorially, as seen from an observator at rest, when 
a charged particle moves at low energies the lines of force 
of its corresponding electric field are isotropic, diverging 
from its center in all directions. This means that the 
field carries a large amount of tidal {E2) components. 
On the other hand, when the particle moves very fast its 
lines of force appear contracted in the direction perpen- 
dicular to its motion due to Lorentz contraction. For the 
observator this field looks like a pulse of plane waves of 
light. But plane waves contain all multipolarities with 
the same weight, and the equivalent photon numbers be- 
come all approximately equal, i.e., uei — nE2 — ^mi, 
and increase logarithmically with the energy for 7^1. 
The difference in the cross sections when 7 S> 1 are 
then approximately equal to the difference in the rela- 
tive strength of the two giant resonances a^'^/a^^ < 0.1. 
The excitation of giant magnetic monopole resonances is 
of less importance, since for low energies umi nsi 
{riMi — (v/cj'^nEi), whereas for high energies, where 
Umi — tlei, it will be also much smaller than the ex- 
citation of electric dipole resonances since their relative 
strength a^^^ /ai^^ is much smaller than unity. 

At very large energies the cross sections for the 
Coulomb excitation of giant resonances overcome the nu- 
clear geometrical cross sections. Since these resonances 
decay mostly through particle emission or fission, this 
indicates that Coulomb excitation of giant resonances is 
a very important process to be considered in relativistic 
heavy ion collisions and fragmentation processes, espe- 
cially in heavy ion colliders. At intermediate energies 
the cross sections are also large and this offers good pos- 
sibilities to establish and study the properties of giant 
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resonances. 

The formalism described in the previous sections has 
also been used in the analysis of the data of Ref. [BB90] , 
in which a projectile of ^^O with an energy of Eiab — 84 
MeV/nucleon excites the target nucleus ^osp]^ |-]^g Qj_ 
ant Dipole Resonance (GDR). The optical potential has 
a standard Woods-Saxon form with parameters given 
in Ref. [BB90]. In order to calculate the inelastic 
cross section for the excitation of the GDR, one can use 
a Lorentzian parameterization for the photoabsorption 
cross section of ^"^Pb [Ve70] , 



(208) 



with chosen to reproduce the Thomas-Reiche-Kuhn 



sum rule for El excitations, and Eq. (2061 for isoscalar 
giant quadrupole resonances. We use the widths FGDn — 
4 MeV and Tcqr = 2.2 MeV for ^o^Pb. 

At ~ 84 MeV/nucleon the maximum scattering angle 
which still leads to a pure Coulomb scattering (assuming 
a sharp cut-off at an impact parameter h = Rp + Rt) 
Inserting this form into Eq. (162) and doing 



4°. 



IS 

the calculations implicit in Eq. (1601 for duEi/dU,, one 
obtains the angular distribution which is compared to 
the data in Fig. |2T] The agreement with the data is 
excellent, provided one adjusts the overall normalization 
to a value corresponding to 93 % of the energy weighted 
sum rule (EWSR) in the energy interval 7—18.9 MeV (see 
section 6.10). Taking into account the ±10% uncertainty 
in the absolute cross sections quoted in Ref. [Bar88], 
this is consistent with photoabsorption cross section in 
that energy range. 

To unravel the effects of relativistic corrections, one 
can repeat the previous calculations unplugging the fac- 
tor 7 = (1 — which appears in the expressions 
(1631 and (1641 and using the non-relativistic limit of 
the functions Geitti of Eq. (125). These modifications 
eliminate the relativistic corrections on the interaction 
potential. The result of this calculation is shown in fig. 



20 (dotted curve). For comparison, the result of a full 
calculation, keeping the relativistic corrections (dashed 
curve), is also shown. One observes that the two results 
have approximately the same pattern, except that the 
non-relativistic result is slightly smaller than the rela- 
tivistic one. In fact, if one repeats the calculation for the 
excitation of GDR using the non-relativistic limit of eqs. 
( 163 1 and ( 164 1, one finds that the best fit to the data is 



obtained by exhausting 113% of the EWSR. This value 
is very close to the 110% obtained in Ref. [Bar88]. 

In fig. [20] we also show the result of a semiclassi- 
cal calculation (solid curve) for the GDR excitation in 
lead, using Eq. (177) for the virtual photon numbers. 



The semiclassical curve is not able to fit the experimen- 
tal data in figure |2T] which shows a perfect agreement 
with the Coulomb excitation with the eikonal approxi- 
mation [BN93]. This is mainly because diffraction ef- 
fects and strong absorption are not included. But the 
semiclassical calculation displays the region of relevance 
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FIG. 20: Virtual photon numbers for the electric dipole mul- 
tipolarity generated by 84A MeV O projectiles incident on 
^"^^Pb, as a function of the center-of-mass scattering angle. 
The solid curve is a semiclassical calculation. The dashed 
and dotted curves are eikonal calculations with and without 
relativistic corrections, respectively. 
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FIG. 21: Differential cross section for the excitation of the 
isovector giant dipole resonance in ^"^Pb by means of ^^0 
projectiles at 84 MeV/nucleon, as a function of the center- 
of-mass scattering angle. Data are from Ref. [BB90]. The 
theoretical calculation is from Ref. [BN93]. 
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FIG. 22; Differential cross section for the excitation of a giant 
dipole resonance in Pb+Pb collisions at 640 MeV/nucleon, as 
a function of the center of mass scattering angle. 
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FIG. 23: Cross section for the excitation of the GDR with- 
out the detection of the decay photon. Data are from Ref. 
[Bee90]. The solid line is a calculation from Ref. [BN93]. 



D. Excitation and photon decay of the GDR 



for Coulomb excitation. At small angles the scattering 
is dominated by large impact parameters, for which the 
Coulomb field is weak. Therefore the Coulomb excitation 
is small and the semiclassical approximation fails. It also 
fails in describing the large angle data (dark-side of the 
rainbow angle), since absorption is not treated properly. 
One sees that there is a "window" in the inelastic scat- 
tering data near = 2 — 3° in which the semiclassical 
and full calculations give approximately the same cross 
section. 

In fig. |22] we show a similar calculation, but for the 
excitation of the GDR in Pb for the collision ^°®Pb + 
2°®Pb at 640 MeV/nucleon. The dashed line is the re- 
sult of a semiclassical calculation. Here we see that a 
purely semiclassical calculation, is able to reproduce the 
quantum results up to a maximum scattering angle 9„i, 
at which strong absorption sets in. This justifies the use 
of semiclassical calculations for heavy systems, even to 
calculate angular distributions. The cross sections in- 
crease rapidly with increasing scattering angle, up to an 
approximately constant value as the maximum Coulomb 
scattering angle is approached. This is explained as fol- 
lows. Very forward angles correspond to large impact 
parameter collisions in which case Lob/^v > 1, the vir- 
tual photon numbers in Eq. ( 163 1 drop quickly to zero. 



and the excitation of giant resonances in the nuclei is not 
achieved. As the impact parameter decreases, increasing 
the scattering angle, Lub/jv < 1 and excitation occurs. 



We now consider the excitation of the target nucleus 
to the giant dipole resonance and the subsequent photon 
decay of that excited nucleus, leaving the target in the 
ground state. Experimentally, one detects the inelasti- 
cally scattered projectile in coincidence with the decay 
photon and demands that the energy lost by the projec- 
tile is equal to the energy of the detected photon. To 
the extent that the excitation mechanism is dominated 
by Coulomb excitation, with the exchange of a single vir- 
tual photon, this reaction is very similar to the photon 
scattering reaction, except that in the present case the in- 
cident photon is virtual rather than real. In this section, 
we investigate whether the connection between these two 
reactions can be formalized. 

We have seen that, under the conditions AE/E ^ 
1, the cross section for excitation of the target nucleus 
factorizes into the following expression (we assume that 
the contribution of the i?l-multipolarity is dominant): 



dndE^ 



1 dn-y 

EL ~dn 



(E^) {E,) 



(209) 



Figure |23] the cross section for the excitation of the 
GDR without the detection of the decay photon. Data 
are from Ref. [Bee90]. The solid line is a calculation 
from Ref. [BN93] using the eikonal wavefunction for the 
Coulomb excitation. 

The usual way to write the cross section d^ac^y/dildE^ 
for the excitation of the target followed by photon decay 
to the ground state is simply to multiply the above ex- 
pression by a branching ratio R-y, which represents the 
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probability that the nucleus excited to an energy will 
emit a photon leaving it in the ground state [Bee90] : 



dndE^ 



1 dn^ 

'K ~dn 



(E^) a^{E^) R^{E^). (210) 



Instead, one can also use the following expression, in com- 
plete analogy with Eq. ( 209 1 : 



dndE-y ^ ''' E^ 



flT) 

-^{E,) a,,{E,), (211) 



where a^j (Ej) is the cross section for the elastic scatter- 
ing of photons with energy Ej. Formally, these expres- 
sions are equivalent in that they simply define the quan- 
tity R~f. However, if one treats literally as a branch- 
ing ratio, then these expressions are equivalent only if 
it were true that the photon scattering cross section is 
just product of the photoabsorption cross section and the 
branching ratio. In fact, it is well-known from the the- 
ory of photon scattering that the relationship between 
the photoabsorption cross section and the photon scat- 
tering cross section is more complicated. In particular, 
it is not correct to think of photon scattering as a two- 
step process consisting of absorption, in which the target 
nucleus is excited to an intermediate state of energy Ej, 
followed by emission, in which the emitted photon has 
the same energy E^. Since the intermediate state is not 
observable, one must sum over all possible intermediate 
states and not just those allowed by conservation of en- 
ergy. Now, if the energy E^ happens to coincide with 
a narrow level, then that level will completely dominate 
in the sum over intermediate states. In that case, it is 
proper to regard the scattering as a two-step process in 
the manner described above, and the two expressions for 
the cross section will be equal. However, for E^ in the nu- 
clear continuum region (e.g., in the region of the GDR), 
this will not be the case. 

The calculation [BN93] is compared to the data in fig. 
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The left part shows the cross section for the ex- 
citation of the GDR without the detection of the de- 
cay photon. The agreement with the data is excellent. 
The right part of fig. |24] shows the cross section for the 
excitation-decay process as a function of E^. Although 
the qualitative trend of the data are well described, the 
calculation systematically overpredicts the cross section 
on the high-energy side of the GDR (solid curve) . If the 
Thompson amplitude is not included in ajj, the calcu- 
lation is in significantly better agreement with the data 
(dashed curve). 



E. Multiphonon resonances 

The Coulomb excitation of harmonic vibrations is for- 
mulated in Appendix F. The formalism can be applied 
to excitation of multiple giant resonances, i.e. giant res- 
onances excited on top of other giant resonances, also 
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FIG. 24: Cross section for excitation followed by 7-decay of 
^°**Pb by projectiles at 84 A MeV. The solid (dashed) 
line includes (excludes) the Thompson scattering amplitude. 
Data are from Ref. [Bee90]. 



called multiphonon giant resonances. Coulomb excita- 
tion of multiphonon giant resonances was formulated in 
the harmonic vibrator model in Ref. [BB86] . The excita- 
tion of multiphonon states may be viewed as the absorp- 
tion by the target (projectile) of several photons from the 
pulse of equivalent photons generated by the relativistic 
projectile (target) ; or it can be also described as multiple 
excitation of a harmonic oscillator with GDR quantum 
energy. 

As explained in Appendix F, the probability for the 
excitation of a iV-phonon state, Ppfib) can be obtained 
from the first-order theory for the probability, P^.o. (6), 

- ^ \Pist{b)f exp i-Pistib)] . (212) 

Pist(5) is calculated as in explained in the previous sec- 
tions. 

The multiphonon states correspond to large-amplitude 
vibrations of neutrons against protons in nuclei. Al- 
though the energy deposit is small (multiples of the en- 
ergy of a GDR state), such large collective motion may 
lead to exotic decays of the nuclei. This is a relatively 
cold fragmentation process, in contrast to the violent 
fragmentation following central collisions of relativistic 
heavy ion collisions, in which high pressures and temper- 
atures are achieved. Figure 25 shows the neutron and 
proton matter distribution in ^•^*U for a = 4 multi- 
phonon GDR state at the moment of maximum displace- 
ment [VGB90]. 

The first Coulomb excitation experiments for the exci- 
tation of the double giant dipole resonance (DGDR) were 



30 




-5 5 

radial distance 

FIG. 25: Proton and neutron density distributions at the mo- 
ment of largest separation between them in a N = A multi- 
phonon state in '^^^ U [VGB90]. 



states. Generally, the GRs in heavy nuclei decay by neu- 
tron emission. One expects that the double, or triple, 
GDR decays mainly in the 2n and 3n decay channel. In 
fact, such a picture has been adopted in [ABS93] with 
success to explain the total cross sections for the neutron 
removal in peripheral collisions. 

Coulomb excitation of multiphonon resonances is a rich 
playground for theories of nuclear response. There are 
still open issues related to the magnitude of the cross 
sections, the widths of the multiphonon states and the 
excitation mechanisms in general. Some of the reviews 
on the subject are Refs. [ABE98,BP99]. 



VIII. INTERMEDIATE ENERGY COLLISIONS 
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FIG. 26: Experimental resuhs for ^'^^Xe projectile excitation 
(at 690 MeV/nucleon) on a Pb target (squares) and a C tar- 
get (circles). The spectrum for the C target is multiplied by 
a factor 2 for better presentation. The resonance energies 
for one and two-phonon giant resonances are indicated. The 
dashed curve reflects the results of a hrst-order calculation for 
the Pb target. The figure is taken from Ref. [Sch93]. 



performed at the GSI facility in Darmstadt/Germany 
[Rit93,Sch93]. In Fig. 26 we show the result of one of 



these experiments [Sch93|7 which looked for the neutron 
decay channels of giant resonances excited in relativistic 
projectiles. The excitation spectrum of relativistic ^^^Xe 
projectiles incident on Pb are compared with the spec- 
trum obtained in C targets. A comparison of the two 
spectra immediately proofs that nuclear contribution to 
the excitation is very small. Another experiment [Rit93] 
dealt with the photon decay of the double giant reso- 
nance. A clear bump in the spectra of coincident photon 
pairs was observed around the energy of two times the 
GDR centroid energy in 208p|-, targets excited with rela- 
tivistic ^°^Bi projectiles. 

Much of the interest on multiphonon resonances relies 
on the possibility of looking at exotic particle decay of 
these states. Since there is more energy deposit in the 
nuclei, other decay channels are open for the multiphonon 



The semiclassical theory of Coulomb excitation in low 
energy collisions accounts for the Rutherford bending of 
the trajectory, but relativistic retardation effects are ne- 
glected [AW75]. On the other hand, in the theory of 
relativistic Coulomb excitation [WA79] recoil effects on 
the trajectory are neglected (one assumes straight-line 
motion) but retardation is handled correctly. In fact, 
the onset of retardation brings new important effects 
such as the steady increase of the excitation cross sec- 
tions with bombarding energy. In a heavy ion colli- 
sion around lOOA MeV the Lorentz factor 7 is about 
1.1. Since this factor enters the excitation cross sec- 
tions in many ways, such as in the adiabacity parameter, 
^(i?) = LOfiR/jv, one expects that some sizable modifi- 
cations in the theory of relativistic Coulomb excitation 
should occur [AB89]. Recoil corrections are not negligi- 
ble either, and the relativistic calculations based on the 
straight-line parametrization should not be completely 
appropriate to describe the excitation probabilities and 
cross sections. The Coulomb recoil in a single collision is 
of the order of = Z\Zie} j vnQV^ , which is half-distance 
of closest approach in a head-on collision, with mg equal 
to the reduced mass of the colliding nuclei. Although 
this recoil is small for intermediate energy collisions, the 
excitation probabilities are quite sensitive to it. This is 
important for example in the excitation of giant reso- 
nances because the adiabacity parameter is of the order 
of one. When ^(6) ^ 1, the excitation probabilities de- 
pends on h approximately like 1/6^, while when ^(6) be- 
comes greater than one they decrease approximately as 
g-2ir{(b)^j,2^ Therefore, when ^ ~ 1 a slight change of h 
may vary appreciably the excitation probabilities. 



A. Recoil and retardation corrections 

In intermediate energy collisions, where one wants to 
account for recoil and retardation simultaneously, one 
should solve the general classical problem of the motion 
of two relativistic charged particles. This has been stud- 
ied in Refs. [AB89,AAB90]. But, even if radiation is 
neglected, this problem can only be solved if one particle 
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has infinite mass [LL79]. This approximation should be 
suflicient if we take, e.g., the collision ^^O + ^'^'^Pb as our 
system. 

In the classical one-body problem, one starts with the 
relativistic Lagrangian 



1 



1/2 



, (213) 



where f and (j) the radial and the angular velocity 
of the particle, respectively. Using the Euler-Lagrange 
equations one finds three kinds of solutions, depending 
on the sign of the charges and the angular momentum in 
the collision. In the case of our interest, the appropriate 
solution relating the coUisional angle </> and the distance 
r between the nuclei is [LL79] 



1 



= A[e cos(VK0) - 1] 



where 



W = 



cLn 



1/2 



A = 



ZiZ2e^E 



cLq 



ZiZ2e^E 



mic^ + ( 



m^cZiZie^ 



(214) 



(215) 

-,1/2 

(216) 



E is the total bombarding energy in MeV, mo is the mass 
of the particle and Lq its angular momentum. In terms 
of the Lorentz factor 7 and of the impact parameter 6, 
E ~ ^maC? and Lq — ^mQvb. The above solution is 
valid if Lq > Z\Z2e^ j c. In heavy ion collisions at in- 
termediate energies one has Lq 3> Z\Z2e^ jc for impact 
parameters that do not lead to strong interactions. It is 
also easy to show that, from the magnitudes of the pa- 
rameters involved in heavy ion collisions at intermediate 
energies, the trajectory (214) can be very well described 
by approximating 



W ^\ 



762 ' 



(217) 



where ao is half the distance of closest approach in a 
head on collision (if the nuclei were pointlike and if non- 
relativistic kinematics were used), and e is the eccentric- 
ity parameter. In the approximation (217) e is related to 
the deflection angle by e = (ao/7) cot d. 

The time dependence for a particle moving along the 
trajectory, Eq. (214), may be directly obtained by solv- 
ing the equation of angular momentum conservation. In- 
troducing the parametrization 



[e cosh X + 1] 



we find 



[x + e sinh x] 



(218) 



(219) 



Using the scattering plane perpendicular to the z-axis, 
one finds that the corresponding components of r may 
be written as 



X ~ a [cosh X -l- e] 



. \/ e'^ — 1 sinh x 



- , 



(220) 



where a — ao/7. This parametrization is of the same 



form as in the non-relativistic case, Eq. (83), except that 
ao replaced by ao/7 = a. 

Non-relativistically the two-body problem is solvable 
by introduction of center of mass and relative motion 
coordinates. Then, the result is equivalent to that of a 
particle with reduced mass too = mprriT / {mp -\- ttlt) 
under the action of the same potential. The particle 
with reduced mass mo is lighter than those with mass 
mp and m^, and this accounts for the simultaneous re- 
coil of them. An exact relativistic solution should re- 
produce this behavior as the relative motion energy is 
lowered. We can use the reduced mass definition of mo 
as usual in the parametrization of the classical trajec- 
tory of Coulomb excitation in intermediate energy col- 
lisions, as outlined above. For a ^^O -I- ^''^Pb collision 
this is not a bad approximation. For heavier systems like 
U-l-U it would be the simplest way to overcome this dif- 
ficulty. But, as energy increases, this approximation is 
again unimportant since the trajectories will be straight- 
lines parametrized by an impact parameter b. An exact 
result was obtained numerically in Ref. [AAB90] using 
the Darwin Lagrangian to determine the classical tra- 
jectory in collisions at intermediate energies. But, the 
parametrization of the classical trajectory as given by 
Eqs. (220) with a reduced mass particle, besides repro- 



ducing both the non-relativistic and the relativistic ener- 
gies, gives a reasonable solution to the kind of collisions 
we want to study. 

From here on the procedure is the same as in the 
derivation of low-energy excitation amplitudes. When- 
ever ao appears in those equations it will have to be re- 
placed by ao/7. The orbital integrals in Eqs. (90) and 



(91), will be done along the trajectory defined by Eqs. 
220 1 . Explicitly, one finds that the orbital integrals in 



Eq. (98) becomes [AB89] 



-71-17/2 



(e -t- ismhx — ^/ — \ cosh %)'' 
{iesmhx + 1)'^^-'- 

(A -I- 1) hx - zhx+i - {-Ye-q coshx ■ hx 
c 



where hxs are spherical Hankel functions, and 

V 

z = — rj (iesinhx + 1) ■ 
c 



(221) 



(222) 
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FIG. 27: Electric dipole number of equivalent photons per 
unit area d?h = 2-Kbdh, with energy of 10 MeV, incident on 
^°*Pb in a collision with at impact parameter b — 15 fm, 
and as a function of the bombarding energy in MeV per nu- 
cleon. The dotted line and the dashed line correspond to 
calculations performed with the non-relativistic and with the 
relativistic approaches, respectively. The solid line represents 
a more correct calculation, as described in the text. 



For magnetic excitations [AB89], 



X e 



(2A-1)!! 

ivx + * ^ ^ Ve^ - 1 cosh x)'^ 
(iesinh x + 1)^ 



(223) 



The excitation amphtudes are still given by Eq. (87 1 
with the expressions for S{ttLM) obtained from Eqs. (90 
91 1 by replacing ao by ao/7 and using e = (Q0 /7) c ot d 
The excitation cross sections are given by Eq. ( |102 1 and 
the virtual photon numbers by Eqs. ( 105]|106[ ). 

As an example, we consider the excitations induced in 
208p|-, jj-^ 16 Q 208pj^ collisions. Since the expression 
(1021, with Eqs. ( 105|106 l, is quite general, valid for all 
energies, the equivalent photon numbers contain all infor- 
mation about the differences among the low and the high 
energy scattering. In Figs. 27 28]and 29 we show dn.^\ g, 

and Ml (Fig. ^ 



27), E2 (Fig 



160 + 208pb with 
They are the equivalent 



for the El (Fig, 

multipolarities, and for the collision 
an impact parameter 6=15 fm 
photon numbers with frequency u; = 10 MeV//i incident 
on ^°8pb. The dotted lines are obtained by using the non- 
relativistic Eq. ( 98 1 , while the dashed lines correspond to 
the relativistic expressions ( |132|133|134| . One observes 
that the relativistic expressions overestimate the equiv- 
alent photon numbers at low energies, while the non- 
relativistic expressions underestimate them at high en- 
ergies. The correct values are given by the solid lines. 
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FIG. 28: Same as Fig. but for the E2 multipolarity. 
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FIG. 29: Same as Fig. but for the Ml multipolarity 



calculated according to Eqs. (2211 and (2231. They re- 
produce the low and the high energy limits, giving an 
improved interpolation between these limits at interme- 
diate energies. These discrepancies are more apparent 
for the El and the E2 multipolarities. In the energy in- 
terval around lOOA MeV neither the low energy theory 
nor the high energy one can reproduce well the correct 
values. This energy interval is indeed very sensitive to 
the effects of retardation and of Coulomb recoil. 

At these bombarding energies, the differences between 
the magnitude of the non-relativistic and the correct 
relativistic virtual photon numbers are kept at a con- 
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stant value, of about 20%, for excitation energies Ex = 
hu) < 10 MeV. However, they increase sharply when one 
reaches the excitation energy of — hoj > 10 MeV. 
The reason is that, for such excitation energies, the adi- 
abaticity factor becomes greater than unity > 1). 
This means that excitation energies of order of 10 MeV 
(like in the case of giant resonance excitation) are in the 
transition region from a constant behavior of the equiva- 
lent photon numbers to that of an exponential e~'^^) 
decay. This is more transparent in Fig. [13] where 
we plot the equivalent photon numbers for Eiab —100 
MeV/nucleon, 6 = 15 fm, and as a function of hw. One 
also observes from this figure that the E2 multipolar- 
ity component of the electromagnetic field dominates at 
low frequencies. Nonetheless, over the range of fku up 
to some tens of MeV, the E2 matrix elements of exci- 
tation are much smaller than the El elements for most 
nuclei, and the E2 effects become unimportant. How- 
ever, such effects are relevant for the excitation of the 
isoscalar E2 giant resonance (GQRis) which have large 
matrix elements. 

Only for the El multipolarity the orbital integral in 
Eq. (221 1 can be done analytically. 
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(224) 



where f- — 1/ sin{§/2), a = 1/137, C — ujao/jv. Ki^ is 
the modified Bessel function with imaginary index, JC^'^ is 
the derivative with respect to its argument. If one com- 
pares this equation with Eq. ( 111 I one notes that besides 
the appearance of 7 factors in a = ao/7, the first term 
inside square brackets is also reduced by a 1 /j^ , as com- 
pared to the non-relativistic expression. For relativistic 
energies, ( ^ 0, e ^ 00, with (e const.. It is t hen e asy 
to see that the above expression reduces to Eq. (1321, in 
terms of the impact parameter dependence. 

The integration over impact parameter can also be 
doen analytically for the El case. One obtains [BB88] 
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where C — toao/jv, 



Co 



1, 



for 2a > bmim , 



R/a — 1, for 2a < h„ 
and ^ = eoC = iubmunliv. 



(225) 



(226) 



It is easy to see that this equation reduces to Eq. ( 136 1 
in the relativistic limit, when Q — > 0, eo — > 00. In 
practice, it is much easier to numerically calculate the 



orbital integrals in Eqs. ( 221]|223 1 directly. A numerical 
code (COULINT) is available for the purposes [Ber04]. 



IX. E&M RESPONSE OF WEAKLY-BOUND 
SYSTEMS 

Coulomb excitation has been one of the main tools to 
study the structure of unstable isotopes in radioactive 
beam facilities. This is a relatively new field and many 
new aspects of Coulomb excitation have been developed 
in the last two decades. We will discuss few of these 
aspects and, in particular, the role of transitions in/to 
the continuum of loosely-bound systems. 



A. Dominance of one resonance 

It often occurs that the Coulomb excitation mecha- 
nism is dominated by transitions between the ground 
state and a continuum resonant state. This part of 
the excitation mechanism can be described exactly in a 
coupled-channels approach, while the other excitations 
can be described perturbatively [Ca96]. Fig. 30 repre- 
sents the procedure. One resonance is coupled to the 



ground state while the remaining resonances are fed by 
these two states according to first order perturbation the- 
ory. The coupling matrix elements involves the ground 

state and a set of doorway states l-D^'^^ > , where n speci- 
fies the kind of resonance and A/i are angular momentum 
quantum numbers. The amplitudes of these resonances 
in real continuum states are 



a(")(e) =< 0(e) IP 



(n) 
A /J, 



(227) 



where (/)(e) denotes the wavefunction of one of the numer- 
ous states which are responsible for the broad structure 
of the resonance. In this equation e — E^ ~ £"„ i where 
Ex is the excitation energy and En is the centroid of 
the resonance considered. 

Let us assume for simplicity that the dominant res- 
onant state has spin and parity = 1^ . However, 
the following results can be easily generalized to all spin- 
parity types. The following coupled-channels equations 
are obtained [Ca96]: 
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X exp 
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X exp^--iE + e)t\ al^^t), (228) 
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FIG. 30: Schematic representation of Coulomb excitation of 
a dominant resonant state and of other states weakly-coupled 
to the ground state. 



A better agreement with the experimental line shapes 
of the giant resonances is obtained by using a Lorentzian 
(L) parametrization for |a^^-'(e)p , i.e., 



[El-ElY + TlEl 



(233) 



where Ex = i?i + e. The energy integral can still be per- 
formed exactly, but now it leads to the more complicated 
result 



de |a(i)(e)l' 



exp 
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(234) 



and 



ih 



da 



(1) 

e,l/j- . 



dt 



(t) = [(«(!)(£) F("i)(i) 



X exp<j -{Ei + e)t)- ao{t). (229) 



Above, (n = 1) stands for the dominant resonance, 
aodenotes the occupation amplitude of the ground state 
and a[^l^ the occupation amplitude of a state located 
at an energy e away from the centroid of the reso- 
nance, and with magnetic quantum number fi ifJ- = 
— 1,0,1 ). We used the short hand notation v/f^\t) = 



Integrating Eq. (229 1 and inserting the result in 



Eq. ( 228 ) , we get the integro-differential equation for the 
ground state occupation amplitude 



d'^ao 
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exp 



-{E, + e){t~t')\ao{t'), 
(230) 



where we used that 



(i — — oo) — 0. To carry out 



the integration over e , we should use an appropriate 
parametrization for the doorway amplitude a'^^(e) . A 
convenient choice is the Breit-Wigner (BW) form which 
yields the square amplitude 
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(231) 



where Fi is chosen to fit the experimental width. In this 
case, this integral will be the simple exponential 



de |a(i)(e)l' exp 



.(Sl-^^l/2)^ 



.{Ei + e)t 



exp < —t 



(232) 



where AC(t) is a non-exponential correction to the de- 
cay. This correction can be shown numerically to be neg- 
ligible. It will therefore be ignored. After integration over 
e , Eq. ( 230 1 reduces to 



d'^ao 



it) 



(Ei-iTi/2)(t-t') \ , ,^ 
X exp<;-i^— ^ -{—^ -^ao{t') 



(235) 

where the factor Si is 5*1 = 1 for BW-shape and 5*1 = 
1 — iTi/2Ei for L-shape. 

We can take advantage of the exponential time- 
dependence in the integral of the above equation, to re- 
duce it to a set of second order differential equations. 
Introducing the auxiliary amplitudes Af^ (t) , given by the 
relation 



ao{t) = l + J2Mt): 



(236) 



with initial conditions A^(t = — oo) — , and taking the 
derivative of Eq. (2361, one obtains 



El 



5,^ 



(01) 



l + Y,A,.{t) 



= 0. 



Mi) 



(237) 



Solving the above equation, we get ao(t). Using this 
amplitude and integrating Eq. (229), one can evaluate 



a"Pi^{t) ■ The probability density for the population of a 
dominant continuum state with energy E^ in a collision 
with impact parameter 5, Pi(b, E^) , is obtained trough 
the summation over the asymptotic (t oo ) contribu- 
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tion from each magnetic substate. This yields 



dB (EL) /dE^, in Eq. (|240l) is given by 



dt' 



(238) 



where \a('^>{E^ - Ei)Y is given by Eq. (|23ip or by 



Eq. ( 233 1 , depending on the choice of the resonance 
shape. 

To first order, the contribution to the excitation of 



another resonance (the one on the top of Fig. (30)) from 
the dominant one is given by 

P2{b,E^) = \a'^''\E^-E2)\''Si 



E 



X exp 



r dt'e^p{^EJ} J2 Wy,'Ht') 

X l' dt"v^''Ht") 

.{Ei-iTi/2){t-t') 



ao{t") 



(239) 



We should point out that Eq. ( 239 1 is not equivalent to 
second-order perturbation theory. This would be true 
only in the limit ao{t) — > 1 . In this approach, aoit) ^ 1, 
since it is modified by the time-dependent coupling to 
the dominant state. This coupling is treated exactly by 
means of the coupled-channels equations. This approach 
is justified due to the (assumed) small excitation ampli- 
tude for the transition 1 — > 2 , since ai(t) <C ao{t). 



Equations similar to ( 238 1 can also be used to calcu- 
late the excitation probabilities directly from the ground- 
state, with the proper choice of energies, widths, and 
transition potentials. The above formalism has been used 
to study the excitation of giant resonances multiphonon 
resonances [Ca96] . It can be easily applied to description 
of pigmy resonances in nuclei far from the stability. 



B. Cluster-like light nuclei 

In first-order perturbation theory, the photo-nuclear 
cross sections are related to the reduced matrix elements, 
for the excitation energy E^ through the relation 



(27r)3(L + l) fE. 



dE^ ^ i[(2L+l)!!]' \ hc) dE.,^ 

(240) 

where dB/dE^, are the reduced matrix elements, or 
Coulomb response functions. 
The response function for electric transitions. 



dB{EL) |(J^||ri(f)|| J- 



dEj 2Ji + 1 

drr^+^ Sp\f^^{r) 



w (E^) , (241) 



where w (E^) is the density of final states (for nuclear 
excitations into the continuum) with energy E^ = Ef — 
E,. 

(EL') 

The transition density 5p\j: (r) will depend upon the 
nuclear model adopted. We will discuss here very simple 
models to calculate this response. To get a better repro- 
duction of experimental data, more microscopic calcula- 
tions such as the Random Phase Approximation (RPA) 
are often necessary. 



1. Cluster response function 

Due to the small binding energy of nuclear halo sys- 
tems composed of two clusters (a = b + c), the most 
important part of the wavefunction is beyond the range 
of the nuclear potential. The external part of the wave- 
function is well described by an Yukawa of the form 



ijbcir) = No 



Tj e 



(242) 



where ry = yj2pLi)cS/h, fihc is the reduced mass of (b + c), 
S is the separation energy, and Nq = e"^"^" + r^r^ is a 
normalization factor which corrects the wavefunction to 
account for the finite range, tq, of the (6 -f c) potential. 
There is no resonance structure in the b + c continuum. 
This is clearly a good assumption for the deuteron and 
also for other neutron halo systems. 

If we further assume that the final state is a plane- 
wave state (i.e., we neglect final state interactions) ipf = 
(q|r) = e*'! '". The response functions for electric multi- 
pole transitions 



dB 
dE^ 



{EX,E,) = SFY,\ {ci\M{ELM)\i;bc{r)} 



2 d3 



(27r)3 
(243) 

where SF is a spectroscopic factor (i.e., probabihty to 
find the system in the state {b + cj). M{ELM) is the 
electric multipole operator corrected to exclude cm. mo- 
tion, i.e.. 



MiELfi) 



1 - 



+ (-1) 



l( Z-I) 



TiYLM 



(244) 



where the sum runs over all protons (p) and neutrons (n) 
in the nucleus. 
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For a cluster-system the sum runs over the (effective) 
charges of the clusters. For example, 



M{E1M) 
M{E2M) 



e I : I TbcYiMirbc) 



Aa 

ZbAl + Z,Al 
Al 



rbcY2M{vbc) (245) 



where Vbc is the relative position of b-c. 

Using (Pq = q^dqd^q = ^/2E^{^ibc/h'^f/'^dEdnq in 
Eq. (2431 and integrating over fi, one finds [BS92] 



dB{EL-E^) 
dK 



, {Lr)\2L+l){SF)N^[0{L)] 
' VS{E^ - ^)^+i/2 



TT 



E: 



2L+2 



where 



OiL) 



ZbA^ ~ (-l)^^cA^ 



A^ 



The maximum of this function occurs at 



1/2 



L + 3/2 



(246) 



(247) 



(248) 



In Figure [31] we show the photonuclear cross sections 
obtained by using Eq. ( 246 1 in the definition ( 240 1 . The 



E2 and E3 photonuclear cross sections are smaller than 
the El by factors of order of 10^ and 10^, respectively. 
The spectroscopic factor and the normalization constant 
{S, and No, in Eq. (2461 were set to unity, for simplicity. 



The have a longer tail than the El photonuclear cross 
section due to the factor E"^^'^ in Eq. (|240|. Thus, we 



conclude that a cluster-like correlation is not as effective 
in producing an enhancement of the low energy part of 
the photonuclear cross sections for higher multipolarities, 
as in the dipole case. 



2. One-neutron halo 

We will now discuss the effect of final state inter- 
actions. We start by considering a one-neutron halo 
nucleus. The initial wavefunction can be written as 
*JM r^^uijj{r)yijM, where Rijj{r) is the radial 
wavefunction and yijM is a spin-angle function [BD04]. 

As before, we use a simple Yukawa form for an s-wave 
initial wavefunction, UQ{r) = NoCxp{—rir), and a p-wave 
final wavefunction, ui{r) = ji{kr) cos 6i — ni(fcr) sin (Si. 
In these equations t] is related to the neutron separation 
energy Sn = ^■^?7^/2/i, n is the reduced mass of the neu- 
tron -|- core system, and hk = \/2^Er, with E^. being the 
final energy of relative motion between the neutron and 
the core nucleus. A^o is the normalization constant of the 
initial wavefunction. 

The transition density is given by r'^Spif (r) — 
eeffAiUi(r)uf(r), where i and / indices include angu- 
lar momentum dependence and Ce// = —eZc/A is the 
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FIG. 31: Photo absorption cross section of ^^Li in the clus- 
ter model. The solid (dashed) [dotted] curve corresponds to 
dipole (quadrupole) [octupole] multipolarity. The quadrupole 
(octupole) values were multiplied by a factor of 10^ (10^) 
[BS92]. 



effective charge of a neutron+core nucleus with charge 
Zc- The El transition integral T^ijf = Jq dr Spif (r) 
for the wavefunctions described above yields 



^s^p — Beff 



2fc2 



°eff 



h'^ 2Er 



2m (5„ 



Er 



COS 5i + sin 5i 



1- 



1] {ff + 3fc2) 



2/c3 



\2K^) ~l/ai+ iiriEr/ti^ 



(249) 



where the effective range expansion of the phase shift, 
^2i+i ^ ^ — 1/a; + r-/fc^/2, was used in the second line 
of the above equation. For Z = 1, ai is the "scattering 
volume" (units of length'^) and ri is the "effective mo- 
mentum" (units of 1 /length) . Their interpretation is not 
as simple as the I = effective range parameters. Typical 
values are, e.g. ai = —13.82 fm~'^ and ri = —0.419 fm^^ 
for n+^He pi/2-wave scattering and ai = —62.95 fm"^ 
and ri = —0.882 fm^^ for n+^He p3/2-wave scattering 
[Arn73]. 

The energy dependence of Eq. ( 249 1 has some unique 
features. As shown in the previous section, the matrix 
elements for electromagnetic response of weakly-bound 
nuclei present a small peak at low energies, due to the 
proximity of the bound state to the continuum. This 
peak is manifest in the response function as 



dB{EL) 2 



(250) 



(EL) 

It appears centered at the energy [BS92] E^^'^' 
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FIG. 33: Jacobian coordinates (x and y) for a three-body 
system consisting of a core (c) and two nucleons (1 and 2). 



FIG. 32: \Is^p\ calculated using Eq. ( |249[ l, assuming e^ff = 
e, A = 11 and S„ = 0.5 MeV, as a function of Er, tie relative 
energy of the fragments. The long dashed curve corresponds 
to ai = —10 fm~^ and vi = —0.5 fm~^ , the dashed curve 
corresponds to ai = —50 fm^"^ and ri = 1 fm~^ , the solid 
curve corresponds to ai = ri = 0, and finally, the dotted 
curve corresponds to ai = — 10 fm~^ and vi — 0.5 fm~^ . 



3. Two-neutron halo 



£+1/2 



L+3/2'-'" ^™ ^ generic electric response of niultipolarity 
L. For El excitations, the peak occurs at Eq ~ 3S'„/5. 
The second term inside brackets in Eq. (2491 is a 



modification due to final state interactions. This mod- 
ification may become important, as shown in figure [32] 
where |Xs^p|^ calculated with Eq. (2491 is plotted as a 
function of Er. Here, for simplicity, we have assumed the 
values Ceff = e, A = 11 and S'„ = 0.5 MeV. This does 
not correspond to any known nucleus and is used to as- 
sess the effect of the scattering length and effective range 
in the transition matrix element. The long dashed curve 
corresponds to ai = —10 fm~'^ and ri — —0.5 fm~^, the 
dashed curve corresponds to ai = —50 fm"-^ and ri — 1 
fm~^, the solid curve corresponds to ai = ri = 0, and 
finally, the dotted curve corresponds to ai — —10 fm^'^ 
and ri = 0.5 fm^^. Although the effective range expan- 
sion is only valid for small values of Er , it is evident from 
the figure that the matrix element is very sensitive to the 
effective range expansion parameters. 

The strong dependence of the response function on the 
effective range expansion parameters makes it an ideal 
tool to study the scattering properties of light nuclei 
which are of interest for nuclear astrophysics. It is impor- 
tant to notice that the one-halo has been studied in many 
experiments, e.g. for the case of ^^Be for which there are 
many data available. In the literature one can find a 
detailed analysis of how the nuclear shell-model can ex- 
plain the experimental data, by fitting the spectroscopic 
factors for several single-particle configurations. 



Weakly-bound nuclei, such He or Li, require a 
three-body treatment in order to reproduce the electro- 
magnetic response more accurately. In a popular three- 
body model, the bound-state wavefunction in the center 
of mass system is written as an expansion over hyper- 
spherical harmonics (HH), see e.g. [Zhu93], 



^'(x,y) = 



P' 



1 

J/2 



E 



JkL (f^5) ® XS 



JM 



(251) 



Here x and y are Jacobi vectors where (see figure 33 1 
r2) and y = 



2{A-2) 



ri+r2 
2 



™), wher 



A is the nuclear mass, ri and T2 are the position of the 
nucleons, and Yc is the position of the core. The hy- 
perradius p determines the size of a three-body state: 

= + y^. The five angles {^^5} include usual angles 
), {Oy,(l)y) which parametrize the direction of the 
unit vectors x and y and the hyperangle 9, related by 
X — psin0 and y = pcosO, where < < n/2. 

The insertion of the three-body wavefunction, Eq. 



(251 1, into the Schrodinger equation yields a set of cou- 



pled differential equations for the hyperradial wavefunc- 



tion ^"^'^5 (p). Assuming that the nuclear potentials be- 
tween the three particles are known, this procedure yields 
the bound-state wavefunction for a three-body system 
with angular momentum J. 

In order to calculate the electric response we need the 
scattering wavefunctions in the three-body model to cal- 
culate the integrals in Eq. (241 ). One would have to use 



final wavefunctions with given momenta, including their 
angular information. When the final state interaction 
is disregarded these wavefunctions are three-body plane 
waves [Pus96]. To carry out the calculations, the plane 
waves can be expanded in products of hyperspherical har- 
monics in coordinate and momentum spaces. However, 
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FIG. 34: Comparison between the calculation of the response 
function (in arbitrary units) with eqs. ( |252[ ) and ( |253[ ), using 
Snn = and 5„c = 0, (dashed line), or including the effects of 
final state interactions (continuous line). The experimental 
data are from Ref. [Shi95]. 



since we are only interested in the energy dependence of 
the response function, we do not need directions of the 
momenta. Thus, instead of using plane waves, we will 
use a set of final states which just include the coordinate 
space and energy dependence. 

For weakly-bound systems having no bound sub- 
systems the hyperradial functions entering the expan- 
sion (2511 behave asymptotically as [Mer74] ^aip) — > 
constant x exp (—rjp) as p — *■ oo, where the two-nucleon 
separation energy is related to rj by S'2„ = fi^if' / (2mjv). 
This wavefunction has similarities with the two-body 
case, when p is interpreted as the distance r between 
the core and the two nucleons, treated as one single par- 
ticle. But notice that the mass mjv would have to be 
replaced by 2mj4 if a simple two-body (the dineutron- 
model [BB88]) were used for "Li or ^He. 

Since only the core carries charge, in a three-body 
model the El transition operator is given by M 
yYiM (y) for the final state. The E\ transition matrix 
element is obtained by a sandwich of this operator be- 
tween <^a (p) /p^^^ and scattering wavef unctions. In Ref. 
[Pus96] the scattering states were taken as plane waves. 
We will use distorted scattering states, leading to the 
expression 



I{E1) 



dxdy ^\[^^ y^xup (y) Ug (x) , 



P 



,5/2 



(252) 



where u,, (y) = ji (py) cos (5„c - ni (py) sin S^c is the 
core-neutron asymptotic continuum wavefunction, as- 
sumed to be a p-wave, and Uq (x) — jo (qx) cos Snn ~ 
no (qx) sin Snn is the neutron-neutron asymptotic con- 
tinuum wavefunction, assumed to be an s-wave. The 
relative momenta are given by q = ^ (qi — q2), and 



The El strength function is proportional to the square 
of the matrix element in Eq. ( |252 1 integrated over all mo- 
mentum variables, except for the total continuum energy 
E,. = (g^ /2m7v. This procedure gives 



dB{El) 
~dE7~ 



constant J \I{El)fE'^cos'^Q 



sin" QdQdflqdflp, 



(253) 



where Q — tan"^ (q/p)- 

The ^Sq phase shift in neutron-neutron scattering is 
remarkably well reproduced up to center of mass energy 
of order of 5 MeV by the first two terms in the effective- 
range expansion fccot(5„„ ~ — l/a„„ -I- r„„fc^/2. Experi- 
mentally these parameters are determined to be a„„ = 
—23.7 fm and r„„ = 2.7 fm. The extremely large (neg- 
ative) value of the scattering length implies that there 
is a virtual bound state in this channel very near zero 
energy. The p-wave scattering in the n-^Li ("'^^Li) sys- 
tem appears to have resonances at low energies [Thoe99] . 
We assume that this phase-shift can be described by the 

resonance relationsin (5„c ~ (r/2)/-^ (^E^ 
with Er = 0.53 MeV and T 

Most integrals in eqs. (2521) and (|253[) can be done an- 



Enf + ryA, 

0.5 MeV [Thoe99]. 



alytically, leaving two remaining integrals which can only 
be performed numerically. The result of the calculation 
is shown in figure ( 34 1 . The dashed line was obtained us- 



ing Snn — and Snc — 0, that is, by neglecting final state 
interactions. The continuous curve includes the effects 
of final state interactions, with Snn and Snc parametrized 
as described above. The experimental data are from Ref. 
[Shi95]. The data and theoretical curves are given in 
arbitrary units. Although the experimental data is not 
perfectly described by either one of the results, it is clear 
that final state interactions are of extreme relevance. 

As pointed out in Ref. [Pus96], the El three-body 
response function of ^"'^Li can still be described by an ex- 
pression similar to Eq. ( 250 1 , but with different pow- 

11/2 



ers. 



Explicitly, dB{El)/dEr cx E^/ {s^l/' + E, 



Instead of Sor, , one has to use an effective S'. 



eff _ 



IS2 



with a ~ 1.5. With this approximation, the peak of the 
strength function in the three-body case is situated at 
about three times higher energy than for the two-body 
case, Eq. (250). In the three-body model, the maximum 

is thus predicted at Eq ' ~ 1.852n, which fits the ex- 
perimentally determined peak position for the ^^Li El 
strength function very well [Pus96]. It is thus apparent 
that the effect of three-body configurations is to widen 
and to shift the strength function dB (El) /dE to higher 
energies. 

It is worthwhile mentioning that the data presented in 
figure [34] and of other experiments, is different in form 
and magnitude of the more recent experiment of Ref. 
[Nak06] . The reason for the discrepancy is attributed to 
an enhanced sensitivity in the experiment of Ref. [Nak06] 
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to low relative energies below Erei = 0.5 MeV com- 
pared to previous experiments. Also, this recent exper- 
iment agrees very well with the nn-correlated model of 
Ref. [EB92] . This theoretical model is different than the 
model presented in this section in many aspects. In prin- 
ciple, the three-body models should be superior, as they 
include the interactions between the three-particles with- 
out any approximation. For example, Ref. [EB92] use a 
simplified interaction between the two- neutrons. On the 
other hand, they include the many-body effects, e.g. the 
Pauli blocking of the occupied states in the core. It is 
not well known the reason why the data of Ref. [Nak06] 
is better described with the model of Ref. [EB92] than 
traditional 3-body models. 





(a) 



(b) 



C. Pigmy resonances 

We have seen that the energy position where the soft 
dipole response peaks depends upon the few body model 
adopted. Except for a two-body resonance in ^°Li, there 
was no reference to a resonance in the continuum. The 
peak in the response function can be simply explained 
by the fact that it has to grow from zero at low energies 
and return to zero at large energies. In few-body, or clus- 
ter, models, the form of the bound-state wavefunctions 
and the phase space in the continuum determine the po- 
sition of the peak in the response function. Few-body 
resonances will lead to more peaks. 

Now we shall consider the case in which a collective 
resonance is present. As with giant dipole resonances 
(GDR) in stable nuclei, one believes that pygmy res- 
onances at energies close to the threshold are present 
in halo, or neutron-rich, nuclei (in English, pigmy = 
pygmy) . This was proposed in Ref. [SIS90] using the hy- 
drodynamical model for collective vibrations. The possi- 



bility to explain the soft dipole modes (figure 34 1 in terms 
of direct breakup, has made it very difficult to clearly 
identify the signature of pygmy resonances in light ex- 
otic nuclei. 

The hydrodynamical model needs adjustments to ex- 
plain collective response in light, neutron-rich, nuclei. 
Because clusterization in light nuclei exists, not all neu- 
trons and protons can be treated equally. The necessary 
modifications are straight-forward and discussed next. 

When a collective vibration of protons against neu- 
trons is present in a nucleus with charge (neutron) num- 
ber Z (N), the neutron and proton fiuids are displaced 
with respect to each other by di = aiR and each of the 
fluids are displaced from the origin (center of mass of the 
system) by dp — Ndi/A and d„ — —Zdi/A. This leaves 
the center of mass fixed and one gets for the dipole mo- 
ment Di = Zedp = aiNZeR/A. The GT (Goldhaber- 
Teller) model [Gt48] assumes that the restoring force is 
due to the increase of the nuclear surface which leads to 
an extra energy proportional to A'^^^. In this model, the 
inertia is proportional to A and the excitation energy is 
consequently given by oc \/ A'^/^/A = A~^l^ . 



FIG. 35: Hydrodynamical model for collective nuclear vibra- 
tions in halo nuclei. The (a) Steinwedel- Jensen (SJ) mode and 
the (b) Goldhaber-Teller (GT) mode are shown separately. 



For light, weakly-bound nuclei, it is more appropri- 
ate to assume that the neutrons inside the core {Ac, Zc) 
vibrate in phase with the protons. The neutrons and 
protons in the core are tightly bound. An overall dis- 
placement among them requires energies of the order of 
10-20 MeV, well above that of the soft dipole modes. The 
dipole moment becomes Di = eA\(Z ^A — Z A^) j A = 
Z^^j}jedi, where di is a vector connecting the center of 
mass of the two fluids (core and excess neutrons). We 
see that the dipole moment is now smaller than before 
because the effective charge changes from NZ/A in the 

case of the GDR to Z^Jj>f = {Z^A - ZA^) /A. This effec- 
tive charge is zero if ZcA — ZA^ and no pigmy resonance 
is possible in this model, only the usual GDR. 

Figure [35] shows a schematic representation of the hy- 
drodynamical model for collective nuclear vibrations in 
a halo nucleus, as considered here. Part (a) of the figure 
shows the Steinwedel- Jensen (SJ) mode [SJ50] in which 
the total matter density of both the core and the halo 
nucleons do not change locally. Only the local ratio of 
the neutrons and protons changes. Part (b) of the fig- 
ure shows a particular case of the Goldhaber-Teller (GT) 
mode, in which the core as a whole moves with respect 
to the halo nucleons. 

For spherically symmetric densities, the transition den- 
sity in the GT mode can be calculated from 6pp — 
Pp (|r — dpi) — pp (r), where pp is the charge density. Us- 
ing di <C R, it is straight-forward to show that Sp^p'' (r) = 
Sp'^p'' (r) YiQ (f ) , where 



(1) (r\ ^ 



<5<^ (r) 



^'""-Z.^^ay-^, (254) 



and po is the ground state matter density. 

In the Steinwedel- Jensen (SJ) mode, the local varia- 

(2) 

tion of the density of protons is found to be 5pp (r) = 



40 



0.2 




FIG. 36: Contour plot of the nuclear transition density in 
the hydrodynamical model consisting of a mixture of GT and 
SJ vibrations. The darker areas represent the larger values 
of the transition density in a nucleus which has an average 
radius represented by the dashed circle. The legend on the 
right displays the values of the transition density within each 
contour limit. 



(2) 

5p-p (r) Yio (f), where 

^Pp^ (^) = \/^^e//«2if J'l [kr] po (r) 



(255) 



where K = 9.93. If the proton and neutron content of the 
core does not change [SIS90] , the effective charge number 

in the SJ mode is given by Z^fj, = Z^{N-Nc) /A{Z+Nc). 

The transition density at a point r from the center- 
of-mass of the nucleus is a combination of the SJ and 
GT distributions and is given by Spp (r) = Spp (r) Yio (f), 
where 



6pp (r) 




(256) 

Changes can be accommodated in these expressions to 
account for the different radii of the proton and neutron 
densities. 

Figure [36] shows the contour plot, in arbitrary units, 
of the nuclear transition density in the hydrodynamical 
model, consisting of a mixture of GT and SJ vibrations. 
The darker areas represent the larger values of the tran- 
sition density in a nucleus which has an average radius 
represented by the dashed circle. In this particular case, 
a Hartree-Fock density calculation for ^^Li, and a radius 
R — 3.1 fm, was used. The parameters ai and a2 were 



chosen so that Z 



(1) 



-i^fjai = ^gj'^Q!2i i-G. a symmetric mix- 
ture of the SJ and GT modes. 



(2) 
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FIG. 37: Hydrodynamical transition densities for ^^Li and 
three different assumptions for the SJ+GT admixtures, ac- 
cording to Eq. ( |256[ ). The dashed curve is for a GT oscillation 
mode, with the core vibrating against the halo neutrons, with 
effective charge number Z''^^^ — 6/11, radius R — 3.1 fm, and 
ai = 1. The dotted curve is for an SJ oscillation mode, with 
effective charge number Z^j}^ = 2/11, and ai = 1. The solid 
curve is their sum. 



Figure 37 shows the transition densities for ^^Li for 
three different assumptions of the SJ+GT admixtures, 
according to Eq. (2561. The dashed curve is for a GT 



oscillation mode, with the core vibrating against the halo 
neutrons, with effective charge number Z^J'^'^j, = 6/11, ra- 
dius R ^ 3.1 fm, and ai — 1. The dotted curve is 
for an SJ oscillation mode, with effective charge num- 
ber Z^^-'r = 2/11, and = 1. The solid curve is their 



^gyj -= 2/11, and a2 
sum. Notice that the transition densities are peaked at 
the surface, but at a radius smaller than the adopted 
"rms" radius i? = 3.1 fm. 

The liquid drop model predicts an equal admixture of 
SJ-I-GT oscillation modes for large nuclei [Mye77]. The 
contribution of the SJ oscillation mode decreases with 
decreasing mass number, i.e. a2 — > as ^ — *■ 0. This 
is even more probable in the case of halo nuclei, where a 
special type of GT mode (oscillations of the core against 
the halo nucleons) is likely to be dominant. For this 
special collective motion an approach different than that 
used in Ref. [SIS90] has to be considered. The resonance 
energy formula derived by Goldhaber and Teller [GT48] 
changes to 



E 



PR 



\2aRmi\fAr 



1/2 



(257) 



where Ar — Ac {A — Ac) / A and a is the length within 
which the interaction between a neutron and a nucleus 
changes from a zero-value outside the nucleus to a high 
value inside, i.e. a is the size of the nuclear surface, ip is 
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the energy needed to extract one neutron from the proton 
environment. 

Goldhaber and Teller [GT48] argued that in a heavy 
stable nucleus (p is not the binding energy of the nucleus, 
but the part of the potential energy due to the neutron- 
proton interaction. It is proportional to the asymme- 
try energy. In the case of weakly-bound nuclei this pic- 
ture changes and it is more reasonable to associate Lp 
to the separation energy of the valence neutrons, S. I 
will use Lp = pS, with a parameter /? which is expected 
to be of order of one. Since for halo nuclei the product 
aR is proportional to S~^, we obtain the proportionality 
EpR cx S. Using Eq. (|257| for "Li 
i? = 3 fm and p ~ S2n 



with 



- 1 fm, 

0.3 MeV, we get Epr = 1.3 



MeV. Considering that the pygmy resonance will most 
probably decay by particle emission, one gets Ej. ^ 1 
MeV for the kinetic energy of the fragments, which is 
within the right ballpark (see figure 34 1. 

Both the direct dissociation model and the hydrody- 
namical model yield a bump in the response function 
proportional to S, the valence nucleon(s) separation en- 
ergy. In the direct dissociation model the width of the 
response function obviously depends on the separation 
energy. But it also depends on the nature of the model, 
i.e. if it is a two-body model, like the model often adopted 
for ^"'^Be or ^B, or a three-body model, appropriate for 
^^Li and ^He. In the two-body model the phase-space 
depends on energy as p{E) cx cPp/dE (x Ve, while in 
the three-body model p (E) ~ E"^. This explains why the 
peak of figure [34] is pushed toward higher energy values, 
as compared to the prediction of Eq. ( 250 1 . It also ex- 



plains the larger width of dB /dE obtained in three-body 
models. In the case of the pigmy resonance model, this 
question is completely open. 

The hydrodynamical model predicts [Mye77] for the 
width of the collective mode F — hv/R, where v is 
the average velocity of the nucleons inside the nucleus. 
This relation can be derived by assuming that the col- 
lective vibration is damped by the incoherent collisions 
of the nucleons with the walls of the nuclear potential 
well during the vibration cycles (piston model). Using 
V — 3vi?/4, where vp — y/2Ep /niN is the Fermi veloc- 
ity, with Ep = 35 MeV and R = 6 fm, one gets F ~ 6 
MeV. This is the typical energy width a giant dipole res- 
onance state in a heavy nucleus. In the case of neutron- 
rich light nuclei v is not well defined. There are two 
average velocities: one for the nucleons in the core, Vc, 
and another for the nucleons in the skin, or halo, of the 
nucleus, v/j. One is thus tempted to use a substitution in 
the form v — y/W/Wh- Following Ref. [BM93], the width 
of momentum distributions of core fragments in knock- 
out reactions, CTc, is related to the Fermi velocity of halo 
nucleons by Yp — ^/Sof /m . Using this expression with 
Uc ^ 20 MeV/c, we get F = 5 MeV (with i? = 3 fm). 
This value is also not in discordance with experiments 
(see figure pi]) . 

Better microscopic models, e.g. those based on random 
phase approximation (RPA) calculations are necessary to 



study pigmy resonances. The halo nucleons have to be 
treated in an special way to get the response at the right 
energy position, and with approximately the right width. 



D. Comparison between cluster breakup and 
pigmy resonance excitation 

One might argue that the total breakup cross section 
would be a good signature for discerning direct dissoci- 
ation versus the dissociation through the excitation of a 
pigmy collective vibration. The trouble is that the en- 
ergy weighted sum rule for both cases are approximately 
of the same magnitude. This can be shown by using the 
electric dipole strength function in the cluster breakup 
model, namely 



dB{El) Jfie^Z^eff VS;AE- S,, 



,3/2 



dE 



= C 



E^ 



(258) 



where E — Er + Sn is the total excitation energy. C is 
a constant of the order of one, accounting for the correc- 
tions to the wavefunction. 

The sum rule for dipole excitations. Si {El) = 



J- dEE^-qp, is 



Si (El) = C 



2^2 



y2 



(259) 



with Z'^f f 



{Z^A - ZAc) l\AA^ {A ~ A^)\. This is the 
same (with C = 1) as Eq. 1 of Ref. [AGB82], which 
is often quoted as the standard value to which models 
for the nuclear response in the region of pigmy resonance 
should be compared to. The response function in Eq. 
(2581, with C = 1, therefore exhausts 100% of the so- 



called cluster sum rule [AGB82]. The total cross section 
for electron breakup of weakly-bound systems is roughly 
proportional to 5*1. This assertion can be easily verified 
by using eqs. (332 1 and (2581, assuming that the loga- 



rithmic dependence of the virtual photon numbers on the 
energy E = E^ can be factored out of the integral in Eq. 
(I332|. 



The dipole strength of the pigmy dipole resonance is 
given by the same equation ( 259 1 . The constant C is still 



of the order of unity, but not necessarily the same as in 
Eq. ( |258[ ) and the effective charges are also different. For 
the Goldhaber- Teller pigmy dipole model the eff'ective 
charge is given by z'^^}^ = (Z^A — ZAc) /A, whereas for 



the Steinwedel-Jensen it is Z^fj = {Z'^ /A){N - N^) /{Z + 
Nc). Assuming that the Goldhaber- Teller mode prevails, 
one gets the simple prediction for the ratio between the 
cross sections for direct breakup versus excitation of a 
pigmy collective mode: 



direct 



ijpigmy 



c 



A 



Ac (A - A,) ■ 



(260) 



For "Li this ratio is llC/18 while for "Be it is llC/10. 
One thus concludes that it is very difficult to identify 
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FIG. 38: Spectrum of the longitudinal component of the 
center-of-mass velocity of ^Li and two neutrons in the frame 
of the incident ^^Li. (b) Spectrum of the longitudinal compo- 
nent of the relative velocity Vg - Yin- The histogram shows 
the result of a Monte Carlo simulation assuming no Coulomb 
acceleration effects. 



if the nuclear response at low energies is due to collec- 
tive excitations of pigmy resonances, or a trivial phase- 
space effect due to the cluster dissociation, which is a 
basic manifestation of the low binding energy. Other ob- 
servables are necessary to make firm conclusions, such as 
correlations of emitted fragments. 



E. Higher-order effects 

Breakup processes in nucleus-nucleus collisions are 
complicated, in whatever way they are treated. They 
constitute at least a three-body problem, which is further 
complicated due to the long range Coulomb force. Exact 
treatments (such as the Fadeev-approach) are therefore 
prohibitively cumbersome. On the other hand, many ap- 
proximate schemes have been developed and these ap- 
proaches have been used with considerable success. 

Higher-order effects in Coulomb breakup were first ob- 
served in Ref. [Ie93]. In Fig. 38 we see that the veloc- 
ity of ^Li fragments are faster, in average, than those of 
the two detached neutrons in a Coulomb break-up of 28 
MeV/nucleon ^^Li incident on lead targets. 

This effect can be understood qualitatively as follows. 
The ^^Li is deccelerated up to a point where it dissoci- 
ates. This occurs around the distance of closest approach. 
Afterwards the ^Li fragments is accelerated, whereas the 
neutrons are not. Since ^Li is lighter than ^^Li its final 
velocity is greater than the incoming beam velocity. The 
neutrons are consequently slower. 

One should expect that perturbation theory fails in 
describing the breakup process when the cross sections 



attain very high values. In fact, the Coulomb break- 
up probability for weakly-bound nuclei calculated with 
first-order perturbation theory is close to unity. This 
can be understood with use of simple arguments. The 
energy transferred by the Coulomb field to the excitation 
of a projectile nucleus, with N neutrons and Z protons, 
incident with velocity u on a target nucleus with charge 
cZt at an impact parameter b is approximately given by 



E* 



2{NZ/A){ZTe^) 



2\2 



niN 



1)2 y2 



(261) 



where tojv is the nucleon mass. For ^^Be projectiles {N = 
7, Z — A) incident on lead at 6 = 15 fm and v ~ 0.5c, one 
gets iJ* « 1 MeV. This energy is more than sufficient to 
break ^^Be apart, since the separation energy a neutron 
from this nucleus is about 0.5 MeV. This means that, at 
small impact parameters the break-up probability is of 
order of unity and a non-perturbative treatment of the 
break-up process should be carried out. 

A higher-order treatment of Coulomb dissociation is 
therefore desirable. However, one faces the difficulty that 
the final states are in the continuum and the coupling 
matrix elements present divergency problems, caused by 
the non-localized behavior of the continuum wavefunc- 
tions. This difficulty is avoided by a discretization of the 
continuum along the fines discussed below, introduced 
in Ref. [BC92]. The method is a useful starting point 
to understand the main features of a coupled-channels 
calculation with states in the continuum. What I will 
describe next is a simple example of what is commonly 
known as the Continuum-Discretized Coupled- Channels 
(CDCC) method, which can include several levels of so- 
phistication. 



1. Discretization of the continuum 

Our basis of time-dependent discrete states are defined 
as (B s the binding energy) 



-iEot/h I 



, with Eq 



-B 



and ^ e-'""''/^ T^{E) \E,em) (262) 



where \E,£m > are continuum wavefunctions of the pro- 
jectile fragments (without the interaction with the tar- 
get), with good energy and angular momentum quantum 
numbers E, £, m. The functions Tj{E) are assumed to be 
strongly peaked around an energy Ej in the continuum. 
Therefore, the discrete character of the states {(pjim > 
(together with \(f)o >) allows an easy implementation of 
the coupled-states calculations. 

We assume that the projectile has no bound excited 
states. This assumption is often the rule for very loosely- 
bound systems. The orthogonality of the discrete states 
( 262 1 is guaranteed if 



j dE T,{E) Tj{E)=S,. 



(263) 
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For the continuum set \E£m > we use, for the sake of 
simpHcity, the plane wave basis 

<r\Eem> = Ui^E{r)Yern{r) 

o \ 3/4 

which obey the normahzation condition (E — h?q^ /2^) 

< E£m\E'£'m' >= 5u,5mm'5{E ~ E'). (265) 

These states arise from the partial wave expansion of 
the plane wave exp(zq.r). Writing the time-dependent 
Schrodinger equation for ^(t) = X^jfm ^i^™ "^^j^™; 
ing the scalar product with the basis states and using 
orthonormality relations, we get the equations (see Ap- 
pendix B) 



ih 



da 



dt 



Vjim-j'l'm' aj'l'm' 6 
j'£'m' 



-i(E'--Ej)t/h 



(266) 



We use the index j = for the ground state |0 > and 
j = 1, 2, . . . for the discrete continuum states. Vjim-j'e'm' 
are the matrix elements < 4'j£m\V\4>j'iim' >■ 

For Tj{E) we consider two different sets of functions. 
Firstly the set Fi (£;),.. . ,Tm{E); 



T,{E) 
T,iE) 



1 



0. 



for (j -l)a < E < ja, 



otherwise. 



(267) 



This set corresponds to histograms of constant height 
l/v^ and width a. The states Tj{E) trivially satisfy the 
orthonormalization condition of Eq. ( 263 ) . They present 



the advantage of leading to simple analytical expressions 
for the coupling matrix elements. On the other hand they 
have discontinuities at the edges, which lead to numerical 
difficulties. The second set consists of the functions 



-n,{E/a) 



The normalization constant 
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(269) 



guarantees that / Xji^) Xji^) dE — 1. The functions Xj 
are peaked at E — rij a and have width a. The integer 
rij = K.j is proportional to the index-j and the propor- 
tionality constant, a small integer K, is a parameter of 
the set which determines the overlap of two consecutive 
functions Xj ^^nd Xj+i- However, this set fails to satisfy 
the orthogonality condition of Eq. ( 266 1 . This shortcom- 



ing can be fixed by the definition of a new set Tj{E) of 
linear combinations 



N 



T,{E) = Y,C,kXk{E), 



(270) 




-0 06 



o 200 Aoo eoo eoc ooo 

r[fm] 



FIG. 39: Radial wave functions for the discretized continuum 
using the histogram set (a) and the continuous set (b). From 
[BC92]. 



with the coefficients Cij determined so that the resulting 
combinations be orthogonal. These coefficients can be 
found by means of an orthogonalization procedure as, 
e.g., the Gram-Schmidt method. The set of Eq. (2701 



has the advantages of being continuously derivable and 
of leading to reasonably simple coupling matrix elements. 

A comparison between basis states (j^jtmif) generated 
with each of these sets [through Eq. ( 262[ )] is made 
in Fig. [39] One observes that the discrete wavefunc- 
tions 4>jirn decrease rapidly enough with r, so that the 
matrix elements < 4>jim\fYi^j\(l)jiiim' > are finite. The 
use of the histograms (2671 for Tj{E) leads to beats 

These beats are 



39 ». 



in (j)jim as displayed in Fig. 
the result of the discontinuous nature of Tj{E) and 
arise from the interference from the borders of the his- 
tograms. Due to this behavior, the numerical evaluation 
of < 4'jim\'''Yifj_\(f)jieim' > is more involved than with the 
second set of F, -functions, (2701. Indeed, as we see from 



Fig. 39 'b) the beats disappear with the use of the basis 
set ( 270 1 . Although the use of plane- wave basis allows 



the derivation of simple results with both sets, this fact 
is of relevance for improved coupled-channels calculations 
in the continuum. 



fc=i 



Using Eq. (262) and the properties of the spherical 
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harmonics one finds 
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where 



= j r^dr j dETj{E) J dE' Tj,{E') 

X u; E{r)ui,,E'{r). (272) 



From (271 1 one deduces that the interaction potential is 
diS^erent from zero only if |€ — £'| = 1, as expected. 

The use of the plane wave basis is especially useful 
because, exploiting the recursion and closure relations of 
the spherical Bessel functions, one obtains the general 
result 



where 



ti^ (l + £' + 2 



(273) 



Fn' = j dqT,{E)T,,{E) 

Gjr = j dqqT,{E) ^r,,{E) (274) 

with E — h^q^ /2^. Explicit forms can be found for each 
basis set: 

(a) Histogram - Applying this relation to the his- 
togram set ( 267 1 , one can show that for j, / ^ 
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FIG. 40: (a) Radial matrix elements, Eqs. ( |271| l and ([2721 
for the transition j ^ j + 1 (daslied-line) , and for the j — > j 
one (solid line). We used £ = and £' = 1. (b) Radial matrix 
elements for the transition j — > j' , keeping Ej = const, and 
varying [BC92]. 



where F(z) is the gamma-function and we simplified the 
notation using n = nj. For j = 0, or / = 0, one finds 



.^00; jl — Ijl; 00 



E 
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As we have seen above, the use of the plane wave basis 
( 264 1 results in the elegant derivation of Ijg. jigi presented 
by Eqs. (p73| and (p74l). Nonetheless, the s-wave {£ = 0) 



otherwise. 

(275) 

For j = or / = 0, only the integral with £, or = 1 is 
necessary, and the result is 



,3/4 



iOO; jl 
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state of Eq. ( 264 1 is not orthogonal to the bound-state 
wave function. To restore orthogonality one has to add 
an extra piece to this function. One expects however 
that this approximation does not affect the results ap- 
preciably since to access this state one needs at least two 
transitions: the ^ jl followed by the j'l — *■ j'O one. 
But the later transition competes with the transition to 
the ground state, jl — > 00, which is the dominant one. 
A more severe restriction is the use of plane waves to de- 
scribe the continuum. A realistic calculation would have 
to use outgoing waves for u^^^C'") which would carry in- 
formation about the final state interactions of the (5+c)- 
system. 

The break-up probability per unit energy interval, 
P{E)^ is given by 



(276) 



pBU 



where Ej = (j — 1/2) a. 

(h) Continuous basis - For the set of continuous energy 
functions ( 270[ ) one finds, for j, j' ^ 



where 
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(280) 



The validity of the dipole approximation for the in- 
teraction potential ( 271 1 to calculate the continuum- 



continuum coupling can only be justified for <C 1. 
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FIG. 41: (a) Coulomb break-up probability, per unit energy 
interval (MeV~^), of ^^Li projectiles incident on lead at 100 
MeV/nucleon and fe = 15 fm, as a function of the final total 
kinetic energy of the fragments, (b) Breakup probability as a 
function of the time, t = vt/c [BC92]. 



But, as shown in Fig. (39), the discretized wavefunc- 
tions extend up to 400 fm. Thus, unless the matrix ele- 
ments for the continuum-continuum coupling, Eq. (272 1, 
have its main contribution from r ^ 20fm, the dipole ap- 
proximation is not valid. The jj-coupling do satisfy this 
requirement. In this case the wave functions have equal 
energies, but different angular momenta. This causes an 
asymptotically (r :s> 1/q) constant phase difference be- 



This leads 
( |272| for large r. 
coupling. In 



tween the wave functions entering in J, 
to cancellations in the integrand of Eq. 
The situation is different for the 
this case the integrand has contributions from larger val- 
ues of r and these contributions increase with the energy. 
With a correct treatment of the multipole expansion of 
the interaction potential the integrals would de- 

crease with E. We expect that the transitions between 
00 1 and — 1 — > 00 dominate the excita- 

tion process, so that the matrix elements between states 
with j j' do not play an important role. Also, 
to minimize the consequence of the breaking down the 
dipole approximation in the continuum-continuum cou- 
pling at J f , one can use a large parameter K (e.g., 
here K = A was used). This leads to small Ij^j' ■ 



In Figure 41 a) we show the break-up probability per 
unit energy interval for the reaction ^^Li -I- Pb at 100 
MeV/nucleon and 6 = 15 fm, calculated from Eq. (279) 
by solving the coupled-differential Eqs. (266) for aum- 



We see that the energy distribution of the fragments is 
peaked at E ^ 0.3 MeV. Therefore, the most relevant 
momentum transfer to the ^^Li nucleus occurs &t q — 
V2/i6cS/fi- 20 fm-i. 

In Figure 41 'b) the solid line represents P^^ , the total 
breakup probability [Eq. ( 279 1 integrated over energy] , 



as a function of the adimensional parameter r — vt/b, 
for 6 = 15 fm. This is obtained by solving the coupled- 
channels Eqs. (266) for a time t and calculating the sum 



to the neglect of all transitions, except for the - 
ones. The solid-line includes all possible transitions. 



The dashed-line corresponds 
-> ji 
The 



break-up probability occurs in a time scale of At ^ b/v. 
As t ^ oo the break-up probability is 40% smaller than 
that calculated by first order perturbation theory. 



2. Dynamical model in space-time lattice 

Another method for the treatment of higher-order 
effects in Coulomb breakup is the solution of the 
Schrodinger equation in a discretized space-time lattice 
[BB93,BB94,EBB95]. This should yield the same results 
as the CDCC model described in the previous section, if 
the model space is large enough. By that we mean that 
the wavefunction expansion in the CDCC method, and 
the lattice discretization in the present model, are robust, 
convergent, and accurate. With present computer power, 
this can only be achieved with rather simplified structure 
models. 

One assumes a potential model for a cluster-like {c-\-x) 
halo nucleus. For example, a Woods-Saxon potential may 
be used for the nuclear interaction. To this potential a 
time-dependent non-relativistic Coulomb interaction is 
added, 



Vcc 



Kit) 



-R(t)| i?(t)- 



(281) 



-rAj;/Aa, and 



where a = c + x is the projectile, Tc 
r^^ r{l - AJAa). 

The potential above can be expanded into multipoles. 
Normally, the dipole part of the expansion is the most 
relevant one. The time dependent wavefunction for the 
relative motion of c + a; is given by 



^{r,t) = -^u^™(r,t)y^„(f). 



((3.9)) 



Inserting this expansion in the Schrodinger equation 
and retaining only the dipole expansion of the potential 
(p8l]) we get 
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This equation can be solved by a finite difference method. 
A truncation on the m values is needed. Only a few an- 
gular momentum states are needed. Denoting a = {t, to), 
the wave function Ua at time t + At is obtained from 
the wave function at time t, according to the algorithm 
[BB93] 



--A(2) 
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Ua{t). 

(284) 



In this equation r = hAt/A^i,x{Ar)'^ and Sua{t) 
E^i'^Wo'W- Also, 



Va.{r) = VN{r) + 



(285) 



The wave functions Ua{r,t) are discretized in a mesh 
in space, with a mesh-size Ar. The second difference 
operator A^^^ is defined as 



A(2)u« = uii+'Ht) + u'^J^'Ht) - 2u(f')(t), 

with U^^^ EE Ua{rj,t). 



(286) 



The I'' operation on the right side of Eq. (2841 is 



trivial. The second one needs special attention (see Ap- 
pendix H). 

The wave function calculated numerically at a very 
large time will not be influenced by the Coulomb field. 
The numerical integration can be stopped there. The 
continuum part of the wave function is extracted by 
means of the relation (and normalized to unity) 



^'c(r,t) = ^-g, I ^ >] 



1- l<* 



^ >| 
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where ^'gs is the initial wave function. 

This wave function can be projected onto continuum 
eigenstates of the potential V]y{r) in order to obtain the 
excitation probability of that state. For illustration lets 
us use the simplifying assumption that the final states 
are plane waves. A projection onto these states is equiv- 
alent to a Fourier transform of the time dependent wave 
function. One gets 



*c(p) = ^Cf™(p)y(!,„(p) 



(288) 



where 



Cemip) = ^i' Jdrr jdpMtir, t) (289) 

where ji is the spherical Bessel function. The probability 
density for an excitation to a final state with energy E is 
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FIG. 42: (a) Energy spectrum for the breakup of ^^Li at 
28 MeV/nucleon. Data points are from [Iek93]. Solid 
(dashed) line includes (does not include) the reacceleration 
effect [EBB95]. (b) Longitudinal relative momentum of the 
neutrons and the ^Li fragments in the breakup of^^Li projec- 
tiles at 28 MeV/nucleon. Data are from Ref. [Iek93]. Solid 
curve is a calculation from Ref. [EBB95]. 
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In first-order perturbation theory this spectrum would 
be given by 



P<^^\b,E) = 



47r (2ZTe^ 
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where the K' s are the modified Bessel functions and x — 
Eb/hv. In this case dB{E\)/ dE is calculated by using the 
ground state and continuum states of the same potential 
for c-\- X. 

As shown in figure |42] the first-order perturbation the- 
ory does not work well for the break-up of ^^Li at 28 
MeV/nuclcon. 

The reacceleration effect is very important in this case 
due to the very large break-up probabilities at small 
impact parameters. This induces higher order dynam- 
ical continuum-continuum interactions which distort the 
spectrum appreciably. One also observes that although 
the dynamical corrections modify the spectrum in the 
right direction the discrepancy with the experiment is 
appreciable. This might be a deficiency of the cluster 
model and has to be studied more carefully. Also, the 
spectrum of the relative longitudinal momentum of the 
fragments is not in a good agreement with the dynamical 
calculations of Ref. [EBB95] , at least for the higher part 
of the spectrum. This is shown in the figure 42 ^b). 



Higher-order effects are manifest in published theoret- 
ical calculations, and seem to have some experimental 
support. However, in most cases the effect is not rele- 
vant. There are also some theoretical results which does 
not show post-acceleration effects. 
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X. RADIATIVE CAPTURE IN STARS 

If you have read to this point, you are now convinced 
that there must be numerous situations where Coulomb 
excitation and dissociation can be used for accessing pre- 
cious information on nuclear structure. Obviously, any 
nucleus is part of a process in stars (perhaps the super- 
heavies are not). Nuclear physics is therefore ultimately 
linked to astrophysics. We also expect that Coulomb 
dissociation plays a role in stars. In fact, it also plays a 
role in Cosmology. The Greisen-Zatsepin-Kuzmin (GZK) 
effect, for example, is nothing more that Coulomb exci- 
tation by (3 K) real photons remnant from the Big Bang 
and by virtual photons in intergalactic magnetic fields 
[GZK66]. 

As we have seen before. Coulomb excitation probes the 
same matrix elements as real photons. It is thus natural 
that many astrophysical reactions in stars involving pho- 
tons in the initial or final channels can be studied on earth 
via Coulomb excitation. In particular. Coulomb dissoci- 
ation is directly related to radiative capture in stars. In 
the next sections we discuss the Coulomb dissociation 
method, its experimental applications to some astrophys- 
ical problems and some of the newest theoretical devel- 
opments and difficulties associated to it. 



A. The Coulomb dissociation method 

Coulomb dissociation is a process analogous to what 
happened to the comet Shoemaker-Levy as it disinte- 
grated during its approximation to Jupiter in 1994 (see 
Figure 43 1. Approximately 1.5 to 2.2 hours after clos- 
est approach, the comet (which was presumably a sin- 
gle body at the time) was broken apart by tidal forces 
into at least 21 pieces. The pieces continued to orbit 
Jupiter with a period of approximately 2 years. Due to 
gravitational forces from the Sun, which changed the or- 
bits slightly on the next approach to Jupiter, the pieces 
impacted the planet. Much stronger tidal forces occur 
when nuclei come close to each other due to their mutual 
electromagnetic field. This leads to their dissociation, 
especially for weakly-bound nuclei. 

The idea behind the Coulomb dissociation method is 
relatively simple. The (differential, or angle-integrated) 
Coulomb breakup cross section for a + A b + c + A 
follows from 



N-^{E^-e-4') .<J:,Xa^b+c{E^), (293) 



where is the energy transferred from the relative mo- 
tion to the breakup, and cr^+o ^ b+ci^i) photo- 
dissociation cross section for the multipolarity ttL and 
photon energy E^. Time reversal allows one to deduce 




FIG. 43: Jupiter and comet Shoemaker-Levy 9, as imaged 
by the Hubble Space Telescope (HST), on May 18, 1994, 
when the giant planet was at a distance of 420 million miles 
(670 million km) from Earth. The gravitational interaction 
of Jupiter with the comet has broken it up into many pieces 
(picture from NASA). Much stronger tidal forces occur when 
nuclei come close to each other due to their mutual electro- 
magnetic fields. 



(E^), i.e. 



2(2j„ + 1) 



7, 



(2jb + l)(2je+l)fc2 



(294) 



the radiative capture cross section b -\- c 



7 from 



where k-y — E^/hc is the photon wavenumber, and 
k = y/2iJ,{Ey — B)/h is the wavenumber for the rela- 
tive motion of b-t-c. Except for the extreme case very 
close to the threshold (A: — > 0), we have <C k, so that 
the phase space favors the photodisintegration cross sec- 
tion as compared to the radiative capture. Direct mea- 
surements of the photodisintegration near the break-up 
threshold do hardly provide experimental advantages and 
seem presently impracticable. On the other hand the co- 
pious source of virtual photons acting on a fast charged 
nuclear projectile when passing the Coulomb field of a 
(large Z) nucleus offers a more promising way to study 
the photodisintegration process as Coulomb dissociation. 

This method was introduced in Ref. [BBR86] and has 
been tested successfully in a number of reactions of in- 
terest to astrophysics. The most celebrated case is the 
reaction ''Be(p,7)*B, first studied in Ref. [Mot94], fol- 
lowed by numerous experiments in the last decade. This 
reaction is important because it produces ^B in the core 
of our sun. These nuclei decay by emitting a high energy 
neutrinos which are one of the best probes of the sun's in- 
terior. The measurement of such neutrinos is very useful 
to test our theoretical solar models. 
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FIG. 44: Electric dipole response function for ^He. The 
shaded area represents the experimental results from a 
Coulomb dissociation experiment [Aum99]. The dashed and 
dotted lines correspond to results from three-body decay mod- 
els from Refs. [Da98,Cob97]. (Courtesy of T. Aumann). 



Another example is the two-neutron capture on ^He 
could perhaps play a role in the post-collapse phase in 
type-II supernovae. The bottleneck in this nucleosyn- 
thesis scenario is the formation of nuclei with A > 9 
from nucleons and a-particles. In principle, the reaction 
■*He(2rt, 7)^He could be relevant in bridging the instabil- 
ity gap at A = 5, although it is believed that this reac- 
tion cannot compete with the (an, 7) process in a type-II 
supernova scenario. Experiments with Coulomb dissoci- 
ation have been used to study this question, as shown 
in the example presented in Figure |44j The figure dis- 
plays the electric dipole response function for ^He. The 
shaded areas represent the experimental results from a 
Coulomb dissociation experiment [Aum99]. The dashed 
and dotted lines correspond to results from three-body 
decay models from Refs. [Dan98,Cob97]. 

In figure |45] we show the measurement of two-body 
correlations in the three-body decay of ^He. The lower 
panels display the ratio between the measured a-n and n- 
n relative-energy spectra (upper panels) and the spectra 
simulated (histograms) according to standard phasespace 
distributions [Aum99]. From the analysis of this exper- 
iment it was found that 10% of the dissociation cross 
section proceeds via the formation of ^He. A rough es- 
timate yields 1.6 mb MeV for the photoabsorption cross 
section for ^He(7,n)^He, which agrees with theoretical 
calculations [Efr96] . From this experiment one concludes 
that the cross sections for formation of ^He and ^He via 
one (two) neutron capture by '^He are not large enough 
to compete with the (an, 7) capture process (for more 
details, see Ref. [Aum06]). Nonetheless, this and the pre- 
viously mentioned examples, show the relevance of the 
Coulomb dissociation method to assess some of the basic 
questions of relevance for nuclear astrophysics. 
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FIG. 45: Measurement of two-body correlations in the three- 
body decay of ''He. The lower panels show the ratio between 
the measured a — n and n-n relative-energy spectra (upper 
panels) and the spectra simulated (histograms) according to 
standard phasespace distributions [Aum99]. (Courtesy of T. 
Aumann). 

B. Higher-order corrections 

The Coulomb dissociation method is specially useful if 
first-order perturbation theory is valid. If not, one can 
still extract the electromagnetic matrix elements involved 
in radiative capture reactions. But a much more careful 
analysis of the high-order effects need to be done (see 
figure |46[ We will discuss briefly a continuum-discretized 
coupled-channels calculation, appropriate for relativistic 
collision energies. 

Coupled-channels calculations are a complication that 
we could live better without. It helps some people to 
get jobs, but they do not make life easier, as they have 
very often unphysical (or better said, unpredictable) re- 
sults. The reason for the uncertain behavior is that, in 
each interaction step, plus and minus signs of matrix el- 
ements accumulate in one or the other direction (like a 
brownian motion). In seminars you might ask why a 
coupled-channels calculation reduced, or increased, the 
cross section. The plain answer is invariably: "it is a 
coupled-channels effect" . Well, who said that hfe is easy? 

1. Continuum-Discretized Coupled- Channels 

The continuum-discretized coupled-channels method 
(CDCC) [Kam87] is one of the most accurate mod- 
els to describe the breakup of halo nuclei taking ac- 
count of higher-order couplings explicitly. The eikonal 
CDCC method (E-CDCC) [Oga03,Oga06], is a deriva- 
tion of CDCC that enables one to efliciently treat the 
nuclear and Coulomb breakup reactions at Eiab > 50 
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MeV/nucleon. An essential prescription described in 
Refs. [Oga03,Oga06] is the construction of hybrid (quan- 
tum and eikonal) scattering amplitudes, with which one 
can make quantum-mechanical (QM) corrections to the 
pure eikonal wavefunctions with a minimum task. These 
corrections are, however, expected to become less impor- 
tant as the incident energy increases. 

First I will describe a simplified CDCC method pre- 
sented in Ref. Ber05] where the inclusion of relativis- 
tic continuum-continuum effects in CDCC calculations 
were introduced. Let us consider the Klein-Gordon (KG) 
equation with a potential Vq which transforms as the 
time-like component of a four- vector [AC79]. For a sys- 
tem with total energy E (including the rest mass M) , the 
KG equation can be cast into the form of a Schrodinger 
equation (with h = c = 1), 



(V^ + ~ U) ^ = 0, 



(295) 



where fc^ = {E^ - M^) and U = 2Vo{2E - Vq). When 
Vo < M, and E ~ M, one gets U = 2MVo, as in the 
non-relativistic case. The condition Vq <C M is met in 
peripheral collisions between nuclei at all collision ener- 
gies. Thus, one can always write U = 2EVo. A further 
simplification is to assume that the center of mass mo- 
tion of the incoming projectile and outgoing fragments is 
only weakly modulated by the potential Vq. 

To get the dynamical equations, one discretizes the 
wavefunction in terms of the longitudinal center-of-mass 
momentum fc^, using the ansatz 



(296) 



In this equation, {z,h) is the projectile's center-of-mass 
coordinate, with b equal to the impact parameter. (^) 
is the projectile intrinsic wavefunction and (fc, K) is the 
projectile's center-of mass momentum with longitudinal 
momentum k and transverse momentum K. There are 
hidden, uncomfortable, assumptions in Eq. (2961. The 



center of mass of a relativistic system of interacting par- 
ticles is not a well defined quantity. Also, the separation 
between the center of mass and intrinsic coordinates is 
not permissible under strict relativistic treatments. For 
high energy collisions we can at best justify Eq. (296) 



for the scattering of light projectiles on heavy targets. 
Eq. ( 296 1 is only reasonable if the projectile and target 



closely maintain their integrity during the collision, as in 
the case of very peripheral collisions. 

Neglecting the internal structure means (j>ka (^) = 1 
and the sum in Eq. ( 296 1 reduces to a single term with 
a = 0, the projectile remaining in its ground-state. It is 



straightforward to show that inserting Eq. ( 296 ) in the 
KG equation 



2EVo) * = 0, 



(297) 



and neglecting W'^So {z,h) relative to ikdzSo {z,h) one 
gets ikdzSo {z, b) = EVqSo {z, b), which leads to the cen- 




FIG. 46: Schematic description of a multi-step excitation, 
with excmsions in the continuum. 



ter of mass scattering solution 



So (z,b) = exp 



-1 /" dz'Vo{z\h) 

'J — OO 



(298) 



with V — k/E. Using this result in the Lippmann- 
Schwinger equation, one gets the familiar result for the 
eikonal elastic scattering amplitude, i.e. 

k f 

/o = -i— / dh exp {iQ ■ b) {exp [ix{h)] - 1} , 

where the eikonal phase is given by x(t>) = '5o(oo,b), 
and Q = K' — K is the transverse momentum transfer. 
Therefore, the elastic scattering amplitude in the eikonal 
approximation has the same form as that derived from 
the Schrodinger equation in the non-relativistic case. 
For inelastic collisions we insert Eq. ( 296 1 in the KG 



equation and use the orthogonality of the intrinsic wave- 
functions (f>k^ (^). This leads to a set of coupled-channels 
equations for Sa- 



k^) Sc^e' 



with the notation la) 



^(a|C/|a') e"=='^ (299) 



Neglecting terms of the 



form V^5q {z,h) relative to ikdzSa {z,h), Eq. (2991 re- 
duces to 



(300) 
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The scattering amplitude for the transition 
given by 



using 



/a(Q) 



ik 
'2^ 



with Sg (b) = Sa {z 



dh exp (iQ • b) [Sa (b) - 6^^] , 

(301) 

= 00, b). The set of equations (300) 



and ( [30T| are the eikonal-CDCC equations (E-CDCC) 
They are much simpler to solve than the low-energy 
CDCC equations because the z and b coordinates de- 
couple and only the evolution on the z coordinate needs 
to be treated non-perturbatively. 

The E-CDCC equations can be used to study the dis- 
sociation of projectiles at high energies. The ener- 
gies transferred to the projectile are small, so that the 
wavefunctions can be treated non-relativistically in the 
projectile frame of reference. In this frame the wavefunc- 
tions will be described in spherical coordinates, i.e. \a) = 
\jlJM) , where j, I, J and M denote the angular mo- 
mentum numbers characterizing the projectile state. It 
is important to notice that the matrix element (a \ Vs \ a') 
is Lorentz invariant. Boosting a volume element from the 
projectile to the laboratory frame means (fi^ — > 

where 7 = (l — ii^) is the Lorentz contraction fac- 
tor. The intrinsic projectile wavefunction is a scalar and 
transforms according to <pa i^x,^y,^z) ~* (t>a (C^^ , , 76 ) > 
while Vq, being the time- like component of a four- vector, 
transforms as Vb (&, z; ^x, Cy, 6) ^ jVo {b, z; ^x, ^y,'y£,z)- 
Thus, redefining the integration variable z in the labora- 
tory as = leads to the afore mentioned invariance. 
We can therefore calculate (a | Vol a') in the projectile 
frame. 

The longitudinal wavenumber ka — {E^ — M^)^/^ also 
defines how much energy is gone into projectile exci- 
tation, since for small energy and momentum transfers 
fc'„ - fc„ = {E'^ - Eg) /v. In this limit, eqs. dSOO]) and 



(301 1 reduce to semiclassical coupled-channels equations, 
if one uses z = vt for a projectile moving along a straight- 
line classical trajectory, and changing to the notation 
Sa {z,b) — aa{t,b), where aa(t,b) is the time-dependent 
excitation amplitude for a collision wit impact param- 
eter b. With these changes, equation (301) also agrees 



with the equivalent semiclassical equation for the inelas- 
tic scattering amplitude. Therefore, for non-relativistic 
collisions and small momentum and energy transfers, the 
above derivation reduces to previous methods used in the 
literature. 

If the state \a) is in the continuum (positive 
proton-|-^Be energy) the wavefunction is discretized ac- 
cording to 



\a;Ea)^ J dE'^T{E'^) \a;E'^\ 



(302) 



where the functions T{Ea) are assumed to be strongly 
peaked around the energy Ea with width AE. For con- 
venience the histogram set is chosen. The inelastic cross 
section is obtained by solving the E-CDCC equations and 



da 



dQ.dE a 



(Q)r ^'\Ea). 



(303) 



The potential Vo contains contributions from the nu- 
clear and the Coulomb interaction. The nuclear po- 
tentials are constructed along traditional lines of non- 
relativistic theory. The potentials are also expanded into 
multipolarities. These potentials are then transformed 
as the time- like component of a four- vector, as described 
above. The multipole expansion of the Coulomb inter- 
action in the projectile frame including retardation. The 
first term (monopole) of the expansion is the retarded 
Lienard-Wiechert potential which does not contribute to 
the excitation, but to the center of mass scattering. Due 
to its long range, it is hopeless to solve Eq. (300) with 



the Coulomb monopole potential, as Sa {z, b) will always 
diverge. This can be rectified by using the regularization 
scheme 

Sa (b) ^ exp [i (2ryln {kb) + Xa (&))] (b) , (304) 

where ?/ = ZpZ^ j^v and the Sa (b) on the right-hand 
side is now calculated without inclusion of the Coulomb 
monopole potential in Eq. ( 300 1 . The purely imaginary 



absorption phase, Xa {b), was introduced to account for 
absorption at small impact parameters. 

The multipole-expansion of the relativistic Coulomb 
potential between the target nucleus (T) with the atomic 
number Zt and the projectile (P), consisting of C and v 
clusters, are given in Ref. [Ber05]: 



^El.-^f {i) (^2^^2,2)3/2 ^ (if M = 0) 



'yZTCCEi f Tb (if fJ- — ±1) 

M = 0) 
(305) 



for the El (electric dipole) field and 



V^E2^i = 



7ZTeeE2 



(62 + -^2,2)5/2 

b^ (if ^ = ±2) 
T^-f^bz (if Ai-±1) (306) 
/273 (27^22 - 62) (if ^^0) 



for t he E2 (electric quadrupole) field. In Eqs. (305) 



(3061, = [Z^iAc/Ap)^ + Zci-A,/Apy]e are ef- 



fective charges for A = 1 and 2 multipolarities for the 
breakup of P — > C-|-v. The intrinsic coordinate of v with 
respect to C is denoted by ^ and 6 is the impact parame- 
ter (or transverse coordinate) in the collision of P and T, 
which is defined by 6 = \/ -V with R = (x, y, z), the 
relative coordinate of P from T in the Cartesian repre- 
sentation. Note that these relations are obtained with so- 
called far-field approximation [EB05], i.e. i? is assumed 
to be always larger than i^. The Coulomb coupling po- 
tentials in E-CDCC are obtained with Eqs. ([305|-p06l 
as shown below. 
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The full E-CDCC equations, including relativistic ef- 
fects and recoil corrections for the three-body reaction 
are given by [OB09]: 



E. 



X V^?)(z) e'(^='"^=)^ (307) 



where c denotes the channel indices {i, £, m}; i > 
(i = 0) stands for the zth discretized-continuum (ground) 
state and £ and m are respectively the orbital angular mo- 
mentum between the constituents (C and v) of the pro- 
jectile and its projection on the z-axis taken to be parallel 
to the incident beam. Note that we neglect the internal 
spins of C and v for simplicity. The impact parameter 
b is relegated to a superscript since it is not a dynami- 
cal variable. The total energy and the asymptotic wave 
number of P are denoted by Ec and Kc, respectively, and 



(308) 



with the Sommerfeld parameter. The local wave num- 
ber Kc^\z) of P is defined by energy conservation as 



Er 



R 



(309) 



where mp is the mass of P and Zpe (^tc) is the charge 
of P (T). The reduced coupling-potential 5'|.^/(z) is given 

by 



R 



where 



(«'c|C/cT- 



-i{m — m) 



(310) 



(311) 



The <& denotes the internal wave functions of P, (pu is 
the azimuthal angle of b and Uct (Uvt) is the potential 
between C (v) and T consisting of nuclear and Coulomb 
parts. In actual calculations we use the multipole expan- 
sion J^'^){z) = ^cc''\^)^ explicit form of which is 
shown in Ref. [Oga06]. 

In order to include the dynamical relativistic effects 
described above, we make the replacement 



(312) 



The factor /a,^ is set to unity for nuclear couplings, while 
for Coulomb couplings we take 




(A= 1,A* = 0) 
{\ = 2,^l = ±l) 
(otherwise) 



(313) 



following Eqs. (305 1 and (306 1. Correspondingly, we use 



R 



1 



+ {-fzY 



(314) 



in Eqs. (309) and (3101. The Lorentz contraction factor 
7 may have channel-dependence, i.e., 7 = Ec/{mpc'^), 
which we approximate by the value in the incident chan- 
nel, i.e., Eo/{mpc^). 

It should be remarked that we neglect the recoil motion 
of T in Eq. (307 1; this can be justified because we con- 



sider reactions in which T is significantly heavier than P 
and we only treat forward-angle scattering in the present 
study [Ber05], as shown below. Note also that in the high 

incident-energy limit 7?.[.^^(z) 1 and Kc''\z) Kc, un- 
less the energy transfer is extremely large. Thus, in this 
limit Eq. (307) becomes Lorentz covariant, as desired. 

Using Eqs. (307l-( [314 l, the dissociation observables in 
reactions of loosely bound nuclei and ^^Be on 208p^ 
targets was calculated [OB09]. 
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FIG. 47: The total breakup cross section for ^B+'^°^Pb at 
250 MeV/nucleon, as a function of the scattering angle of the 
cm. of the projectile after breakup. The solid and dashed 
lines show the results of the full CC calculation with and 
without the dynamical relativistic effects, respectively. 



Figure 47 shows the total breakup cross section of °B 
by 208p|3 250 MeV/nucleon, as a function of the scat- 
tering angle 9 of the center-of-mass (cm.) of the projec- 
tile after breakup. The solid and dashed lines represent 
the results of the E-CDCC calculation with and with- 
out the dynamical relativistic effects, respectively; in the 
latter we set 7=1 instead of the proper value 1.268 in 
Eqs. (312|-(314). Note that in all calculations shown in 
this work we use relativistic kinematics. One sees that 
the dynamical relativistic correction gives significantly 
larger breakup cross sections for 6 < 1.3 degrees; the dif- 
ference between the two around the peak is sizable, i.e. 
of the order of 10-15 %. 

We show in Fig. |48] how the coupled-channel calcula- 
tions affect the role of the dynamical relativistic correc- 
tion. The left (right) panel corresponds to the calculation 
with both nuclear and Coulomb breakup (only Coulomb 
breakup) . 

In each panel the solid (dotted) and dashed (dash- 



52 



20 



' ^ (a) nuclear and Coulomb breakup 

_ y fullCCrel. 

' ' • full CC nonrel. 

st-order rel. 
1 st-order nonrel. 




e (deg) 



20 



(b) Coulomb breakup only 

full CC rel. 

full CC nonrel. 

1 st-order rel. 

\ 1 st-order nonrel. 




e (deg) 

FIG. 48: The total breakup cross sections for *B+^°^Pb at 250 
MeV/nucleon with nuclear and Coulomb breakup (top panel) 
and only Coulomb breakup (bottom panel). The solid (dot- 
ted) and dashed (dash-dotted) lines show the results of the 
full CC (hrst-order perturbative) calculation with and with- 
out the relativistic correction, respectively. 



dotted) lines show the results of the full CC (first-order 
perturbative) calculation with and without the dynam- 
ical relativistic correction, respectively. One sees from 
the left panel that relativistic corrections are slightly 
quenched when the first-order calculation is adopted. 
This is also the case when only Coulomb breakup is in- 
cluded, as shown in the right panel of Fig. |48] We stress 
here that due to continuum-continuum couplings rela- 
tivistic effects are non-linear (and non-trivial to inter- 
pret). Thus, as shown in Fig. |48] one cannot infer the 
effect of relativistic corrections by simply carrying out 
first-order calculations. Moreover, it should be remarked 
that the full CC and first-order calculations give quite 
different breakup cross sections. 

One concludes that the effects of relativistic correc- 



tions of the nuclear and Coulomb coupling potentials 
on the breakup cross sections of the weakly bound pro- 
jectiles and "Be by ^ospb targets at 250 and 100 
MeV/nucleon. The relativistic corrections modify ap- 
preciably the breakup cross sections, at the level of 15 
% (10%), in collisions at 250 (100) McV/nuclcon. This 
change is found to be due mainly to the modification of 
the Coulomb potential. We have also shown that rela- 
tivistic corrections infiuence continuum-continuum tran- 
sitions. 



XI. WHAT YOU WANTED TO KNOW, BUT 
WERE AFRAID TO ASK 



A. Close and far fields 

So far, in all derivations of Coulomb excitation we re- 
stricted to distant collision, where the intrinsic radial 
coordinates are smaller than the minimum distance be- 
tween the colliding nuclei. In calculations of the Coulomb 
dissociation of proton halo nuclei, it is of interest also 
to consider close collisions, since the density of the va- 
lence proton can extend to very large distances. More- 
over, higher-order processes may also play a role, as 
suggested by the nonrelativistic calculations of the ^B 
breakup. This may be of relevance for the calculation of 
continuum-continuum matrix elements. 

We derive a general expression for the multipole ex- 
pansion of the electro-magnetic interaction from a fast 
charge particle, which can be employed in higher-order 
dynamical calculations of the Coulomb excitation of halo 
nuclei [EB05]. 

The relativistic expression for the electro-magnetic in- 
teraction from a fast charged particle, moving on a 
straight line trajectory in the z-direction and impact pa- 
rameter b with respect to a target nucleus, i.e. r' = 
(6, 0, vt), and acting on a proton with coordinates r is 

Km(r,r') = ZTe2((/.- ^[p<^+0p]), (315) 

where (/> is the Lienard-Wiechert potential 
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^|b-b'|2 + 72(2_^/)2' 



(316) 



and p is the momentum operator associated with r. 

We start out with the Fourier representation of cj) which 
can be written as (see appendix of Rcf [BB88]) 



47r 
(27r)3 

47r 
(2^ 



g'iq-(r-r') 




giq-(r-r') 



1 



/32 C0S2( 



(317) 



where (3 = v/c. Except for the last factor, this is just the 
ordinary Coulomb potential. 
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To proceed, we introduce the multipole expansion of 
the last term in Eq. pl7| (see Eq. 8.825 of Ref. [GR80]) 



1 



RL^L+A{a,b) = 0, for a>b, and A = 2,4,... (327) 
The exphcit, non-zero expression is [EB05] 



1-/32 cos2(6') 



= ^ (2A+l)r'gA(r')PA(cos(0)), 



(318) 



A even 

and insert the plane wave expansion 

e^q- = 47r^ ,uiqr) Yl^MYLM{f) (319) 

LM 



R 



L,L+A 



{a,b) 



2L + 1 + A 



_p(^+iA-i)(l„2(a/6)2), 



for the two plane waves in Eq. (|317]). Thus we obtain 

</, = 4^ ^ RLL'ir, r') ^ Al^mW >LM(r) Yl,j,,{f'), 

LL' m 

(320) 

where 

RLL'{r,r') = - / dq jL{qr) jL'iqr') 

^'^ JL+i/2{qr) JL'+i/2{qr'), 



(328) 

which is valid when a < b, and A = 0, 2,4, .... 

It is convenient to exploit the properties ( 326|327 l and 
write expressions for that are valid for distant and close 
collisions, i. e. r < r' and r > r', respectively. The 
expression one obtains for distant collisions is 



A'dist 



and 



'rr' Jo 9 



(321) 



(322) 



= J2 (2A+1) p'^QxiP'^) {L'MXO\LM) (LOAO|i'0). 

A even 

(323) 

The last two Clebsch-Gordan coefficients is the matrix 
element {YLAi\Pxicos{9))\YL' m) ■ 

We note that different multipoles, (LM) and (L'M), 
are mixed in the expansion ( 320 ) . Let use therefore first 



take a loolc at the contribution from diagonal terms with 
L' = L. Here the integral (321 1 is particularly simple 



RLL{r,r') 



1 



2i + l ri'+i 



(324) 



where r< = min(r, r') and r> = max(r, r'). Inserting 
this into Eq. ( |320| we obtain the diagonal contribution 

L > m 

(325) 

This has a structure similar to the non-relativistic ex- 
pression, except for the factor Allm{P)- In fact, in the 
non-relativistic limit we have All'a/(/3 = 0) = S^l' and 
the non-relativistic expression is recovered. 

The off-diagonal (i. e. L' ^ L) contributions to 
are more complicated. According to Eq. (323), and 



the symmetry in L and L', we only need expressions for 
i' = L A, A = 2,4, ... One finds that [EB05] 



RL.L+A{a,b) = RL+Ax{b,a) 



(326) 



A=0,2,. 

X Al,L+A,M 5^l+a,j\/(^')> 



LM 



(329) 



whereas the contribution for close collisions is 



t^closc 



A=0.2,. 



LM 



L-A.L,m Y]^_jy T 



(330) 



The last sum over A is finite since L — A must be non- 



negative. The two expressions (329) and (330) reduce to 
for A 



Eq. (325 



The a 



0. 



ibove results have been recently used in Ref. 
[OB09b] to study the Coulomb breakup of ^B. It was con- 
cluded that the close field contributes very little to the 
breakup cross section. As there exists very few proton 
halo nuclei as loosely-bound as ®B, we conclude that the 
usual multipole expansion with distant fields describes 
Coulomb excitation satisfactorily. 



B. Comparison to electron scattering 

In the plane wave Born approximation (PWBA) the 
cross section for inelastic electron scattering is given by 
[Ba62,EG88] 

da _ Sttc^ fp'\\^j + '^^P ■ P' + "^^^"^ 



{hcY \p 



X \F,j{q;CL)\' 



EE' - c2 (p • q) (p' • q) - m'^c^ 
(9^ - ^o)^ 



\F,j{q;ML)\^+\F,,{q;EL)\^ 



(331) 



where Ji {Jf) is the initial (final) angular momen- 
tum of the nucleus, {E,\>) and (£^',p') are the initial 
and final energy and momentum of the electron, and 

(E-E') (p-pQ 



(90, q) 



he ' h 



is the energy and momen- 



tum transfer in the reaction. Fij (q; HL) are form factors 
for momentum transfer q and for Coulomb (C), electric 
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{E) and magnetic (M) multipolarities, H = C, E, M, re- 
spectively. 

Here we will only treat electric multipole transitions. 
Moreover, we will treat low energy excitations such that 
E, E' ^ hcqa, which is a good approximation for electron 
energies E ~ 500 MeV and small excitation energies 
AE = hcqo ~ 1 — 10 MeV. These are typical values 
involved in the dissociation of nuclei far from the stability 
line. 

Using the Siegert's theorem [Sie37,Sa51], one can show 
that the Coulomb and electric form factors in eq. (331) 
are proportional to each other. Moreover, for very for- 
ward scattering angles {9 <C 1) the PWBA cross section, 
eq. 331 can be rewritten as 



da 



dndE^ 



E 

L 



dN, 



(EL) 



{E, Ej 



dndE^ 



^) {E^) , (332) 
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where ct^^^' (^7)1 with E~^ = hcqo, is the photo- nuclear 



cross section for the iJL-multipolarity, given by eq. ( 103 1 
which we rewrite as [BB88] 



(27r)^(i + l) ( E. 
L[(2L+l)!f 



7 

he 



2L-1 



dB {EL) 
dE^ 



(333) 

In the long- wavelength approximation, the response func- 
tion, dB{EL) /dEj, in eq. p33t is given by 



dB (EL) 
dE^ 



\{Jf\\YumJ^: 
2 J, -I- 1 



drr^+^ Splf^'>{r) 



(334) 



where w {E^) is the density of final states (for nuclear 
excitations into the continuum) with energy Ej = Ef — 

Ei. The the transition density 5p\^ (r) will depend 
upon the nuclear model adopted. 



For L>1 one obtains from Eq. (331 1 that 



dN, 



(EL) 



{E,E^,( 



dndE-y 



AL a 
L+IE 



cos2 (6'/2)sin"3(6'/2) 
1 + (2£;/M^c2) sin^ {0/2) 
2 



2E 

E~ 



L 

T T sm , „ 

L + l V 2 



2E . 

sm , 

E^ \2 



tan 



2L-1 



(335) 



One can also define a differential cross section inte- 



FIG. 49: Virtual photon spectrum for the El, E2 and E3 
multipolarities in electron scattering off arbitrary nuclei at 
Ee — 100 MeV and maximum scattering angle of 5 degrees. 



Eqs. ( 332|335 1 show that under the conditions of 
the proposed electron-ion colliders, electron scattering 
will offer the same information as excitations induced 
by real photons. The reaction dynamics information is 
contained in the virtual photon spectrum of eq. (3351, 



while the nuclear response dynamics information will be 
contained in Eq. (3341. This is akin to a method devel- 



oped long time ago by Fermi [Fe24] and usually known as 
the Weizsaecker- Williams method [We34] (see section on 

historical note. The quantities dN^^^^ / dU,dE^ can be 
interpreted as the number of equivalent (real) photons 
incident on the nucleus per unit scattering angle and 
per unit photon energy E^. Note that £'0 (monopole) 
transitions do not appear in this formalism. As imme- 
diately inferred from Eq. (334 1, for L ~ Q the response 



function dB{EL) /dEj vanishes because the volume in- 
tegral of the transition density also vanishes in the long- 
wavelength approximation. But for larger scattering an- 
gles the Coulomb multipole matrix elements (CL) in Eq. 
(3311 are in general larger than the electric {EL) multi- 



poles, and monopole transitions become relevant [Sch54]. 



(EL') ( E'2.') 

grated over angles. Since a\ does not depend on the spectrum yields the ratios dNe /dN, 



In figure 49 we show the virtual photon spectrum for 
the El, E2 and E3 multipolarities for electron scatter- 
ing off arbitrary nuclei at E,, = 100 MeV. These spectra 
have been obtained assuming a maximum scattering an- 
gle of 5 degrees. An evident feature deduced from this 
figure is that the spectrum increases rapidly with decreas- 
ing energies. Also, at excitation energies of 1 MeV, the 

500 and 



(El) 



scattering angle, this can be obtained from eq. (335) 



by integrating dNe /dQ.dE^ over angles, from 6',nin = 



dNi^^^ /dNi^^^ 



100. However, although dN^J^^^ /dE^ 
increases with the multipolarity L, the nuclear response 
decreases rapidly with L, and El excitations tend to 
dominate the reaction. For larger electron energies the 



E^/ E to a maximum value 9^ 
experimental setup. 



which depends upon the 



55 



ratios TV^^^) /Ar(£^i) and (^3) /TV^^^) decrease rapidly. 

In figure |50] we show a comparison between the El 
virtual photon spectrum, dNg/dEj, of 1 GeV electrons 
with the spectrum generated by 1 GeV/nucleon heavy 
ion projectiles. In the case of Coulomb excitation, the 
virtual photon spectrum was calculated according to Eq. 
( |136[ ). For simplicity, we use for the strong interaction 
distance i? = 10 fm. The spectrum for the heavy ion 
case is much larger than that of the electron for large 
projectile charges. For ^ospj-, pj-ojectiles it can be of the 
order of 1000 times larger than that of an electron of the 
same energy. As a natural consequence, reaction rates 
for Coulomb excitation are larger than for electron ex- 
citation. But electrons have the advantage of being a 
cleaner electromagnetic probe, while Coulomb excitation 
at high energies needs a detailed theoretical analysis of 
the data due to contamination by nuclear excitation. As 
one observes in figure |50| the virtual spectrum for the 
electron contains more hard photons, i.e. the spectrum 
decreases slower with photon energy than the heavy ion 
photon spectrum. This is because, in both situations, 
the rate at which the spectrum decreases depends on 
the ratio of the projectile kinetic energy to its rest mass, 
E/mc^, which is much larger for the electron (m = rrie) 
than for the heavy ion (m — nuclear mass). 

To obtain an effective luminosity per unit energy, the 
equivalent photon number is multiplied by the experi- 
mental luminosity, LeAi i-e. '^^e^ ~ ^^aJ^- The 
number of events per unit time, N^, is given by the 
integral Nt = J a{E~^)dLeff, where cr{E^) is the pho- 
tonuclear cross section. Assuming that the photonuclear 
cross section peaks at energy Eq and using the Thomas- 
Reiche-Kuhn (TRK) sum rule, Eq. (204 1, we can approx- 
imate this integral by Nt — x 6 x 10~^^^, where 
dL^ff/dE is expressed in units of cm^^s^^MeV^^. The 
giant resonances exhaust most part of the TRK sum rule 
and occur in nuclei at energies around £"0 = 15 McV. For 
1 GeV electrons dN{E^ = Eo)/dEj « 6 x 10~VMeV. 
With a luminosity of L^a = 10^^ cm~^s~^, one gets 
"^^^^^ w 6 X 10^^ cm~^MeV~^s~^ and a number of events 
Nr ^ 4:X 10~^NZ/A w 10-^ A s"^ Thus, for medium 
mass nuclei, one expects thousands of events per day. 

These estimates increase linearly with the accelerator 
luminosity, L^^a, and show that studies of giant reso- 
nances in neutron-rich nuclei is very promising at the 
proposed facilities [Haik05,Sud01]. Only a small frac- 
tion, of the order of 5%-10%, of the TRK sum-rule goes 
into the excitation of soft-dipole modes. However, these 
modes occur at a much lower energy, Er ~ 1 MeV, where 



the number of equivalent photons (see figure 50 1 is at 



least one order of magnitude larger than for giant reso- 
nance energies. Therefore, inelastic processes leading to 
the excitation of soft dipole modes will be as abundant 
as those for excitation of giant resonances. However, one 
has to keep in mind that it is not clear if experiments 
with very short-lived nuclei will be feasible at the pro- 
posed electron-ion colliders. 
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FIG. 50: Comparison between the virtual photon spectrum 
of 1 GeV electrons (dashed line), and the spectrum generated 
by a 1 GeV/nucleon heavy ion projectile (solid line) for the 
El multipolarity, as a function of the photon energy. The 
virtual photon spectrum for the ion has been divided by the 
square of its charge number. 



C. Non-inertial effects in Coulomb excitation 

Extremely large accelerations occur when atomic nu- 
clei collide. For instance, two lead nuclei in a head on col- 
lision with a center of mass kinetic energy of 500 MeV, 
reach a closest distance of 19.4 fm before they bounce 
back and move outward. At this distance each nucleus 
accelerates with an intriguing ^ 10^^ m/s^. Very few 
other physical situations in the Universe involve nuclei 
undergoing such large accelerations, usually related to as- 
trophysical objects, as in the vicinity of neutron stars and 
black holes, where huge gravitational fields exist. Here 
we will discuss the effects of large accelerations and large 
gravitational fields, and their possible infiuence on nu- 
clear reactions in the laboratory and in astrophysical en- 
vironments. Nuclear reactions are crucial for the forma- 
tion of stellar structures and their rates could be affected 
by various factors. 

There are two systems of reference which are often used 
to describe the effects of the collision: (a) the center-of- 
mass (cm) system of the two nuclei and (b) the system 
of reference of the excited nucleus. System (b) is appro- 
priate to use when the intrinsic properties of the excited 
nucleus is described in some nuclear model. A typical 
example is the case of Coulomb excitation. One assumes 
that the nuclei scatter and their cm wave functions are 
described by Coulomb waves due to the Coulomb repul- 
sion between the nuclei. Then one considers the residual 
effect of the Coulomb potential on the motion of the nu- 
cleons inside the nuclei. This is done by expanding the 
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Coulomb potential in multipoles and using the high order 
terms (higher than first order) as a source of the excita- 
tion process. In this approach one illustrates the privi- 
leged role of the cm of the nuclear system: the net effect 
of the external forces is to (i) accelerate all the particles 
together, along with the cm of the system, and (ii) to 
change the intrinsic quantum state of the system through 
the spatial variation of the interaction within the system. 
Thus the theoretical treatment of accelerated many-body 
systems is well under control in non-relativistic dynam- 
ics. 

In the non-relativistic case, the separation of variables 
into intrinsic motion and relative motion between the cm 
of each nucleus is a simple algebraic procedure. A prob- 
lem arises when one wants to extend the method to de- 
scribe intrinsic excitations of relativistic many-body sys- 
tems. Very few works exist in the literature addressing 
this problem. The reason is that for nuclear reactions in 
the laboratory, the effect is expected to be very small, 
a common belief which must be tested. Another other 
reason is that in stellar environments where the gravi- 
tational fields are large, huge pressures develop, "crush- 
ing" atoms, stripping them from their electrons, and ul- 
timately making nuclei dissolve into their constituents. 
Effects of nuclear excitation are not relevant in the pro- 
cess. 

Nuclei participating in nuclear reactions in a gaseous 
phase of a star follow inertial trajectories between colli- 
sions with other nuclei. Such trajectories are "free fall" 
trajectories in which all particles within the nucleus have 
the same acceleration. That is surely true in the non- 
relativistic case, but not in the relativistic one because 
retardation effects lead to corrections due to the nuclear 
sizes. The central problem here is the question regarding 
the definition of the center of mass of a relativistic many 
body system. We follow Refs. [LL89,Pry48,Few84], with 
few modifications, to show that a correction term pro- 
portional to the square of the acceleration appears in the 
frame of reference of the accelerated system. The results 
shown here are from Ref . [BIIK09] . 



1. Hamiltonian of an accelerated many-body system 

Starting with a Lagrangian of a free particle in an in- 
ertial frame and introducing a coordinate transformation 
into an accelerated frame with acceleration A, a "ficti- 
tious force" term appears in the Lagrangian when written 
in coordinates fixed to the accelerated frame. Thus, in 
an accelerated system the Lagrangian L for a free parti- 
cle can be augmented by a (non-relativistic) interaction 
term of the form ~mAz, that is 



1 



-mc 



-mv 



mAz, 



(336) 



In the relativistic case, the first step to obtain the 
Lagrangian of a many body system in an accelerated 
frame is to setup an appropriate measure of space-time 
in the accelerated frame, i.e. one needs to find out 
the proper space-time metric. The free-particle action 
S = —mc J ds requires that ds = (c — v'^ /2c + Az)dt, 
which can be used to obtain ds^ . To lowest order in 1 /c^ 
one gets 



ds' 



Az 



dt' ~ dx' - dy2 _ ^^2 ^ g^^d^'^d^t^, 

(337) 

where vdt — dr was used, with d^^^ = {cdt, dx, dy, dz) and 

= (500,-1,- 1,-1), 300= {l + Az/c^y. The indices 
/i run from to 3. Eq. ( |337[ ) gives a general form for the 
metric in an accelerated system. This approach can be 
found in standard textbooks (see, e.g. Ref. [LL89], § 87). 

From the definition for the Hamiltonian, iJ = p- v — L, 
with p = dL/dv = rnv/ goo — w^/c^, and using the ac- 



tion with the metric of Eq. ( 337 1 , after a straightforward 
algebra one finds 



H = 



goomc 



9m 



= c^/ gon (p^ + m^c^). 



(338) 



Expanding H in powers of 1/c^, one obtains 



H = 



2m 



1 - 



P 



^2 



-mAz 1 



P 



2m? c^ 



-0 



1 

(339) 

This Hamiltonian can be applied to describe a system 
of particles with respect to a system of reference moving 
with acceleration A, up to order \/c?. For an acceler- 
ated nucleus the obvious choice is the cm system of the 
nucleus. But then the term carrying the acceleration cor- 
rection averages out to zero in the center of mass, as one 
has (^j miAzi — 0). There is an additional small contri- 
bution of the acceleration due to the term proportional to 
p' . Instead of exploring the physics of this term, one has 
to account for one more correction as explained below. 

The above derivation of the Hamiltonian for particles 
in accelerated frames does not take into account that the 
definition of the cm of a collection of particles is also 
modified by relativity. This is not a simple task as might 
seem at first look. There is no consensus in the literature 
about the definition of the cm of a system of relativistic 
particles. The obvious reason is the role of simultaneity 
and retardation. Ref. [Pry48] examines several possi- 
bilities. For a system of particles it is found convenient 
to define the coordinates of the center of mass as the 
mean of coordinates of all particles weighted with their 
dynamical masses (energies). The relativistic (covariant) 
generalization of center of mass is such that the coordi- 
nates must satisfy the relation [Pry48] 



where z is the particle's coordinate along the direction of 
acceleration of the reference frame [LL89] . 



pO 



q^ = 



(340) 
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where the coordinates of the ith particle with respect 
to the center of mass are denoted by zf and the total 

momentum vector by — ^^pf • Ref. [Pry48] chooses 



eq. ( |340[ ) as the one that is most qualified to represent 
the definition of cm of a relativistic system, which also 
reduces to the non-rclativistic definition of the center of 
mass. We did not find a better discussion of this in the 
literature and we could also not find a better way to 
improve on this definition. 



The above definition, Eq. (3401, leads to the compact 
form, to order 



E 



500 



1 



2c2 



- 0, 



(341) 



where = (xi,yi,Zi) is the coordinate and Vi is the ve- 
locity of the ith particle with respect to the cm. 

For a system of non-interacting particles the condition 
in Eq. (341 1 implies that, along the direction of motion. 



rriiZi 



2c2 



T2~ 



(342) 



Hence, the Hamiltonian of Eq. (3391 for a collection of 
particles becomes 



H 



E 



2m, 



Pt 



A 



1 



2 ^m,zi+U{n)+0{-^), 
(343) 

where we have added a scalar potential U (r^), which 
would represent a (central) potential within an atom, a 
nucleus, or any other many-body system. 

Notice that the term proportional to —raAz com- 
pletely disappears from the Hamiltonian after the rel- 
ativistic treatment of the cm. This was also shown 
in Ref. [Few84]. It is important to realize that non- 
inertial efi^ects will also carry modifications on the inter- 
action between the particles. For example, if the particles 
are charged, there will be relativistic corrections (mag- 
netic interactions) which need to be added to the scalar 
i) = Y^.^iQiQilVi-^oV As shown in 



potential U (r^ 
Ref. [Few84], the full treatment of non-inertial effects 
together with relativistic corrections will introduce addi- 
tional terms prop ortional to A and A? in the Hamilto- 
nian of Eq. (3431, to order l/c^. Thus, a more detailed 



account of non-inertial corrections of a many-body sys- 
tem requires the inclusion of ^-corrections in the inter- 
action terms, too. We refer the reader to Ref. [Few84] 
where this is discussed in more details. Here we will only 
consider the consequences of the acceleration correction 



term in Eq. (3431 



^2 



E 



(344) 



2. Reactions in the laboratory 

For a nuclear fusion reaction, at the Coulomb radius 
(distance of closest approach. Re) the nuclear accelera- 
tion is given by 



Ac 



(345) 



where tuq = mjv^i^2/ (^i +^2) is the reduced mass 
of the system and mjv is the nucleon mass. For typical 
values, E = 1 MeV, Zi = Z2 = 10, and Ai = A2 = 
20, one obtains Rc = ZiZ2e'^/E = 144 fm and Ac = 
6.2 X 10^^ m/s^. This is the acceleration that the cm of 
each nucleus would have with respect to the laboratory 
system. 

The cm. of the excited nucleus is the natural choice for 
the reference frame. This is because it is easier to adopt 
a description of atomic and nuclear properties in the cm 
frame of reference. Instead, one could also chose the 
cm of the colliding particles. This later (inertial) system 
makes it harder to access the acceleration effects, as one 
would have to boost the wave functions to an accelerated 
system, after calculating it in the inertial frame. This is 
a more difficult task. Therefore we adopt the cm refer- 
ence frame of the excited nucleus, using the Hamiltonian 
of the previous section. This Hamiltonian was deduced 
for a constant acceleration. If the acceleration is time- 
dependent, the metric of Eq. ( |337 l also changes. Thus , 
in the best case scenario, the Hamiltonian of Eq. (3431 



can be justified in an adiabatic situation in which the rel- 
ative velocity between the many-body systems is much 
smaller than the velocity of their constituent particles 
with respect to their individual center of masses. 



3. Reactions involving halo nuclei 

The nuclear wave-function of a (s-wave) loosely-bound, 
or "halo" , state can be conveniently parameterized by 



a exp (— ar) 
2tt r 



(346) 



where the variable a is related to the nucleon separation 
energy through S — h^a^ /2mN ■ In first order perturba- 
tion theory the energy shift of a halo state will be given 

by 



N 
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85 



ZiZ2e^h 

R^TTLqC 



(347) 



Assuming a small separation energy 5 = 100 kcV, and 



using the same numbers in the paragraph after Eq. ( 345 1 , 
we get A£;^„ = 0.024 eV, which is very small, except 
for states very close to the nuclear threshold, i.e. for 
5^0. But the effect increases with Z"^ for symmetric 
systems (i.e. Zx — Z2= Ai/2). It is thus of the order of 



AE. 



N 



1 — 10 eV for larger nuclear systems. 
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There might exist situations where this effect could be 
present. For instance, the triple-alpha reaction which 
bridges the mass = 8 gap and forms carbon nuclei in 
stars relies on the lifetime of only 10""'^'' s of *Be nuclei. 
It is during this time that another alpha-particle meets 
^Be nuclei in stars leading to the formation of carbon 
nuclei. This lifetime corresponds to an energy width of 
only 5.57 ± 0.25 eV [Ti04]. As the third alpha parti- 
cle approaches ®Be, the effects of linear acceleration will 
be felt in the reference frame of *Be. This will likely 
broaden the width of the ®Be resonance (which peaks 
at Eh = 91.84 ± 0.04 KeV) and consequently its hfe- 
time. However, this line of thought could be wrong if 
one assumes that the third alpha particle interacts indi- 
vidually with each of the two alpha particles inside ^Be, 
and that the effects of acceleration internal to the ^Be 
nucleus arise from the different distances (and thus ac- 
celerations) between the third alpha and each of the first 
two. This effect has not been discussed elsewhere and 
deserves further investigation, if not for this particular 
reaction maybe for other reactions of astrophysical inter- 
est involving very shallow nuclear states. 



4- Nuclear transitions 

Many reactions of astrophysical interest are deduced 
from experimental data on nucleus-nucleus scattering. 
Important information on the position and widths of res- 
onances, spectroscopic factors, and numerous other quan- 
tities needed as an input for reaction network calculations 
in stellar modeling are obtained by the means of nuclear 
spectroscopy using nuclear collisions in the laboratory. 
During the collision the nuclei undergo huge accelera- 
tion, of the order of ~ 10^* m/s^. Hence, non-inertial 
effects will be definitely important. 

A simple proof of the statements above can be obtained 
by studying the Coulomb excitation. The simplest treat- 
ment that one can use in the problem is a semi-classical 
calculation. Following the same procedure leading to Eq. 
( 56 1 , we can calculate the contribution of the Hamilto- 



Applying Eq 
MeV, we find 



(p550| to the lead-lead collision at 500 

2 



(0.0093^)2. This yields very 

small results for the effect of non-inertial effects in single 
particle transitions {X ~ 1), but can become appreciable 
for the excitation of collective states such as the giant 
resonances, for which X ^ 1. 



XII. CODES 

GO SI A - A non-relativistic semiclassical coupled- 
channel Coulomb excitation least-squares search 
code, GOSIA, has been developed to analyze the 
large sets of experimental data. Manual and code 
available from: 



lhttp://www.pas. rochester.edu/~cline/Gosia/ 



nian of Eq. (344 1. In this case, A = Z^Z^e jni^r and 



the pcrtubativc potential = Vnin is given by 



DWEIKO - A relativistic coupled-channels code for 
Coulomb and nuclear excitations of El, E2, E3,M1, 
and M2 multipolarities. 

Manual: C.A. Bertulani, CM. Campbell, T. Glas- 
macher. Computer Physics Communications 152 
(2003) 317. 

Code available from the Computer Physics Com- 
munications Program Library: 

|http://www.cpc.cs.qub.ac.uk/cpc/| 

COULINT - A code which incorporates retardation 
corrections and recoil at intermediate energy col- 
lisions (~ 100 MeV/nucleon), appropriate for 
Coulomb excitation at rare isotopes facilities, such 
as GANIL, GSI, NSCL and RIKEN. Under de- 
velopment. Preliminary versions already available 
[Ber04]. 



XIII. APPENDICES 

Here I discuss a few of the basic principles of quantum 
mechanics which are necessary to understand this review. 



ZiZ2e^\ XrriN 
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(348) 



where we assume that X nucleons participates in the 
transition. One then finds 



mo / Ibis^hvoc^ 



(349) 



The ratio between the two transition probabilities (the 
above one and Eq. (56)) is 



/SXmNZiZi 



(350) 



A. Decay rates 

A quantum system, described by a wavefunction that is 
an Hamiltonian eigenfunction, is in a well defined energy 
state and, if it does not suffer external influences, it will 
remain indefinitely in that state. But this ideal situation 
does not prevail in excited nuclei, or in the ground state of 
an unstable nuclei. Interactions of several types can add 
a perturbation to the Hamiltonian and the pure energy 
eigenstates no more exist. In this situation a transition 
to a lower energy level of the same or of another nucleus 
can occur. 

An unstable state normally lives a long time compared 
to the fastest nuclear processes, e. g., the time spent by 
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a particle with velocity near that of the light to cross 
a nuclear diameter. In this way we can admit that a 
nuclear state is approximately stationary, and to write 
for its wave function: 



(351) 



|^'(r,t)p(iF is the probability to find the nucleus in the 
volume dV and, if the state described by ^' decays with 
a decay constant A, it is reasonable to write 



|vI/(r,t)p = |vI/(r,0)|V 



-At 



(352) 



To obey simultaneously (3511 and (3521, W must be 



a complex quantity with imaginary part —Xh/2. We can 
write: 



(353) 



which shows that the wave function ( 351 1 does not rep- 
resent a well defined stationary state since the exponen- 
tial contains a real part — Ai/2. However, we can write 
the exponent of (351 1 as a superposition of values corre- 



sponding to well defined energies E (for t>0): 



-{iEo/h+\/2)t 



(354) 



Functions connected by a Fourier transform relate to 
each other as 



lit) 



1 



+00 



\/27r 7-00 

\ /--l-oo 



(355) 
(356) 



V2^ J- 

This allows to establish the form of the amplitude A{E): 

AiE) = ^l^iE-Eo)/h^X/2]t (357) 

27rft Jo 

where the lower limit indicates that the stationary sate 



was created at the time t = 0. The integral (357) is of 
easy solution, giving 



A{E) = 



h\/2-2m{E ~ E^y 



(358) 



The probability of finding a value between E and E - 
dE in an energy measurement is given by the product 



A*iE)AiE) 



1 



(359) 



h^X^/A + ATT^{E~ Eoy 
and this function of energy has the form of a Lorentzian, 



with the aspect shown in figure 51 Its width at half- 
maximum is r = fi,A = H/t. The relationship 



tT : 



(360) 
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FIG. 51: 


Form of the distribution (3,5fM. 



between the half-life and the width of a state is directly 
connected to the uncertainty principle and shows that the 
longer time a state survives the greater is the precision 
with which its energy can be determined. In particular, 
only to stable states one can attribute a single value for 
the energy. 

The decay constant A was presented as the probability 
per unit time of occurrence of a transition between quan- 
tum states, and its values were supposed to be known 
experimental data. Now we will show that a formula to 
evaluate the decay constant can be obtained from the 
postulates of perturbation theory. 



B. Coupled-channels 

We can describe an unstable state by the addition of 
a perturbation to a stationary state. Formally we can 
write 



H = Ho + V, 



(361) 



where H is the Hamiltonian of the unstable state, com- 
posed by the non-perturbed Hamiltonian Hq and a small 
perturbation V. The Hamiltonian Hq satisfies an eigen- 
value equation 



Holpn ^ Enlpn, 



(362) 



whose eigenfunctions form a complete basis in which the 
total wavefunction , that obeys 



= ih 



can be expanded: 

* = ^a„(t)V;„e-^^"*/^ 



(363) 



(364) 



Using ( 361 ) and ( 364 1 in ( 363 1 , and with the aid of ( 362 1 



we obtain: 

ih J2 a„^«e-^^"*/'' = ^a„^„e-*^"*/^ (365) 
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with d„ = dan{t)/dt. Using the orthogonalization prop- 
erties of the tpn, let us multiply ( 365 1 by ^pl and integrate 
it in the coordinate space. From this, results the coupled- 
channels equations 



(366) 



where we introduced the matrix element (dr is the vol- 
ume element) 



(367) 



Vkn = / 'iPlVipn dr. 



C. Transition rates and cross sections 

Let us make the following assumptions about the per- 
turbation V: it begins to act at time t = 0, when the 
unperturbed system is described by an eigenstate ^pm- It 
stays at a very low value and after a short time interval 
becomes zero ai t — T. These assumptions allow us to 
say that the conditions 

a„i = 1, and a,„ = if (to ^ n) (368) 

are rigorously verified for t < Q and also works approxi- 
mately for t > 0. Thus, (366 1 has only one term and the 



value of the amplitude is obtained from 



ak 



Vurae'^'^--''-'^'/^ dt, 



(369) 



whose value must be necessarily small by the assumption 
that followed (368). The above approach is also known 
as first order perturbation theory. The integral of (369) 



gives 



ak 



Vkm [1 



exp [i 



Ek-E„ 



T)] 



Ek — Em 

We need now to interpret the meaning of the ampli- 
tude flfc. The quantity a^Ofc measures the probability of 
finding the system in the state k. This characterizes a 
transition occurring from the initial state to to the state 
fc, and the value of a^Ofe divided by the interval T should 
be a measure of the decay constant relative to the 
state k. The total decay constant is obtained by the sum 
over all states: 

— EKP 



A 



T 



(371) 



Let us now suppose that there is a large number of 
available states k. We can in this case replace the sum- 
mation in (371 1 by an integral. Defining p{E) as the den- 
sity of available states around the energy Ej^., we write 

\ak\^p{Ek)dEk 



A 



+ 00 



\Vkn 



2h 



or] 



{Ek - EmY 



p{Ek)dEk. 

(372) 



The function sm^x/x^ only has significant amplitude 
near the origin. In the case of (372), if we suppose that 
Vkm and p do not vary strongly in a small interval of the 
energy Ek around £",„, both these quantities can be taken 
outside of the integral, and we obtain the final expression 



'p{Ek) 



(373) 



Eq. ( 373 1 is known as the golden rule n- 2 (also known 



as Fermi golden rule), and allows to determine the decay 
constant if we know the wavefunctions of the initial and 
final states. 

We now describe how we can obtain a reaction cross 
section for the system A(a,b)B from Fermi's Golden rule. 
For the application of ( 373[ ) the entrance channel will be 
understood as an initial quantum state, constituted of 
particles of mass TOq hitting a target A of mass rriA- The 
final quantum state is the exit channel, where particles of 
mass TOfc move away from the nucleus B of mass to b ■ The 
essential hypothesis for the application of the golden rule 
is that the interaction responsible for the direct reaction 
can be understood as a perturbation among a group of 
interactions that describe the system as a whole. 

If we designate as A the transition rate of the initial 
quantum state to the final quantum state, we can write 
a relationship between A and the total cross section: 



A = = <y 



kah 



(374) 



where Va is now the velocity of the particles of the inci- 
dent beam. If dX is the part relative to the emission in 
the solid angle dH., then 



da 



-^dX 
kah 



(370) and, using (ISTSl) 



da 



rria 27r 
knh h 



\Vs^^dp{Ef), 



(375) 



(376) 



where the indexes i and / refer to the initial and final 
stages of the reaction. 

The density of final states can be easily obtained by 
counting the available number of states for the projectile 
in a large box within an energy interval dE. For a box 
with volume , one finds [BD04] 



P- 



dN 
dE 



prjihS 



(377) 



If we assume an isotropic distribution of momentum for 
the final states, the fraction dp that corresponds to the 
solid angle d^l is 



pnihS^ dO. kbrribS^ 
dp — - — ^TTT^ — — — — JrvrrUSi . 



2n^h^ 47r 



(378) 



Using (378) in (376), we obtain an expression for the 



differential cross section: 

da rriambki, 
dn ^ {2T:h'^Yka 



(379) 
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involving the matrix element 



I nnx^-^*ir^p)V^a^Axi+Hrc.)dT. (380) 



"if a, "^b, ^A, are the internal wavefunctions of the nu- 
clei a, b, A and B. xi^"*, X^~^ are the scattering wavefunc- 
tions of the relative momentum in the entrance channel 
a and in the exit channel /?. The volume was not writ- 



ten in ( 379 1 because it can be chosen as the unity, since 



the wavefunctions that appear in ( 380 1 are normalized where 



to 1 particle per unit volume. V is the perturbation po- 
tential that causes the "transition" from the entrance to 
the exit channel. It can be understood as an additional 
interaction to the average behavior of the potential and, 
in this sense, can be written as the difference between 
the total potential in the exit channel and the potential 
of the optical model in that same channel. 



D. Electromagnetic potentials 

The Lorentz equation, or the equation of motion of a 
charge in an electromagnetic field, is given by 



d mv V , 

— , =^g(E+- xB . 



(381) 



We now introduce the electromagnetic potentials (</>, A) 
by means of the relations 



B = V X A 



E 



-V0- 



1 dA 



(382) 



In terms of these potentials, we get 



dt ^1 - v2/c2 



= q 
= V 



-V0-if +-.(V X A) 
c ot c 



c 



q dA 

c~dt 



where we have used the convective derivative 
dA[r{t),t] _ ( d 



dt 



dt 



v.V A[r(t),i] 



(383) 



(384) 



V 



V • A 



We can rewrite ((383)) as 
d I rnv q ^ 

(385) 

This equation has the form of the Euler-Lagrange 
equation 



dt 



(VrC) - VrC = 0. 



The suitable Lagrangian that reproduces Eq. (3831 
when replaced in Eq. ( 385 ) is given by 



£(r,v) = — rnc^\/l — v^/c^ -\ — v • A(r, <) — g(/)(r, i). 

^ (386) 



The canonical momentum is 
mv q 



Vv/: 



Vl - v2/c2 



A(r,i) 



TOV/ V 1 — v2/c2 



A(r,t) 
(387) 



(388) 



is the kinetic momentum and qA/c is the momentum 
carried by the electromagnetic field. 
The Hamiltonian is 



mc 



Vl -v2/c2 



We can rewrite (|389|) as 
n{r,p) = c 



q 



p--A(r,t) +{mcy 
c 



1/2 



For non-relativistic particles, 

I Ti I I q 



and 



Ti (r,p) = mc 



- A I << mc 
c 



2 , (P- 9A/c)2 



2m 



(389) 

q(j){r,t). 
(390) 

(391) 
(392) 



The second term has as part of its contribution the 
quantity 

{qAf I 2mc^ 

which is relevant only in processes where two photons are 
involved and may be ignored. The remaining terms yield 
the Hamiltonian for the electromagnetic interaction. 



v.A 



(393) 



where the rest mass and kinetic energy of the particle 
was subtracted from Eq. ( 392 1 . 



For systems involving a charge density p (r, t) and cur- 
rent density j(r,t), Tiint can be generalized to 



ij.A 

c 



(394) 



E. 7-ray emission following excitation 

After the excitation, the nuclear state |//) can decay 
by gamma emission to another state \Ig)- Complications 
arise from the fact that the nuclear levels are not only 
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populated by Coulomb excitation, but also by conversion 
and 7-transitions cascading down from higher states (see 
figure [7|a)). To compute the angular distributions one 
must know the parameters A/ (i — j) and e/ (i — > j) 
for I > 1 [AW66], 

j) = j) A? ^ j) , (395) 

where ai is the total /-pole conversion coefficient, and 

1/2 



?7r(/+ 1) 1 /tj\2;+ 



l[{2l + iy.lf ti ^c. 
X (2/j + 1)"^/^ i'^(')X(7rZ) li) ,(396) 

with s(/) = / for electric (tt = E) and s(Z) = Z + 1 for 
magnetic (vr = M) transitions. The square of A^; is the 
?-pole 7-transition rate (in sec~^). 

The angular distributions of gamma rays following the 
excitation depend on the frame of reference used. In our 
notation, the z-axis corresponds to the beam axis, and 
the statistical tensors are given by (we use the notation 
of [WB65,AW66]) 



"k^yj) (2/1 + 1) 



E (-1) 



If+Mf 



Mf=-{M'j + K),M'j 

,M' !)E«W(^^i)«W/(^i). 

(397) 



-M 



Ah 



where / is the state from which the gamma ray is emitted, 
and 1 denotes the initial state of the nucleus, before the 
excitation. To calculate the angular distributions of the 
gamma rays one needs the statistical tensors for k — 
0,2,4 and -k < n < k (see [WB65,AW66]). 

Instead of the diagram of figure |7][a), we will con- 
sider here the much simpler situation in which the 7- 
ray is emitted directly from the final excited state / to 
a lower state g, which is observed experimentally (see 
figure [7|b)). The probability amplitude for this process 
is 

a^^f^g = E "^-/ (^s^flka \H^\ IjMf) , (398) 

Ml 

where {IgMgka \H^ \ IfMf) is the matrix element for the 
transition f g due to the emission of a photon with 
momentum k and polarization a. The operator Hj ac- 
counts for this transition. The angular dependence of the 
7-rays is given explicitly by the spherical coordinates 9 
and (p of the vector k. 

Since the angular emission probability will be normal- 
ized to unity, we can drop constant factors and write it 



as (an average over initial spins is included) 



w{e) = 



Mi, M„ 



E 

Mi, M„ 



E° 

M, 



I„M„-ka\H^\IfMf, 



The transition operator can be written as 



{nuc) 
Im 



(7) 

Im 1 



(399) 



(400) 



Ira 



Im 



where the first operator in the sum acts between nuclear 
states, whereas the second operator acts between photon 
states of well defined angular momentum, I, m. 

Expanding the photon state jkcr) in a complete set 
\lm) of the photon angular momentum, and using the 
Wigner-Eckart theorem (angular momentum notation of 
Ref . [Ed60] , one gets 



(J.Mgka \H^\ IfMf) = (l^^l'"^) {W 



ffilrn) 

7 



IfMf 



I In 



f - ^g 
Mf m Mg 



ka|?TO>//J|i/W||/ 



(401) 



One can rewrite |kcr) in terms of |z(t), i.e., in terms 
of a photon propagating in the z-direction. This is ac- 
complished by rotating jkiT) to the z-axis, using of the 
rotation matrix [BS94], f^^/, i.e., 

(ka|;m) = D'rnm' ^ k) {za\lm') . (402) 

m' 

Expanding the photon field in terms of angular mo- 
mentum eigenfunctions, one can show that [FS65,EG88] 



(zallrmr) 



21+1 



2 ^<Tm 



for TT = E 
for TT ~ M. 



One has now to express the operator ofZl/ in Eq. (1400 1 



(403) 



in terms of the electric and magnetic multipole parts of 
the photon field. This problem is tedious but straight- 



forward [BR53]. Inserting eqs. (402 1 and (403 1 in Eq 



(4011, yields (neglecting constant factors) 

{IgMgka \H,\IfMf) = Y {-lY'-'^''' 

xy'(2/+l)(2/, + l)(^^" ^ ^ 

X D^^„ (z — > k) [Ab; -I- al^Mi] 
where A^; is given by Eq. p96| . 



'9 " ^/ 

Mf m -Mf 



(404) 
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Inserting Eq. (4041 into Eq. (399) one gets a series 
of sums over the intermediate values of the angular mo- 
menta 



h I If 



xy'(2;+l)(2Z' + l)(2/, + l)(^^^ ^ 

/ ..M,+M'f-l~l' ( Ig I' If \ 

^'^> [Mg m' -M'j ) 

X Dl^\Di,V AiA;, , (405) 



where A; = A^;; + I Ami- The product A;A^, is always 
real since (— l)'*^') = 11 (the parity). 

Assuming that the particles are detected symmetri- 
cally around the z-axis one can integrate over (pparticie, 
what is equivalent to integrating, or averaging, over cf)^. 
This yields the following integral 



d(j)D. 



I 

ma 



2J + 1 ^ / J I' \ f I J I' 



\^ TO — ?71 J \ (T — (T 



Pj (cos ( 



(406) 



To simphfy further Eq. ( 405 1 we use (see Ref. [AW75] , 
p. 441, Eq. II.A.61) 



^ , /, / If \[ I9 I' If 

Mg m ~Mf ) \ Mg to' -M'f 



M. 



= (-l)"''"^' E (-1)''+"'"*^^ {2k + 1) 

(If I'f k \ ( I I' k 
^ {Mf -M} KjXl'f If Ig 



I V k 

TO —to' K 



(407) 



and 



=(-1,1) 



2 ( 1 -1 j ^'^'^^i' ' ^ ^^^^ 



(408) 



, for j = odd, 
where use has been made of the parity selection rule 



I (~1) 5 for electric transitions 
^ I ( — 1)'^^ , for magnetic transitions . 



Eq. ( |405| becomes 
WiO)= J2 (-1)'"'+'^+*'^ 



Qi ^/flj ^f, 



,k,K.,M f ,M 



1,1' ,m,m' 



X J {21 + 1) {21' + 1) {21 f + 1) (2j + 1) (2fc + 1) 



fc 



/f V k 



Mf -M'j K ) \ Mf -to' k 

z J r \ / z J /' 
1 -1 y TO -to' 

X A,A;',Pj (cos 6*) . 





I' 






If 





(409) 



Using 



TO —to' 



? j I' 
1 -1 



(410) 



one gets 

W{e)^ J2 (_l)W+'= ^{21 + 1) (2^' + 1) {2If + 1) 

fc = cvcn,Mj, 

x(2fc+i)k_/(i;^ -^/, 





I' 






If 





or in a more compact form 
W{e) ^ ^ {-lf'^\a,^f\'Fu{l^',Ig,If) 

fc = cvcn,Mj,A/_^, 



if, -M^ ) ^/2fcTlF, (cos 0) A, a;, , 

(411) 



where 

Fu{l.l'jgJf) = {-lY^-'^-' 



XyJ{2l + 1) {21' + 1) {21 f + 1) (2fc + 1) 
I I' k \ ( I I' k 



1 -1 



If If 1 9 



(412) 



The angular distribution of 7-rays described above is 
in the reference frame of the excited nucleus. To obtain 
the distribution in the laboratory one has to perform the 
transformation 



9l = arctan 



sin 



7 [cos 9 + P] 



(413) 



and 



W{eL) {1 + [3cos9fw{e) , (414) 

where 7 = {\ — [3'^)^'^/'^ and (3 = v/c, where v is the 
nucleus-nucleus relative velocity. The photon energy in 
the laboratory is Ef - = jEP^ (1 + Pcos0). 
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F. Second-order perturbation theory 

To second-order, the amplitude for a two-step excita- 
tion to a state |2 > via intermediate states |1 > is given 

by 



nd _ \ " 1 



dt e 



I iUigt' 



dt e"^"' V2iit) 



vMt') , 



(415) 



where V2i{t) is a short notation for the interaction po- 
tential inside brackets of the integral of Eq. ((415 1 ) for 
the transition |1 >^ |2 >. 

Using the integral representation of the step function 

1 f°° e^*''(*^*') 



one finds 

2nd _ 
"20 — 



- lim 



i6 



-dq 



1, if <> t' 
0, if < < t', 
(416) 



I X!"2f(^2i) a}S*(wio) 



27r 



-oo g 



a^f{uj2i-q) al'o*{uJio + q) 



where V stands for the principal value of the integral. For 
numerical evaluation it is more appropriate to rewrite the 



principal value integral in Eq. (4171 as 

r 



-oo 9 
' dq 

q 



a^f (W2i -g) a}g*(wio + (?) 



(^^21 + g) aig*(wio -(?) 



(418) 



To calculate a^**(a;) for negative values of cu, we note 
that the interaction potential can be written as a sum of 
an even and an odd part. This implies that a^''*(— cj) = 

~'a^'\uj)\*. 



G. Coulomb excitation of a harmonic oscillator 

Let us consider the Coulomb excitation of a linear har- 
monic oscillator. For dipole excitations the interaction 
Hamiltonian has the form 



ntm^xfit) + zg{t) 



(419) 



where f(t) and g(t) are time-dependent functions and x 
and z are the intrinsic nuclear coordinates perpendicular 



and parallel to the beam, respectively. This field induces 
independent oscillations in the perpendicular and parallel 
directions. 

Introducing the occupation numbers — 0, 1, 2, • • • 
and Hz — 0, 1, 2, • • • the total number of occupied 
states is 



N ^ rix + riz 



(420) 



Considering only the oscillations along one of the di- 
rections, e.g., along the x-direction, the full Hamiltonian 
is 



n = no~ xf{t) 



(421) 



where the minus sign is arbitrary and is introduced to 
simplify the equations. 

The solutions for the intrinsic Hamiltonian, Hq = 
mw^x^/2, are obtained in terms of the Hermite poly- 
nomials (here, n — ) 

/mcL;\i/4 1 / fmuj \ ( mujx 

(422) 

Since this forms a complete basis, the wavefunction at 
time t can be written as 



(417) ^"(2;, ^) = E [~*^ -'^Z^) V'n(a;) (423) 



n=0 



with the condition a„(0) = J„^o • Inserting this expansion 
into the time-dependent Schrodinger equation 



, 9 , , , 



^'5' 1 22 
2m ox'^ 2 



i){x,t) 

(424) 

and using the identity (which can be derived from 
Hn+i{x) ^ 2xHn{x) - 2nHn-i{x) ) 



Xlpnix) 



f ft 



V 2mLjj 



1/2 



[\/?i'0n-i(a;) + Vn+ li/;„+i(x)] 

(425) 



we get 



dar. 



ih exp [—iuj {n + 1 /2) t] —^(t)ipn{x) — 



n=0 



1/2 



/(t) VexpHw(n+l/2)t] a„(i) 

X [Vr^V'„-i(a;) + Vn+ lV'„+i(a;)] ■ (426) 
Using the orthogonality conditions of ipn one obtains 



where 



(t) = z/o(t) [e-*'^*\Ama„+i(t) +e-^'^*V^a„_i(t)] , 

(427) 



foit) = {2mhLu)-'/^ fit). 



(428) 



65 



The last sum yields , where 

\m'^\x.{t)f + \xAt)f 

Thus, one gets the Poisson formula 



nh(3 




P = — exp[- O] 



P = O exp[- O] 



FIG. 52: The harmonic oscillator model for the excitation of 
collective states in nuclei. 



The above equation allows to obtain a„(i) by iteration, 
starting with n — 0. The solution is 



a„(t) = i"^^exp 



where 



But, 
t 







Jo 



(429) 



(430) 







dt'Ut')e-^-'x{t')= I dt'^{t')x{t')=^\x{t)\' 

(431) 



and the probability to excite the nth state is 



P„(^)=^exp 



(432) 



For a nucleus, n can also be considered as the number 
of phonons for the vibration along the x-direction. The 
probability to have vibrations along the z-direction 
and rij. vibrations along the x-direction is 



^(,)Jx.(t)|--|x.(t)|--^^^ 



Ix.(i)l'-Ix.(t)l'_ 
(433) 



The probability that N-phonons will be excited is 

N 



Pn — ^ Pn^.n^ — ^ Pn^,N-n^ 

^ \Xz{i)\ \xi{t)\ 



n^=0 

X exp 
1 



n,\{N~n,)\ 

[-\x.{t)\'-\x.{t)\\ 



exp 



-\xMV'\xAt)Y 



n,\{N-n,)\ 



(434) 



P^ = ^ I0(t)l'^ exp 



(435) 



(436) 



for dipole excitations of a harmonic oscillator. In Refs. 
[AW66,AW75] it was shown that the same result can be 
obtained for quadrupole excitations. Thus, the above 



result is very general (see Figure 52 1 



The Poisson equation is very useful. It shows that all 
one needs to obtain the probability to excite the nth state 
is to calculate the probability to excite the = 1 state 
with first-order perturbation, -P/j^j** = 10(01'^ i ^iid use 
it in the Poisson formula. Although this only holds for 
a harmonic oscillator system, it has been largely used in 
many theoretical considerations of more complicated sys- 



tems. The factor exp 



accounts for the loss of 



occupation probability due to the excitation to all states 
from a given one. It also is responsible for conservation 
of unitarity, so that X]Af=o P^ ~ ^ ■ 



H. Schrodinger equation in a space-time lattice 



The operations [] Uq, in the r.h.s. of Eq. (2841 is easy. 
The operation [] ^ u'^ is more complicated. The problem 
is to find the vector u in the equation 



where 



A u 



(437) 



u is a vector composed with the u^^^ = u{rj,t) compo- 
nents of the wave-function u{r, t) . In Eq. 284 ^ is a 



tri-diagonal operator (matrix). In matrix notation 



/ ^1 \ 

V2 



(A, 
A. 





Ao Ai 



2 ^2 ^2 • 

A^ A3 A+ 



V 



An I 



(438) 



V «JV / 

This involves the following relations 

A':[ui^i+AiUi+Afui+i= Vi 
Assuming a solution of the form 

Mj+i = ttiUj -I- Pi (440) 



(439) 



and inserting in ( 439 1 we find the recursion relations for 
at and pi : 



7, Ai , = 7, {A+(3, - V,) (441) 
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where 



A,. 



(442) 



At the end of the lattice we assume un — vn . This 
implies that a^^^i — and Pn-i — ^jv • 



For the problem defined by Eq. (284|) we have 

A] 



A,= ^ + 2+|^e\ and 



A+ = -l. (443) 



We can now determine at and Pi by running Eqs. 
( 441 1 backwards from i = N — 2 down to i = 1 . Then 
we use Eq. ( 440 1 running forward from i — 2 to N , 



assuming that ui at the other extreme of the lattice is 
given by m = Vi/Ai . 



Another way to solve the problem (438) for u is by a 
LU-decomposition, followed by a forward and backward 
substitution. This method does not need to involve the 
Dirichlet condition. 



Let us assume that the A matrix in ( 438 1 can be writ 



ten as a product of B and U matrix, where 



A = LU = 



/ bi ci 

02 62 C2 

as 63 C3 

V 



(444) 





1 










0:2 


1 










as 


1 












V 









L 



Then, AU = V - 
elements of y are 



L{Uu) — V , or Ly — V 



(445) 
The 



/ 2/1 \ 

2/2 



KVN J 



/3i 71 

/32 72 

/32 73 

0.. 



/ wi \ 



(446) 



This can be solved by backwards substitution 

VN 



UN 



PiUi + 7iUj+i = yi 



or 



Ui = {yi - iiUi + 1) / 
The other matrix equation 



(447) 
(448) 



/I 0. 

ai 1 0. 

ag 1 . 

04 1 



(yi \ 

2/2 



\ 2/JV / 



V2 



\un / 



(449) 



can be solved by forward substitution 

yi=vi, aiyi-i + yi = Vi , (450) 



yi= Vi-Uiyi^i, 



((D.15)) 



Now, we need to find Ui and Pi as a function of the 
original elements of A 



/I ... \ 
a2 1 . . 
ag 1 . . 



//3i 71 
/32 72 
/33 73 



/ V 



/ 5i ci . . \ 
62 62 C2 . . 
ag 63 C3 . . 



Pn ) 



■ ■\ 




(Pi 




ai 





. 




Also, 






P2 


72 


. 






u ^ 








P3 


73 ■ 






'■ J 












) 


or 



V . . 

which implies 

bi^Pi: &2 = a27i+/32, 

or 

Pi= b,-'y^_j_ai , 



, bi = 7i_iai + Pi 



(451) 



P^- 



(452) 



Thus, knowing a^^ bi and Ci , one can go upwards with 
this set of equations to solve the problem. 
In our particular case, 

a, - c, = -1. (453) 

Thus, the above equations simplify to 

1 



P 



bi, 



2/1 = vi, 



P^= h- 

y%= Vi+ 



UN 



VN 

Pn'' 



Ui = 



Pr^l 
{y^ + U,+i) 

P^ 



i = 2,...,N 
i^2,...,N 



i = A^-l,...,l. 



(454) 
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